Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,002,530 G→A 100% G729D (GGT→GAT)  xdhA → xanthine dehydrogenase, molybdenum binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,002,5300GA98.8% 218.0 / ‑5.9 64G729D (GGT→GAT) xdhAxanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (36/26);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGCCGATGCCATCAACCTGCTG  >  NC_000913/3002385‑3002673
                                                                                                                                                 |                                                                                                                                               
ccGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGatgt                                                                                                                                              >  8:7187/1‑149 (MQ=255)
      aTGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaat                                                                                                                                        <  2:240264/149‑1 (MQ=255)
      aTGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaat                                                                                                                                        <  7:345101/149‑1 (MQ=255)
      aTGCCGGATCTGCCACAACTGGAAAGCGCGTGCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaat                                                                                                                                        <  5:202747/149‑1 (MQ=255)
          cGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGcaatcaa                                                                                                                                     >  1:182770/1‑148 (MQ=255)
                  ccACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGcc                                                                                                                            >  6:24047/1‑149 (MQ=255)
                  ccACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGcc                                                                                                                            >  7:101389/1‑149 (MQ=255)
                    acaACTGGAATGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGCGCCACCAATAATTCCTGTTGCCGCTCCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCACTCAATACACTGCCGc                                                                                                                          >  1:361418/1‑149 (MQ=255)
                    acaACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGc                                                                                                                          >  8:84148/1‑149 (MQ=255)
                    acaACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGc                                                                                                                          >  6:301441/1‑149 (MQ=255)
                      aaCTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCt                                                                                                                         >  1:17547/1‑148 (MQ=255)
                         tGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACg                                                                                                                     >  2:130020/1‑149 (MQ=255)
                                    ttCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGtt                                                                                                          <  1:395964/149‑1 (MQ=255)
                                    ttCGTCGAAATAAATGAGCCGCAATCGGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGtt                                                                                                          <  5:380088/149‑1 (MQ=255)
                                       gtcgAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTtata                                                                                                       >  7:311941/1‑149 (MQ=255)
                                       gtcgAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTtata                                                                                                       <  1:10688/149‑1 (MQ=255)
                                          gAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGCTGTTGCAATCAATACACTGCCGCTGCCGCCAAAACGGTTATATg                                                                                                     >  6:207937/1‑148 (MQ=255)
                                           aaaTCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTGGCCATCATTACCCTGCCGCGGACGCAAACACGGTTATATg                                                                                                     >  6:111757/1‑147 (MQ=255)
                                           aaaTCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGaa                                                                                                   <  7:84148/149‑1 (MQ=255)
                                             aTCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAaga                                                                                                 <  1:130020/149‑1 (MQ=255)
                                                aaTGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGtt                                                                                              >  3:232832/1‑149 (MQ=255)
                                                aaTGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGtt                                                                                              >  8:366486/1‑149 (MQ=255)
                                                aaTGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGtt                                                                                              >  5:82965/1‑149 (MQ=255)
                                                aaTGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTACTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGtt                                                                                              >  6:219168/1‑149 (MQ=255)
                                                  tGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTGCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTcc                                                                                            <  6:292900/149‑1 (MQ=255)
                                                  tGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTcc                                                                                            <  2:17547/149‑1 (MQ=255)
                                                  tGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTcc                                                                                            <  2:182770/149‑1 (MQ=255)
                                                        gCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTgg                                                                                      >  7:356359/1‑149 (MQ=255)
                                                             ccGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGga                                                                                 >  2:396514/1‑149 (MQ=255)
                                                                  tACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGGGGTGAAGATGTCTACCTATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGgattgat                                                                            <  5:301441/149‑1 (MQ=255)
                                                                  tACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACCCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGCCGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGgattgat                                                                            >  7:189252/1‑149 (MQ=255)
                                                                      gACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGa                                                                        <  8:356359/149‑1 (MQ=255)
                                                                      gACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGa                                                                        <  8:311941/149‑1 (MQ=255)
                                                                       aCATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAg                                                                       >  5:329646/1‑149 (MQ=255)
                                                                          tAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGAt                                                                    <  8:215077/149‑1 (MQ=255)
                                                                              tCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAAc                                                                 >  7:298764/1‑148 (MQ=255)
                                                                                  tGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAAcatcat                                                            >  1:118321/1‑149 (MQ=255)
                                                                                          ccaccaATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGAtt                                                    <  8:101389/149‑1 (MQ=255)
                                                                                           caccaATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGAttt                                                   >  7:23862/1‑149 (MQ=255)
                                                                                             ccaATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGAtttt                                                  >  2:316354/1‑148 (MQ=255)
                                                                                             ccaATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGAtttt                                                  <  1:316354/148‑1 (MQ=255)
                                                                                              caATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGc                                                <  6:329646/149‑1 (MQ=255)
                                                                                              caATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGc                                                >  8:133050/1‑149 (MQ=255)
                                                                                               aataatTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGaa                                                                                                   <  8:308564/97‑1 (MQ=255)
                                                                                               aataatTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGaa                                                                                                   >  7:308564/1‑97 (MQ=255)
                                                                                                 taatTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGCTGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGcttc                                             <  2:243156/149‑1 (MQ=255)
                                                                                                   atTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGcttctt                                           <  7:7187/149‑1 (MQ=255)
                                                                                                        tGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGcttctt                                           <  1:397682/144‑1 (MQ=255)
                                                                                                        tGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGcttctt                                           >  2:397682/1‑144 (MQ=255)
                                                                                                         gTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATc                                     >  4:216273/1‑149 (MQ=255)
                                                                                                            gCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATcgc                                   <  2:118321/148‑1 (MQ=255)
                                                                                                                    tATTCGTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGAGTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAAcc                          >  5:85628/1‑149 (MQ=255)
                                                                                                                         gTAAGGCGGTGAAGATGGCTACCGATGTTGCAATCAATACTCTGCCGCTGTCGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGc                     <  2:361418/149‑1 (MQ=255)
                                                                                                                         gTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGc                     >  1:267134/1‑149 (MQ=255)
                                                                                                                         gTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGc                     >  7:198856/1‑149 (MQ=255)
                                                                                                                         gTAACGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGc                     >  6:179064/1‑149 (MQ=255)
                                                                                                                           aaCGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCAGCCGCCAAAACGGTTTTATGAAGAGTTCCAGCTGGCAGGATTGATTTGAGGATAACATCATGTGTGAATTTGCTTCATACCATCGCGCAGCAACCCTAGCCg                   >  8:157395/1‑149 (MQ=255)
                                                                                                                            aCGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGCCGa                  >  7:59624/1‑149 (MQ=255)
                                                                                                                            aCGCGGTGAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGCCGa                  >  6:283665/1‑149 (MQ=255)
                                                                                                                                   gAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGaa                                                                                                   <  4:35472/61‑1 (MQ=255)
                                                                                                                                   gAAGATGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATAAGaa                                                                                                   >  3:35472/1‑61 (MQ=255)
                                                                                                                                       aTGGCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCCACCCTTGCCGATGCCATCAAcc       >  6:355456/1‑149 (MQ=255)
                                                                                                                                         ggCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTg                      >  1:322391/1‑132 (MQ=255)
                                                                                                                                         ggCTACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTg                      <  2:322391/132‑1 (MQ=255)
                                                                                                                                            tACCGATGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGCCGATGCCATCAACctgctg  >  6:78529/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGACGCCAAAACGGTTATATGAAGAGTTCCATCTGGCAGGATTGATTTGAGGATAACATCATGTTTGATTTTGCTTCTTACCATCGCGCAGCAACCCTTGCCGATGCCATCAACCTGCTG  >  NC_000913/3002385‑3002673

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑