Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 4,029,087 A→T noncoding (197/1542 nt) rrsA → 16S ribosomal RNA of rrnA operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,029,0870AT98.4% 124.2 / NA 63noncoding (197/1542 nt)rrsA16S ribosomal RNA of rrnA operon
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (29/33);  total (29/34)
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.92e-01

CATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCCATGCCGCGTGTATGAAGAAGG  >  CP009273/4028855‑4029307
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                            
cATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGt                                                                                                                                                                                                                       <  2:314559/240‑1 (MQ=11)
                                  ttAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGa                                                                                                                                                                                     <  1:160840/240‑1 (MQ=14)
                                                   agcATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGAtt                                                                                                                                                                              <  1:2154/228‑1 (MQ=14)
                                                    tCATGGCTCAGATTGGACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAgg                                                                                                                                                                    <  1:6248/239‑1 (MQ=14)
                                                       tGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGa                                                                                                                                                                                     <  2:214206/219‑1 (MQ=14)
                                                         gCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGGGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGAGCCCAGATGGGATTAGCta                                                                                                                                                                         >  1:264513/1‑229 (MQ=14)
                                                                   aaCGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGCCGGGTGAGAAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGc                                                                                                                                                                                                             >  2:328297/1‑183 (MQ=14)
                                                                           cggcAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGa                                                                                                                                                                                <  2:277858/204‑1 (MQ=14)
                                                                                           tGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGCTGGAGGGGGATAACTACTGGCAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGt                                                                                                                             >  1:280650/1‑239 (MQ=14)
                                                                                           tGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTc                                                                                                                            >  1:399654/1‑240 (MQ=14)
                                                                                            gCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGcg                                                                                                                                             >  2:476176/1‑222 (MQ=14)
                                                                                                    aaCGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATg                                                                                                                                                                                   >  1:458606/1‑176 (MQ=14)
                                                                                                                 agaagCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCGGagag                                                                                                                      >  2:177819/1‑224 (MQ=14)
                                                                                                                  gaagCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCacac                                                                                                       <  2:434530/238‑1 (MQ=14)
                                                                                                                     gcttgcttCTTTGCTGACGAGTTGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCAt                                                                                                                                                                                                   <  2:228786/143‑1 (MQ=17)
                                                                                                                               ttGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                                                                                          <  1:371764/238‑1 (MQ=14)
                                                                                                                               tnGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc                                                                                                            <  2:132722/220‑1 (MQ=14)
                                                                                                                                tGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAAc                                                                                                 >  1:276351/1‑230 (MQ=14)
                                                                                                                                   tggcGAGTGGCGGACGGGTGAGTACTGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCt                                                                                                                                                       <  2:27375/170‑1 (MQ=9)
                                                                                                                                    gACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                                                                                          <  2:280650/233‑1 (MQ=14)
                                                                                                                                      cGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACGGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCa                                                                                   <  2:202685/238‑1 (MQ=14)
                                                                                                                                       gAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                                                                                          <  2:204524/230‑1 (MQ=14)
                                                                                                                                         gttgCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                                                                                          <  1:310743/225‑1 (MQ=14)
                                                                                                                                         gTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGc                                                                                                            <  2:423441/210‑1 (MQ=14)
                                                                                                                                                 cGGGGGAGTAATGTCTGGGAAACATCGTGATGGAGGGGGATAACTACGGGTAACGGTAGCTAATACCACATAACGTCGCAAGACCAAAGAGGGGTCCCTTCTGGCCTCTTGCCAGCGGATGTGCTTAGATAGGATTAGCt                                                                                                                                                                          <  2:56902/140‑1 (MQ=11)
                                                                                                                                                    gTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCAACCCTGGAACTGAGACACGGTCCAGACTCCTACGGGAggg                                                                   >  2:225503/1‑239 (MQ=14)
                                                                                                                                                           aTGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGa                                                                                 <  2:497694/219‑1 (MQ=14)
                                                                                                                                                           aTGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTg                                                            <  1:112279/240‑1 (MQ=14)
                                                                                                                                                             gTCTGGGCAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCGCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGacac                                                                                         <  2:350220/209‑1 (MQ=14)
                                                                                                                                                              tatGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCACAAGACCAATGCGGGGTACCTTCGGGCCTCTTGCCATCGTATGTGACCAGATGGGATAAGCTAGTAgg                                                                                                                                                                    >  2:385464/3‑133 (MQ=11)
                                                                                                                                                              tCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTggg                                                          <  1:354150/239‑1 (MQ=14)
                                                                                                                                                              tCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAACTGGTCTGAGGGGATGACCAGCCACACCGGAACTgaggc                                                                                           >  2:375503/1‑204 (MQ=9)
                                                                                                                                                              tCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGCATTAGCt                                                                                                                                                                          >  2:150289/1‑127 (MQ=17)
                                                                                                                                                                  ggAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGTATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGa                                                        <  1:247847/237‑1 (MQ=14)
                                                                                                                                                                   gAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGaa                                                       >  1:86245/1‑237 (MQ=14)
                                                                                                                                                                        tGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTCGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGac                                                                                           >  1:663/1‑196 (MQ=14)
                                                                                                                                                                        tGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTg                                                                                                     <  2:64042/186‑1 (MQ=14)
                                                                                                                                                                         gCCTGATGGAGGGTGAGTACTACTGGTAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTTGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCAAACTGGAACTGAGACACGGTCCATACTCCTACGGGAGGCAGCAGTggg                                                          <  2:340675/228‑1 (MQ=11)
                                                                                                                                                                             gATGGCGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTg                                                  <  2:183174/232‑1 (MQ=14)
                                                                                                                                                                             gATGGAGGTGGATAACTACTGGAAAGGGTAGCTAATACCGCAGAGCGTCGCAAGACAAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGac                                                                                           <  2:338964/191‑1 (MQ=255)
                                                                                                                                                                             gATGGAGGGGGATAACGACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGgcagca                                                               <  1:301325/219‑1 (MQ=14)
                                                                                                                                                                                 gAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGgcgc                                      <  1:198356/240‑1 (MQ=14)
                                                                                                                                                                                                   aaaCGGTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCCATg                    <  1:476176/240‑1 (MQ=14)
                                                                                                                                                                                                       ggTAGCTAATACCGCATAACGTCGCACGACCAATGAGGGGTACCTTCgg                                                                                                                                                                                                               >  2:359624/1‑49 (MQ=11)
                                                                                                                                                                                                       ggTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGt                                                                                                                                                                                                                       <  2:519952/41‑1 (MQ=255)
                                                                                                                                                                                                       ggTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGt                                                                                                                                                                                                                       >  1:519952/1‑41 (MQ=12)
                                                                                                                                                                                                       ggTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCaa                                    >  1:382917/1‑220 (MQ=14)
                                                                                                                                                                                                       ggTAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCCATGCcgc                >  1:210893/1‑240 (MQ=14)
                                                                                                                                                                                                         tAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAg                                                                                                                                                                                      <  1:60639/72‑1 (MQ=21)
                                                                                                                                                                                                         tAGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAg                                                                                                                                                                                      >  2:60639/1‑72 (MQ=255)
                                                                                                                                                                                                          aGCTAATACCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCCATGCCGCgtg             >  1:393735/1‑240 (MQ=14)
                                                                                                                                                                                                                 aCCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGaca                                                                                          <  2:328080/156‑1 (MQ=14)
                                                                                                                                                                                                                 aCCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAg                                                                     >  1:268038/1‑177 (MQ=14)
                                                                                                                                                                                                                 aCCGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCCATGCCGCGTGTATgaag      >  1:572323/1‑240 (MQ=14)
                                                                                                                                                                                                                  tcGCATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGcga                                                                                                                                            <  2:488808/104‑1 (MQ=17)
                                                                                                                                                                                                                     cATAACGTCGCAATACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTcc                                                                                    >  2:131328/1‑158 (MQ=14)
                                                                                                                                                                                                                     cATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGAt                                                                                                                                                                                              >  2:71943/1‑52 (MQ=255)
                                                                                                                                                                                                                     cATAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCCATGCCGCGTGTATGCAGAAgg  >  1:122751/1‑240 (MQ=14)
                                                                                                                                                                                                                      aTAACGTCGCAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCCATGc                   >  1:162993/1‑222 (MQ=14)
                                                                                                                                                                                                                          cGTCGCAAGACCAATGACGGGTACCTTCGGGCCTCTTGCCATCGGGTGTGCCCAGATGGGATTAGCTCGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCGCGCCTGAGCGGATGACCAGCCACCCTGGAAc                                                                                                 >  2:177545/1‑140 (MQ=11)
                                                                                                                                                                                                                               cAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAg                                                                                                                                                <  2:551746/88‑1 (MQ=18)
                                                                                                                                                                                                                               cAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAg                                                                                                                                                >  1:551746/1‑88 (MQ=25)
                                                                                                                                                                                                                               cAAGACCAATGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACCCTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGCATATTGCACAATGGgcgc                                      >  2:166228/1‑194 (MQ=14)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                            
CATTACGAAGTTTAATTCTTTGAGCGTCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGTACCTTCGGGCCTCTTGCCATCGGATGTGCCCAGATGGGATTAGCTAGTAGGTGGGGTAACGGCTCACCTAGGCGACGATCCCTAGCTGGTCTGAGAGGATGACCAGCCACACTGGAACTGAGACACGGTCCAGACTCCTACGGGAGGCAGCAGTGGGGAATATTGCACAATGGGCGCAAGCCTGATGCAGCCATGCCGCGTGTATGAAGAAGG  >  CP009273/4028855‑4029307

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