Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,207,552 G→T 100% M172I (ATG→ATT aceA → isocitrate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,207,5520GT97.4% 725.7 / 9.0 228M172I (ATG→ATTaceAisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base T (109/113);  total (110/118)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 229 total aligned reads displayed.

ATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCT  >  CP009273/4207406‑4207696
                                                                                                                                                  |                                                                                                                                                
aTCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTaa                                                                                                                                                >  2:1955740/1‑149 (MQ=255)
 tCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATGACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTaaa                                                                                                                                               <  2:2296928/149‑1 (MQ=255)
  caacaaCACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAg                                                                                                                                              >  1:3333423/1‑149 (MQ=255)
   aacaacACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGc                                                                                                                                             >  2:13112/1‑149 (MQ=255)
   aacaacACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGGTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGc                                                                                                                                             >  1:3716206/1‑149 (MQ=255)
    acaacaCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCg                                                                                                                                            >  1:641764/1‑149 (MQ=255)
    acaacaCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCg                                                                                                                                            >  2:3446358/1‑149 (MQ=255)
     caacaCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCga                                                                                                                                           >  1:2477314/1‑149 (MQ=255)
         acCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCgatgat                                                                                                                                       <  2:3209074/149‑1 (MQ=255)
          ccTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGAtt                                                                                                                                      <  1:1087023/149‑1 (MQ=255)
           cTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTg                                                                                                                                     >  1:3116619/1‑149 (MQ=255)
             tCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGaa                                                                                                                                   <  2:4019354/149‑1 (MQ=255)
              ccGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAg                                                                                                                                  <  1:4037281/149‑1 (MQ=255)
                gtcgtGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGcc                                                                                                                                >  2:2570972/1‑149 (MQ=255)
                   gtGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGt                                                                                                                             >  1:1066602/1‑149 (MQ=255)
                   gtGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGt                                                                                                                             >  2:1587481/1‑149 (MQ=255)
                        gATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGc                                                                                                                        >  1:1939372/1‑149 (MQ=255)
                         aTCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCg                                                                                                                       >  2:180757/1‑149 (MQ=255)
                           cAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGc                                                                                                                     <  2:3716206/149‑1 (MQ=255)
                            aGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCa                                                                                                                    >  2:1753625/1‑149 (MQ=255)
                              aTCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGt                                                                                                                  <  1:1151753/149‑1 (MQ=255)
                                 cAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCa                                                                                                               <  1:1997742/149‑1 (MQ=255)
                                  aaTGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCAc                                                                                                              <  1:512745/149‑1 (MQ=255)
                                   aTGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACt                                                                                                             <  1:1770892/149‑1 (MQ=255)
                                   aTGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACt                                                                                                             >  1:2092134/1‑149 (MQ=255)
                                    tGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACtt                                                                                                            >  2:3028299/1‑149 (MQ=255)
                                    tGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGAGTACTTCCTGCCGAGCGTTGCCGATGCGGAAGCCGGGGGTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGAGGAGTGAAGCCGGTGCAGCGGCAGTTCAcgt                                                                                                            >  2:3296289/1‑147 (MQ=255)
                                     ggTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTc                                                                                                           >  1:4098497/1‑149 (MQ=255)
                                      gTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCg                                                                                                          <  2:647299/149‑1 (MQ=255)
                                       tCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGa                                                                                                         <  1:2758503/149‑1 (MQ=255)
                                        ccGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGaa                                                                                                        >  2:1782393/1‑149 (MQ=255)
                                          gcgGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGa                                                                                                      >  2:822145/1‑149 (MQ=255)
                                             ggCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCa                                                                                                   >  2:796219/1‑149 (MQ=255)
                                              gCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAg                                                                                                  >  2:2131527/1‑149 (MQ=255)
                                               cATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGc                                                                                                 <  1:2177407/149‑1 (MQ=255)
                                               cATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGc                                                                                                 >  2:1480314/1‑149 (MQ=255)
                                                aTTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCt                                                                                                <  1:791076/149‑1 (MQ=255)
                                                aTTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCt                                                                                                <  1:3601706/149‑1 (MQ=255)
                                                 ttGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTg                                                                                               <  2:1968226/149‑1 (MQ=255)
                                                 ttGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTg                                                                                               <  2:1967024/149‑1 (MQ=255)
                                                     gCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGt                                                                                           <  1:3466411/149‑1 (MQ=255)
                                                     gCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGt                                                                                           >  2:2236359/1‑149 (MQ=255)
                                                         ggCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGt                                                                                       >  1:3843412/1‑149 (MQ=255)
                                                         ggCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGt                                                                                       >  1:2662831/1‑149 (MQ=255)
                                                            gATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaa                                                                                    <  2:408437/149‑1 (MQ=255)
                                                              tCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaaga                                                                                  >  1:3448173/1‑149 (MQ=255)
                                                                 gcgcTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAAt                                                                               >  2:661893/1‑149 (MQ=255)
                                                                   gcTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGc                                                                             >  2:3443770/1‑149 (MQ=255)
                                                                    cTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCg                                                                            <  2:30800/149‑1 (MQ=255)
                                                                     tATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCgg                                                                           <  1:2379800/149‑1 (MQ=255)
                                                                      aTGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGt                                                                          >  1:2723987/1‑149 (MQ=255)
                                                                       tGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTc                                                                         <  1:1855541/149‑1 (MQ=255)
                                                                          cGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTcaca                                                                      <  2:1866427/149‑1 (MQ=255)
                                                                            aTTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATg                                                                    <  1:1962696/149‑1 (MQ=255)
                                                                            aTTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATg                                                                    <  1:2682253/149‑1 (MQ=255)
                                                                             ttACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATgg                                                                   <  2:297596/149‑1 (MQ=255)
                                                                             ttACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATgg                                                                   <  1:455861/149‑1 (MQ=255)
                                                                              tACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATggg                                                                  <  2:3911380/149‑1 (MQ=255)
                                                                              tACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATggg                                                                  >  2:2089032/1‑149 (MQ=255)
                                                                                     tGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCaaa                                                           >  2:2408491/1‑149 (MQ=255)
                                                                                     tGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCaaa                                                           >  1:3352924/1‑149 (MQ=255)
                                                                                      gCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAg                                                          >  2:3738112/1‑149 (MQ=255)
                                                                                       ccGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGt                                                         <  2:1449521/149‑1 (MQ=255)
                                                                                        cGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtt                                                        <  2:3645608/149‑1 (MQ=255)
                                                                                        cGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtt                                                        >  2:3184712/1‑149 (MQ=255)
                                                                                          aTCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtttt                                                      >  1:2913801/1‑149 (MQ=255)
                                                                                          aTCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtttt                                                      <  2:3397814/149‑1 (MQ=255)
                                                                                          aTCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtttt                                                      >  2:3270131/1‑149 (MQ=255)
                                                                                          aTCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtttt                                                      >  1:2529082/1‑149 (MQ=255)
                                                                                          aTCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtttt                                                      <  2:1049898/149‑1 (MQ=255)
                                                                                          aTCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtttt                                                      >  1:539916/1‑149 (MQ=255)
                                                                                            cGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAg                                                    <  1:3801845/149‑1 (MQ=255)
                                                                                             gTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGt                                                   <  2:1425592/149‑1 (MQ=255)
                                                                                             gTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGt                                                   <  2:3185731/149‑1 (MQ=255)
                                                                                                gCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGTCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGACATGCTGTCACATGTGCGGCAAAGTTTTAGTGcc                                                <  2:774316/149‑1 (MQ=255)
                                                                                                gCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCTATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGcc                                                <  1:2600482/149‑1 (MQ=255)
                                                                                                  cGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCaa                                              <  1:3965496/149‑1 (MQ=255)
                                                                                                   gATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAAc                                             >  2:773047/1‑149 (MQ=255)
                                                                                                        ggAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAgg                                        >  1:2434258/1‑149 (MQ=255)
                                                                                                            gCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGc                                    >  2:1594525/1‑149 (MQ=255)
                                                                                                            gCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGc                                    >  1:161861/1‑149 (MQ=255)
                                                                                                             ccGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCt                                   <  1:1542567/149‑1 (MQ=255)
                                                                                                              cGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTa                                  <  2:1061534/149‑1 (MQ=255)
                                                                                                              cGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTa                                  >  2:2460502/1‑149 (MQ=255)
                                                                                                              cGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTa                                  <  1:2570972/149‑1 (MQ=255)
                                                                                                              cGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTa                                  <  2:1657113/149‑1 (MQ=255)
                                                                                                               ggTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTAt                                 <  2:4098499/149‑1 (MQ=255)
                                                                                                                    tGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAGGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGa                            <  2:318917/149‑1 (MQ=255)
                                                                                                                      gCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGaaa                          <  1:1042158/149‑1 (MQ=255)
                                                                                                                      gCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGaaa                          >  1:3823234/1‑149 (MQ=255)
                                                                                                                          tgtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTgg                      <  1:2807817/149‑1 (MQ=255)
                                                                                                                           gtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGt                     <  1:417326/149‑1 (MQ=255)
                                                                                                                           gtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGt                     >  1:2075631/1‑149 (MQ=255)
                                                                                                                           gtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGt                     >  2:853422/1‑149 (MQ=255)
                                                                                                                               tGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTcgcg                 >  1:1664968/1‑149 (MQ=255)
                                                                                                                                gAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCgcgg                >  2:1840766/1‑149 (MQ=255)
                                                                                                                                      cTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTc          <  1:661893/149‑1 (MQ=255)
                                                                                                                                          gAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGc      >  2:520807/1‑149 (MQ=255)
                                                                                                                                           aaCTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCa     >  1:3799822/1‑149 (MQ=255)
                                                                                                                                              tgattaAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCt  >  1:3575374/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCT  >  CP009273/4207406‑4207696

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: