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,5520GT98.6% 981.8 / 2.5 289M172I (ATG→ATTaceAisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base T (134/151);  total (134/155)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-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 293 total aligned reads displayed.

ATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCTG  >  CP009273/4207406‑4207697
                                                                                                                                                  |                                                                                                                                                 
aTCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTaa                                                                                                                                                 >  1:517954/1‑149 (MQ=255)
 tCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTaaa                                                                                                                                                >  1:2848652/1‑149 (MQ=255)
 tCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTaaa                                                                                                                                                <  2:1710737/149‑1 (MQ=255)
 tCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTaaa                                                                                                                                                <  2:3787203/149‑1 (MQ=255)
 tCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTaaa                                                                                                                                                >  1:712578/1‑149 (MQ=255)
  caacaaCACCTTCCGTCGTGCCGATCAGATCCAAGGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAg                                                                                                                                               <  1:1222465/149‑1 (MQ=255)
   aacaacACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGc                                                                                                                                              >  1:2580454/1‑149 (MQ=255)
    acaacaCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCg                                                                                                                                             >  2:2983706/1‑149 (MQ=255)
    acaacaCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCg                                                                                                                                             <  2:1852727/149‑1 (MQ=255)
     caacaCCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCga                                                                                                                                            >  2:2166641/1‑149 (MQ=255)
       acacCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCgatg                                                                                                                                          >  2:4020258/1‑149 (MQ=255)
       acacCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCgatg                                                                                                                                          >  2:3376306/1‑149 (MQ=255)
       acacCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCgatg                                                                                                                                          <  2:2281858/149‑1 (MQ=255)
       acacCGTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCgatg                                                                                                                                          <  2:4679225/149‑1 (MQ=255)
          ccTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGAtt                                                                                                                                       <  1:2499629/149‑1 (MQ=255)
           cTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTg                                                                                                                                      <  2:1699910/149‑1 (MQ=255)
             tCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGaa                                                                                                                                    <  2:282135/149‑1 (MQ=255)
                gtcgtGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGTCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCTATGAGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGcc                                                                                                                                 <  1:3752715/149‑1 (MQ=255)
                gtcgtGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGcc                                                                                                                                 <  2:4188416/149‑1 (MQ=255)
                   gtGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGt                                                                                                                              >  1:1088254/1‑149 (MQ=255)
                    tGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTg                                                                                                                             >  2:3523866/1‑149 (MQ=255)
                     gCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGc                                                                                                                            <  2:966187/149‑1 (MQ=255)
                         aTCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCg                                                                                                                        >  2:2493874/1‑149 (MQ=255)
                              aTCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGt                                                                                                                   <  2:1117944/149‑1 (MQ=255)
                               tCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGtt                                                                                                                  <  2:4718018/149‑1 (MQ=255)
                                 cAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCa                                                                                                                >  2:768144/1‑149 (MQ=255)
                                  aaTGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCAc                                                                                                               >  1:1022769/1‑149 (MQ=255)
                                  aaTGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCAc                                                                                                               <  1:21659/149‑1 (MQ=255)
                                   aTGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACt                                                                                                              >  2:1004079/1‑149 (MQ=255)
                                    tGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACtt                                                                                                             >  1:1867096/1‑149 (MQ=255)
                                    tGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACtt                                                                                                             >  2:3179805/1‑149 (MQ=255)
                                     ggTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTc                                                                                                            >  1:371864/1‑149 (MQ=255)
                                                 ttGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTg                                                                                                >  1:3214546/1‑149 (MQ=255)
                                                 ttGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTg                                                                                                >  1:382645/1‑149 (MQ=255)
                                                     gCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGt                                                                                            <  2:3474731/149‑1 (MQ=255)
                                                     gCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGt                                                                                            >  1:427487/1‑149 (MQ=255)
                                                        gggCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAg                                                                                         >  1:2673733/1‑149 (MQ=255)
                                                        gggCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAg                                                                                         >  1:1801552/1‑149 (MQ=255)
                                                          gCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTg                                                                                       <  2:4602571/149‑1 (MQ=255)
                                                            gATGCGCGCTATGTCGAGTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGAGGTGTCCTGAATGCCTTTGAACTGATTATAGCGATGATTGATGCAGGTGCAGCTGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTgaa                                                                                     <  2:180102/149‑1 (MQ=255)
                                                                      aTGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGt                                                                           >  2:4104283/1‑149 (MQ=255)
                                                                      aTGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGt                                                                           >  1:622409/1‑149 (MQ=255)
                                                                       tGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTc                                                                          <  1:2186327/149‑1 (MQ=255)
                                                                       tGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTc                                                                          >  1:3305842/1‑149 (MQ=255)
                                                                        gTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTca                                                                         >  1:1017129/1‑149 (MQ=255)
                                                                        gTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTca                                                                         <  1:4148350/149‑1 (MQ=255)
                                                                        gTCGATTACTTCCTGCCGATCGGTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTca                                                                         <  1:2114693/149‑1 (MQ=255)
                                                                         tCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTcac                                                                        >  2:3538071/1‑149 (MQ=255)
                                                                          cGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTcaca                                                                       <  2:3386597/149‑1 (MQ=255)
                                                                           gATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACAt                                                                      >  2:2448058/1‑149 (MQ=255)
                                                                           gATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACAt                                                                      >  2:1633374/1‑149 (MQ=255)
                                                                            aTTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATg                                                                     <  1:1416206/149‑1 (MQ=255)
                                                                            aTTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATg                                                                     <  2:679709/149‑1 (MQ=255)
                                                                             ttACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATgg                                                                    <  2:533378/149‑1 (MQ=255)
                                                                               aCTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGggc                                                                  <  2:307341/149‑1 (MQ=255)
                                                                                        cGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtt                                                         >  1:241521/1‑149 (MQ=255)
                                                                                          aTCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGtttt                                                       >  1:2258262/1‑149 (MQ=255)
                                                                                            cGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAg                                                     <  1:2522690/149‑1 (MQ=255)
                                                                                             ggtGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGt                                                    <  2:1266928/147‑1 (MQ=255)
                                                                                             gTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGt                                                    <  1:3149939/149‑1 (MQ=255)
                                                                                             gTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTGTTAGt                                                    >  2:489617/1‑149 (MQ=255)
                                                                                                    aTGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACt                                             >  2:579345/1‑149 (MQ=255)
                                                                                                      gCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCa                                           <  1:2444685/149‑1 (MQ=255)
                                                                                                      gCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCa                                           <  1:2881496/149‑1 (MQ=255)
                                                                                                      gCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCa                                           <  1:4589407/149‑1 (MQ=255)
                                                                                                       cGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAg                                          <  2:3216236/149‑1 (MQ=255)
                                                                                                        ggAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAgg                                         >  1:3629313/1‑149 (MQ=255)
                                                                                                           aGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAg                                      >  1:3938641/1‑149 (MQ=255)
                                                                                                           aGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAg                                      >  2:2739966/1‑149 (MQ=255)
                                                                                                            gCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGc                                     <  1:2222864/149‑1 (MQ=255)
                                                                                                            gCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGc                                     <  2:2874112/149‑1 (MQ=255)
                                                                                                               ggTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTAt                                  <  2:3028989/149‑1 (MQ=255)
                                                                                                                gTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTAtt                                 <  1:2739293/149‑1 (MQ=255)
                                                                                                                  tttGGCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCa                               <  2:2111247/149‑1 (MQ=255)
                                                                                                                     ggCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGaa                            >  1:890023/1‑149 (MQ=255)
                                                                                                                      gCGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGaaa                           <  1:929475/149‑1 (MQ=255)
                                                                                                                       cGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAAc                          <  2:2424657/149‑1 (MQ=255)
                                                                                                                       cGGTGTCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAAc                          <  2:3654338/149‑1 (MQ=255)
                                                                                                                         gtgtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTg                        >  1:3479627/1‑149 (MQ=255)
                                                                                                                         gtgtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTg                        >  2:4373884/1‑149 (MQ=255)
                                                                                                                          tgtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTgg                       <  2:4042520/149‑1 (MQ=255)
                                                                                                                          tgtCCTGAATGCCTGTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTgg                       <  2:766115/149‑1 (MQ=255)
                                                                                                                           gtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGt                      >  1:2804997/1‑149 (MQ=255)
                                                                                                                           gtCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGt                      >  2:4656651/1‑149 (MQ=255)
                                                                                                                            tCCTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTc                     >  1:2919289/1‑149 (MQ=255)
                                                                                                                              cTGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTcgc                   <  2:1503530/149‑1 (MQ=255)
                                                                                                                               tGAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTcgcg                  >  2:332739/1‑149 (MQ=255)
                                                                                                                                gAATGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCgcgg                 <  1:349226/149‑1 (MQ=255)
                                                                                                                                 aaTGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCgcggc                >  2:3812137/1‑149 (MQ=255)
                                                                                                                                 aaTGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCgcggc                >  1:3865359/1‑149 (MQ=255)
                                                                                                                                   tGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGgcgc              >  1:1506707/1‑149 (MQ=255)
                                                                                                                                   tGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGgcgc              >  2:4351639/1‑149 (MQ=255)
                                                                                                                                   tGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGgcgc              <  1:3576210/149‑1 (MQ=255)
                                                                                                                                   tGCCTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAAATGGTCGCGgcgc              >  1:4256939/1‑149 (MQ=255)
                                                                                                                                      cTTTGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAACTGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTc           >  1:517094/1‑149 (MQ=255)
                                                                                                                                       tttGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCCGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCt          <  2:2673733/149‑1 (MQ=255)
                                                                                                                                       tttGAACTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCt          >  2:1938942/1‑149 (MQ=255)
                                                                                                                                           aaCTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCa      <  1:982314/149‑1 (MQ=255)
                                                                                                                                             cTGATTAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGc    <  1:3376306/149‑1 (MQ=255)
                                                                                                                                               gattaAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCtg  >  1:4362104/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACTTCCTGCCGATCGTTGCCGATGCGGAAGCCGGTTTTGGCGGTGTCCTGAATGCCTTTGAACTGATGAAAGCGATGATTGAAGCCGGTGCAGCGGCAGTTCACTTCGAAGATCAGCTGGCGTCAGTGAAGAAATGCGGTCACATGGGCGGCAAAGTTTTAGTGCCAACTCAGGAAGCTATTCAGAAACTGGTCGCGGCGCGTCTGGCAGCTG  >  CP009273/4207406‑4207697

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: