Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,422,304 G→T 9.5% intergenic (‑183/‑291) rrsD ← / → yrdA 16S ribosomal RNA of rrnD operon/bacterial transferase hexapeptide domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,422,3040GT 9.5% 604.2 / 28.8 175intergenic (‑183/‑291)rrsD/yrdA16S ribosomal RNA of rrnD operon/bacterial transferase hexapeptide domain protein
Reads supporting (aligned to +/- strand):  ref base G (64/94);  new base T (6/10);  total (71/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 177 total aligned reads displayed.

CGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGCACAAGTATTTTTTTG  >  CP009273/3422206‑3422453
                                                                                                  |                                                                                                                                                     
cGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAttttttt                                                                                                    <  1:1164015/150‑1 (MQ=31)
cGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGAttttttt                                                                                                    <  1:1164019/150‑1 (MQ=31)
  ttAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTtct                                                                                                  >  1:2551820/1‑150 (MQ=32)
    aaGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTGAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTctct                                                                                                >  1:1149050/1‑150 (MQ=18)
    aaGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTctct                                                                                                >  1:113340/1‑150 (MQ=33)
    aaGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTctct                                                                                                <  2:2491051/150‑1 (MQ=33)
      gAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTc                                                                                              <  2:724316/150‑1 (MQ=35)
       aaTCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCGCTTCg                                                                                             >  2:2178235/1‑150 (MQ=21)
         tCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgtt                                                                                           <  1:158881/150‑1 (MQ=38)
          ccGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCgttg                                                                                          <  2:2996767/150‑1 (MQ=38)
             tATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTc                                                                                       <  1:756441/150‑1 (MQ=40)
                  tCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGAcc                                                                                  >  1:1460538/1‑150 (MQ=255)
                  tCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGAcc                                                                                  >  1:1462176/1‑150 (MQ=255)
                   cGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCTTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCa                                                                                 >  2:1045504/1‑150 (MQ=255)
                   cGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCTACCa                                                                                 >  2:3515835/1‑150 (MQ=255)
                   cGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCa                                                                                 >  2:3515865/1‑150 (MQ=255)
                   cGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCa                                                                                 <  1:3002906/150‑1 (MQ=255)
                   cGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCAGATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCa                                                                                 >  2:3515873/1‑150 (MQ=255)
                    gAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCAt                                                                                <  1:3886729/150‑1 (MQ=255)
                        gCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCtg                                                                            >  2:1375328/1‑150 (MQ=255)
                        gCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACAATCAtg                                                                            >  2:3715674/1‑150 (MQ=255)
                                 gATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGtt                                                                   >  1:1762831/1‑150 (MQ=255)
                                  aTTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGttt                                                                  >  1:3210882/1‑150 (MQ=255)
                                   ttGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTc                                                                 <  1:829751/150‑1 (MQ=255)
                                   ttGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTc                                                                 <  1:981427/150‑1 (MQ=255)
                                   ttGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTc                                                                 <  1:830557/150‑1 (MQ=255)
                                   ttGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTc                                                                 <  1:830344/150‑1 (MQ=255)
                                     gTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCAc                                                               <  2:1766450/150‑1 (MQ=255)
                                       cTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGt                                                             >  1:3384423/1‑150 (MQ=255)
                                         gATATATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTg                                                           <  1:2444159/150‑1 (MQ=255)
                                          aTAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTgt                                                          <  1:453899/150‑1 (MQ=255)
                                           tAAATTGTTAAAGAGCTGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCTTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTgtc                                                         <  2:4086031/150‑1 (MQ=255)
                                           tAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTgtc                                                         <  1:3210836/150‑1 (MQ=255)
                                           tAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTgtc                                                         <  2:2771259/150‑1 (MQ=255)
                                                tGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGtctc                                                    <  2:1296325/150‑1 (MQ=255)
                                                   tAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAAc                                                 <  2:658577/150‑1 (MQ=255)
                                                   tAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAAc                                                 >  1:4448528/1‑150 (MQ=255)
                                                   tAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAAc                                                 >  2:148644/1‑150 (MQ=255)
                                                        agCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCTTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAgg                                            >  2:1496449/1‑150 (MQ=255)
                                                        agCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAgg                                            <  1:3404086/150‑1 (MQ=255)
                                                        agCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAgg                                            >  2:2667009/1‑150 (MQ=255)
                                                          cAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGgcg                                          <  2:30127/150‑1 (MQ=255)
                                                           aGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGgcgc                                         >  1:2232975/1‑150 (MQ=255)
                                                           aGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGgcgc                                         >  1:575015/1‑150 (MQ=255)
                                                            gTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCa                                        >  1:4227574/1‑150 (MQ=255)
                                                             ttGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCAt                                       <  2:2968892/150‑1 (MQ=255)
                                                             ttGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCAt                                       >  1:3085665/1‑150 (MQ=255)
                                                             ttGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCAt                                       >  1:903765/1‑150 (MQ=255)
                                                             ttGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCAt                                       >  1:2428609/1‑150 (MQ=255)
                                                             ttGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGATGTCGTCTCAACGGAGGCGCAt                                       >  1:2428611/1‑150 (MQ=255)
                                                              tGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCagt                                      >  1:2656251/1‑148 (MQ=255)
                                                              tGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCCtt                                      >  1:2657445/1‑150 (MQ=255)
                                                               gCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATta                                     <  2:2149815/150‑1 (MQ=255)
                                                               gCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATta                                     <  2:4391780/150‑1 (MQ=255)
                                                               gCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATta                                     <  2:943398/150‑1 (MQ=255)
                                                                cGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATtat                                    <  1:3918145/150‑1 (MQ=255)
                                                                cGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATtat                                    <  2:458438/150‑1 (MQ=255)
                                                                cGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATtat                                    <  2:458096/150‑1 (MQ=255)
                                                                cGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATtat                                    <  1:4256059/150‑1 (MQ=255)
                                                                cGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATtat                                    <  1:294562/150‑1 (MQ=255)
                                                                 gACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATtata                                   >  2:2842013/1‑150 (MQ=255)
                                                                 gACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATtata                                   <  1:3761042/150‑1 (MQ=255)
                                                                    gcgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAggg                                >  2:3083603/1‑150 (MQ=255)
                                                                    gcgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAggg                                >  2:3083382/1‑150 (MQ=255)
                                                                    gcgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAAGTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAggg                                >  2:3083715/1‑150 (MQ=255)
                                                                     cgGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCTTTATAGGGa                               >  2:4110996/1‑150 (MQ=255)
                                                                      ggCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGAt                              >  2:2081744/1‑150 (MQ=255)
                                                                      ggCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGAt                              >  2:2081213/1‑150 (MQ=255)
                                                                       gCTTTCAGCTCACTGTCGCGAGGTGGCTTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATc                             <  1:1302674/150‑1 (MQ=255)
                                                                       gCTTTCAGCTCACTGTCGCGAGGTGGCTTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATc                             <  1:1302234/150‑1 (MQ=255)
                                                                       gCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATc                             <  2:3293603/150‑1 (MQ=255)
                                                                       gCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATc                             <  2:2273953/150‑1 (MQ=255)
                                                                       gCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATc                             <  2:3293653/150‑1 (MQ=255)
                                                                       gCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATc                             <  2:3695024/150‑1 (MQ=255)
                                                                        cTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTCCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATcc                            <  1:3977115/150‑1 (MQ=255)
                                                                         tttCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATccc                           >  1:1048738/1‑150 (MQ=255)
                                                                         tttCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATccc                           >  1:1048798/1‑150 (MQ=255)
                                                                         tttCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTACGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATccc                           >  2:3182774/1‑150 (MQ=255)
                                                                            cAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAAt                        <  2:4297472/150‑1 (MQ=255)
                                                                            cAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAAt                        <  1:2314348/150‑1 (MQ=255)
                                                                             aGCTCACTGTCGCGAGGTGGCTTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAAtt                       <  1:565757/150‑1 (MQ=255)
                                                                             aGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCaagt                       <  2:3574051/150‑3 (MQ=255)
                                                                              gCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAAttt                      <  2:1718373/150‑1 (MQ=255)
                                                                              gCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAAttt                      <  1:1340863/150‑1 (MQ=255)
                                                                               cTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAAtttt                     <  2:1011954/150‑1 (MQ=255)
                                                                               cTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAAtttt                     <  2:1011873/150‑1 (MQ=255)
                                                                               cTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAAtttt                     <  2:3843279/150‑1 (MQ=255)
                                                                               cTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAAtttt                     <  1:531274/150‑1 (MQ=255)
                                                                               cTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAAtttt                     <  1:530237/150‑1 (MQ=255)
                                                                                 cACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAAtttttt                   <  1:680919/150‑1 (MQ=255)
                                                                                 cACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAAtttttt                   <  1:3432627/150‑1 (MQ=255)
                                                                                  aCTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTg                  >  2:629617/1‑150 (MQ=255)
                                                                                    tGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGca                >  1:3422673/1‑150 (MQ=255)
                                                                                        gcgAGGTGGCTTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGCACAAg            <  2:242801/150‑1 (MQ=255)
                                                                                        gcgAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGCACAAg            <  1:2042873/150‑1 (MQ=255)
                                                                                        gcgAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGCACAAg            <  1:2042893/150‑1 (MQ=255)
                                                                                           aGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGCACAAGTAt         >  1:696851/1‑150 (MQ=255)
                                                                                               ggCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGCACAAGTAttttt     >  2:3756427/1‑150 (MQ=255)
                                                                                                  gTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGCACAAGTATTTTTTTg  >  2:3949625/1‑150 (MQ=255)
                                                                                                  gTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGCACAAGTATTTTTTTg  >  2:1692502/1‑150 (MQ=255)
                                                                                                  |                                                                                                                                                     
CGTTAAGAATCCGTATCTTCGAGTGCCCACACAGATTGTCTGATAAATTGTTAAAGAGCAGTTGCGACGCGGCTTTCAGCTCACTGTCGCGAGGTGGCGTATATTACGCTTTCCTCTTTCAGAGTCAACCCTGAATTTCAGGATTTTTTTCTCTTCGTTGTTCCGACCATCCTGTGAAGTGTTTCACGTTGTCGTCTCAACGGAGGCGCATTATAGGGATCCCAATTTTTTGCACAAGTATTTTTTTG  >  CP009273/3422206‑3422453

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: