Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,886,725 G→T 86.9% I256I (ATC→ATA yoaA ← putative ATP‑dependent helicase, DinG family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,886,7250GT86.9% 628.7 / 98.2 258I256I (ATC→ATAyoaAputative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base G (19/15);  new base T (100/124);  total (119/139)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

TGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTT  >  CP009273/1886579‑1886870
                                                                                                                                                  |                                                                                                                                                 
tGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTAt                                                                                                                                                 <  2:458690/149‑1 (MQ=255)
 gAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATg                                                                                                                                                >  2:1459455/1‑149 (MQ=255)
 gAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATg                                                                                                                                                >  1:2125443/1‑149 (MQ=255)
 gAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATg                                                                                                                                                >  1:264284/1‑149 (MQ=255)
  aaTTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATgg                                                                                                                                               >  1:3451065/1‑149 (MQ=255)
  aaTTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATgg                                                                                                                                               <  1:3611551/149‑1 (MQ=255)
   aTTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGt                                                                                                                                              <  1:1961450/149‑1 (MQ=255)
    tttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTg                                                                                                                                             >  1:3020655/1‑149 (MQ=255)
     ttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGa                                                                                                                                            <  2:4498541/149‑1 (MQ=255)
         ggATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTc                                                                                                                                        >  2:2174591/1‑149 (MQ=255)
         ggATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTc                                                                                                                                        >  2:3592626/1‑149 (MQ=255)
         ggATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTc                                                                                                                                        >  2:3641063/1‑149 (MQ=255)
          gATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCt                                                                                                                                       <  1:4484030/149‑1 (MQ=255)
          gATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCt                                                                                                                                       >  2:405165/1‑149 (MQ=255)
           aTTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCtt                                                                                                                                      >  2:1746624/1‑149 (MQ=255)
            ttAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCttt                                                                                                                                     <  2:2400939/149‑1 (MQ=255)
             tAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCtttt                                                                                                                                    <  1:4005787/149‑1 (MQ=255)
               gCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGc                                                                                                                                  >  2:181227/1‑149 (MQ=255)
               gCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGc                                                                                                                                  >  2:266208/1‑149 (MQ=255)
                cTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGcc                                                                                                                                 <  2:3371113/149‑1 (MQ=255)
                cTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGcc                                                                                                                                 <  2:1199270/149‑1 (MQ=255)
                 tAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTTTGGTGATGTCTTTTGCCa                                                                                                                                <  2:606893/149‑1 (MQ=255)
                 tAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCa                                                                                                                                <  2:2636361/149‑1 (MQ=255)
                  aaCAGCTCACGCAGGTAACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAg                                                                                                                               >  2:839926/1‑149 (MQ=255)
                     aGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTc                                                                                                                            <  1:1419767/149‑1 (MQ=255)
                       cTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGa                                                                                                                          <  1:3664654/149‑1 (MQ=255)
                           cGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAg                                                                                                                      >  2:243340/1‑149 (MQ=255)
                               ggTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGt                                                                                                                  >  2:1981665/1‑149 (MQ=255)
                                gTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTc                                                                                                                 <  1:292386/149‑1 (MQ=255)
                                   aCCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGAc                                                                                                              <  2:3713799/149‑1 (MQ=255)
                                     cgcgATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTg                                                                                                            <  1:85588/149‑1 (MQ=255)
                                     cgcgATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTg                                                                                                            <  1:1650778/149‑1 (MQ=255)
                                        gATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGgag                                                                                                         >  1:287254/1‑149 (MQ=255)
                                        gATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGgag                                                                                                         <  1:166744/149‑1 (MQ=255)
                                        gATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGgag                                                                                                         >  1:1333063/1‑149 (MQ=255)
                                        gATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGgag                                                                                                         <  1:1254373/149‑1 (MQ=255)
                                         aTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGgaga                                                                                                        >  1:136923/1‑149 (MQ=255)
                                          tAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGagag                                                                                                       <  2:3487283/149‑1 (MQ=255)
                                           aaCCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGt                                                                                                      >  2:2795850/1‑149 (MQ=255)
                                           aaCCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGt                                                                                                      <  1:2285301/149‑1 (MQ=255)
                                                  cTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGtgactgac                                                                                               <  2:495318/149‑1 (MQ=255)
                                                   tCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGAcc                                                                                              <  2:3073620/149‑1 (MQ=255)
                                                   tCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGAcc                                                                                              >  2:1016276/1‑149 (MQ=255)
                                                    cACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCa                                                                                             >  1:3259603/1‑149 (MQ=255)
                                                         aGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATa                                                                                        <  2:949407/149‑1 (MQ=255)
                                                         aGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATa                                                                                        <  1:811159/149‑1 (MQ=255)
                                                             gCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActgg                                                                                    <  2:2910290/149‑1 (MQ=255)
                                                               aGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggct                                                                                  >  1:4197854/1‑149 (MQ=255)
                                                               aGACGAAAATCCTGCGCACGCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATActggct                                                                                  <  2:2853165/149‑1 (MQ=255)
                                                                gACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCACAATActggctg                                                                                 >  2:4401298/1‑149 (MQ=255)
                                                                   gAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCa                                                                              <  1:3185308/149‑1 (MQ=255)
                                                                    aaaaTCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCaa                                                                             >  1:1455265/1‑149 (MQ=255)
                                                                    aaaaTCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCaa                                                                             <  2:4073335/149‑1 (MQ=255)
                                                                    aaaaTCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCaa                                                                             >  2:3024653/1‑149 (MQ=255)
                                                                           tGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCgg                                                                      <  1:4352325/149‑1 (MQ=255)
                                                                           tGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCgg                                                                      <  1:3182212/149‑1 (MQ=255)
                                                                           tGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCgg                                                                      <  2:1163343/149‑1 (MQ=255)
                                                                           tGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCgg                                                                      <  2:4424348/149‑1 (MQ=255)
                                                                               cACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGc                                                                  <  2:3886907/149‑1 (MQ=255)
                                                                               cACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGc                                                                  <  1:534163/149‑1 (MQ=255)
                                                                                aCTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCt                                                                 >  1:2454405/1‑149 (MQ=255)
                                                                                        cAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCt                                                         <  2:2850426/149‑1 (MQ=255)
                                                                                        cAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCt                                                         >  1:1170887/1‑149 (MQ=255)
                                                                                          aGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTtc                                                       <  1:4003690/149‑1 (MQ=255)
                                                                                               aTCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGa                                                  >  2:1945380/1‑149 (MQ=255)
                                                                                                tCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGaa                                                 >  1:1911374/1‑149 (MQ=255)
                                                                                                  aGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGa                                               <  1:2654554/149‑1 (MQ=255)
                                                                                                  aGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCCTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGa                                               <  1:454761/149‑1 (MQ=255)
                                                                                                   gcgcACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTTGGCGAGGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGAt                                              <  2:113630/149‑1 (MQ=255)
                                                                                                   gcgcACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGTGCTTCGTCGAAGAt                                              >  2:1707901/1‑149 (MQ=255)
                                                                                                    cgcACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATc                                             <  1:950046/149‑1 (MQ=255)
                                                                                                      cACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCAt                                           >  1:2437583/1‑149 (MQ=255)
                                                                                                        cTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGa                                         >  1:1042309/1‑149 (MQ=255)
                                                                                                          tttGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACg                                       >  2:2755130/1‑149 (MQ=255)
                                                                                                           ttGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGt                                      >  1:3865166/1‑149 (MQ=255)
                                                                                                             gTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTcc                                    >  1:4102218/1‑149 (MQ=255)
                                                                                                               aaCTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGc                                  <  2:715991/149‑1 (MQ=255)
                                                                                                                aCTGCTGGGTGTCTTTTTATTCGGTGCGGTTGGCTATGGTGAGGTCTTTTGCCTGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCATTGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCt                                 <  2:1873275/149‑1 (MQ=255)
                                                                                                                aCTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGTCTGGAGAGTGACTGACCATAATACTGGCTGGCAATGTCCTGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCt                                 <  1:3791948/149‑1 (MQ=255)
                                                                                                                aCTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCt                                 <  1:2032708/149‑1 (MQ=255)
                                                                                                                 ctgctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCtt                                <  1:3441031/149‑1 (MQ=255)
                                                                                                                  tgctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTc                               <  2:79805/149‑1 (MQ=255)
                                                                                                                  tgctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTc                               <  1:2031619/149‑1 (MQ=255)
                                                                                                                   gctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTcc                              >  2:4604643/1‑149 (MQ=255)
                                                                                                                    ctgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCg                             <  2:1671170/149‑1 (MQ=255)
                                                                                                                     tgGGTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCgg                            >  2:1279087/1‑149 (MQ=255)
                                                                                                                       ggTGTCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGa                          >  1:2266887/1‑149 (MQ=255)
                                                                                                                       ggTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGa                          <  1:1126210/149‑1 (MQ=255)
                                                                                                                        gtgtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGAt                         >  2:2288032/1‑149 (MQ=255)
                                                                                                                        gtgtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGAt                         <  1:3463487/149‑1 (MQ=255)
                                                                                                                         tgtCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATc                        >  1:248449/1‑149 (MQ=255)
                                                                                                                           tCTTTTAATTCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAg                      >  1:1181854/1‑149 (MQ=255)
                                                                                                                             ttttAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGtt                    >  2:3825862/1‑149 (MQ=255)
                                                                                                                                    tCGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCaaa             >  1:3210850/1‑149 (MQ=255)
                                                                                                                                     cGGTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAAt            <  2:4285798/149‑1 (MQ=255)
                                                                                                                                      ggTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATc           <  1:4042689/149‑1 (MQ=255)
                                                                                                                                      ggTGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATc           <  1:3638743/149‑1 (MQ=255)
                                                                                                                                        tGCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCa         >  2:2558526/1‑149 (MQ=255)
                                                                                                                                         gCGGTAGGCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAACTACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCAc        <  2:4265340/149‑1 (MQ=255)
                                                                                                                                               ggCTATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCtt  >  2:1594000/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACTGGAGAGTGACTGACCAAAATACTGGCTGGCAATGTCCGGTAGCTGGTGGGCTTCGTCGAAGATCATGACGTCCGCTTCCGGGATCAGTTCGCCAAATCCACTCTCTT  >  CP009273/1886579‑1886870

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: