Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,158,047 G→T 20.6% noncoding (1461/1542 nt) rrsB → 16S ribosomal RNA of rrnB operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,158,0470GT20.6% 242.7 / 52.6 107noncoding (1461/1542 nt)rrsB16S ribosomal RNA of rrnB operon
Reads supporting (aligned to +/- strand):  ref base G (51/34);  new base T (11/11);  total (62/45)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 110 total aligned reads displayed.

TGGGAGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTA  >  CP009273/4158000‑4158196
                                               |                                                                                                                                                     
tGGGAGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGt                                                 >  1:161519/1‑150 (MQ=11)
 gggAGTGTGTTGCAAAAGAAGTTGGTAGCTTAACCTTCGTGAGGGCTCTTACAACTTTGTGTTTCATGACTTGGGTGAAGTCGTAACAATGTAACCGTAGTTGAACCTGCGTTTTGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTa                                                <  1:3267986/150‑1 (MQ=11)
 gggAGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTa                                                <  1:3268005/150‑1 (MQ=14)
 gggAGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTa                                                >  1:1874098/1‑150 (MQ=31)
 gggAGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTa                                                <  2:3424055/150‑1 (MQ=31)
  ggAGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAg                                               <  2:2879674/150‑1 (MQ=32)
  ggAGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAg                                               >  2:3903512/1‑150 (MQ=32)
   gAGTGGGTTGCAAATGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGt                                              <  1:4395088/150‑1 (MQ=17)
   gAGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGt                                              <  2:2956313/150‑1 (MQ=32)
   gAGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGt                                              <  1:4395033/150‑1 (MQ=32)
    aGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATTACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTg                                             >  1:1531277/1‑150 (MQ=17)
    aGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTg                                             >  2:1620161/1‑150 (MQ=32)
    aGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTg                                             >  1:4212790/1‑150 (MQ=32)
      tGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCt                                           >  2:3586623/1‑150 (MQ=32)
      tGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCt                                           >  1:2466835/1‑150 (MQ=32)
        ggTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTca                                         >  2:1970066/1‑150 (MQ=17)
         gTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcac                                        >  2:4079037/1‑150 (MQ=17)
         gTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcac                                        >  1:301327/1‑150 (MQ=32)
         gTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcac                                        >  2:2248922/1‑150 (MQ=32)
          ttGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcaca                                       <  2:2225375/150‑1 (MQ=17)
          ttGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAATGAAGCGTACTTTGTAGTGCTcaca                                       >  1:679446/1‑150 (MQ=17)
          ttGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTGGTAGTGCTcaca                                       >  1:679507/1‑150 (MQ=17)
          tagCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGTTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcaca                                       >  1:679520/3‑150 (MQ=17)
           tGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcacac                                      <  1:4357055/150‑1 (MQ=32)
           tGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcacac                                      <  1:3520881/150‑1 (MQ=32)
           tGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcacac                                      <  1:3521527/150‑1 (MQ=32)
            gCAAACGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcacaca                                     <  2:2101018/150‑1 (MQ=17)
            gCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcacaca                                     >  1:3099389/1‑150 (MQ=32)
            gCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTcacaca                                     >  2:4736949/1‑150 (MQ=32)
             cAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAg                                    >  2:3834036/1‑150 (MQ=17)
             cAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAg                                    >  2:4458934/1‑150 (MQ=17)
             cAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAg                                    >  2:152525/1‑150 (MQ=32)
             cAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAg                                    <  1:4647957/150‑1 (MQ=32)
             cAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAg                                    >  2:2751206/1‑150 (MQ=32)
               aaaGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGAt                                  >  1:4177029/1‑150 (MQ=32)
               aaaGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGAt                                  <  1:2144154/150‑1 (MQ=32)
                aagaagTAGGTAGCTTAACCTTCGGGAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGAtt                                 <  2:2077366/150‑1 (MQ=17)
                aagaagTAGGTAGCTTAACCTTCGGGAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGAtt                                 <  2:2077843/150‑1 (MQ=17)
                aagaagTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGAtt                                 <  2:1306437/150‑1 (MQ=32)
                aagaagTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGAtt                                 <  1:790212/150‑1 (MQ=32)
                aagaagTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGAtt                                 >  2:3524029/1‑150 (MQ=32)
                aagaagTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGAtt                                 >  2:375389/1‑150 (MQ=32)
                aagaagTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGAtt                                 >  2:4845921/1‑150 (MQ=32)
                aagaagTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGAtt                                 >  2:953021/1‑150 (MQ=32)
                aagaagTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGAtt                                 >  1:4745015/1‑150 (MQ=32)
                aagaagTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTAATTTGTAGTGCTCACACAGAtt                                 >  2:953310/1‑150 (MQ=17)
                 agaagTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTg                                <  1:1706700/150‑1 (MQ=32)
                 agaagTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTg                                >  1:1059137/1‑150 (MQ=32)
                   aagTAGGTAGCTTAACCTTCGGGAGGGAGATTACAAGTTAGTGATACATGACTTGGGTGAAGTCGAAACAAGGTAAAAGTAGGGGAACCTGAGGTTGGATCACCTACTTACCTTAAAGAAGCGTACTTTGTATTGCACACACAGATAGta                              >  2:2890515/1‑149 (MQ=11)
                    agTAGGTAGCTTAACCTTCGGGAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                             >  1:3542663/1‑150 (MQ=17)
                    agTAGGTAGCTTAACCTTCGGGAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTCTGTAGTGCTCACACAGATTGTCt                             >  1:3542728/1‑150 (MQ=17)
                    agTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCt                             <  2:3314226/150‑1 (MQ=32)
                     gtaggtagCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTg                            >  2:2753890/1‑150 (MQ=32)
                     gtaggtagCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTg                            >  2:2753924/1‑150 (MQ=32)
                     gtaggtagCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTg                            <  1:918396/150‑1 (MQ=32)
                     gtaggtagCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTg                            <  2:5009191/150‑1 (MQ=32)
                      taggtagCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGa                           >  1:208241/1‑150 (MQ=32)
                       aggtagCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGAt                          <  1:5081637/150‑1 (MQ=32)
                          tagCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGa                       >  1:4816543/1‑150 (MQ=34)
                           agCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGaa                      >  1:3762354/1‑150 (MQ=34)
                           agCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGaa                      <  2:399918/150‑1 (MQ=34)
                            gCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGaaa                     >  2:105278/1‑150 (MQ=34)
                            gCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGaaa                     <  1:373860/150‑1 (MQ=34)
                             cTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAg                    <  2:2280276/150‑1 (MQ=34)
                              ttAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGt                   >  1:1114740/1‑150 (MQ=34)
                              ttAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGt                   >  1:1017155/1‑150 (MQ=34)
                              ttAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGt                   >  2:2508490/1‑150 (MQ=34)
                              ttAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGt                   <  2:3625678/150‑1 (MQ=34)
                              ttAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGt                   >  1:1017205/1‑150 (MQ=34)
                               tAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACGTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTg                  >  2:1797038/1‑150 (MQ=21)
                                aaCCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGa                 <  1:1316701/150‑1 (MQ=34)
                                aaCCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGa                 >  1:2571615/1‑150 (MQ=34)
                                  ccTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGaaa               >  1:1864531/1‑150 (MQ=34)
                                  ccTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGaaa               >  1:1864881/1‑150 (MQ=34)
                                  ccTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGaaa               >  1:1865076/1‑150 (MQ=34)
                                   cTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGaaaa              <  2:2412196/150‑1 (MQ=34)
                                   cTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGaaaa              >  2:3774100/1‑150 (MQ=34)
                                    ttCGGGAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAg             >  1:1627428/1‑150 (MQ=21)
                                     tCGGGAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGc            >  1:4538604/1‑150 (MQ=21)
                                     tCGGGAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGc            >  1:4538420/1‑150 (MQ=21)
                                     tCGGGAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGc            >  1:4538981/1‑150 (MQ=21)
                                     tCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGc            <  1:337158/150‑1 (MQ=34)
                                      cGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCa           >  1:3192490/1‑150 (MQ=34)
                                      cGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCa           >  1:3188961/1‑150 (MQ=34)
                                       gggAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCaa          <  1:993361/150‑1 (MQ=21)
                                       gggAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCaa          >  2:364191/1‑150 (MQ=34)
                                        ggAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAg         <  1:1838786/150‑1 (MQ=34)
                                        ggAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAg         <  2:4936227/150‑1 (MQ=34)
                                        ggAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACAGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAg         >  1:3766259/1‑150 (MQ=21)
                                         gAGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAgg        >  1:2808452/1‑150 (MQ=21)
                                         gAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAgg        >  2:3132740/1‑150 (MQ=34)
                                         gAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAgg        <  2:3894745/150‑1 (MQ=34)
                                         gAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAgg        <  1:3444111/150‑1 (MQ=34)
                                          aGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGc       <  2:3284469/150‑1 (MQ=21)
                                          aGGGCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGc       <  2:3284414/150‑1 (MQ=21)
                                           gggCTCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCg      <  1:2739643/150‑1 (MQ=21)
                                           gggCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCg      <  1:4489755/150‑1 (MQ=34)
                                            ggCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGt     >  2:3958113/1‑150 (MQ=34)
                                              cgcTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGttt   <  1:3236465/150‑1 (MQ=35)
                                               gcTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTa  <  1:1113264/150‑1 (MQ=37)
                                               |                                                                                                                                                     
TGGGAGTGGGTTGCAAAAGAAGTAGGTAGCTTAACCTTCGGGAGGGCGCTTACCACTTTGTGATTCATGACTGGGGTGAAGTCGTAACAAGGTAACCGTAGGGGAACCTGCGGTTGGATCACCTCCTTACCTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTA  >  CP009273/4158000‑4158196

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: