Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,867,454 2 bp→AA coding (97‑98/993 nt) rpoS ← RNA polymerase sigma factor RpoS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,867,4540TA98.6% 269.0 / NA 72Q33M (CAG→TTG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (33/38);  total (33/39)
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.98e-01
*NC_0009132,867,4550GA98.6% 267.6 / NA 72Q33M (CAG→TTG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (32/39);  minor base C (1/0);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCC  >  NC_000913/2867211‑2867697
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                  
aGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTctac                                                                                                                                                                                                                                              <  1:600763/251‑1 (MQ=255)
                         cAATTTTTACCACCAGCCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCa                                                                                                                                                                                                                     <  2:164144/251‑1 (MQ=255)
                              tttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTc                                                                                                                                                                                                                 <  1:616018/250‑1 (MQ=255)
                               ttACCACCAGACGCAAGGTACTCTCGACCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGCGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTcc                                                                                                                                                                                                                <  1:506/250‑1 (MQ=255)
                                  ccaccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGCCATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACCCGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCACTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtc                                                                                                                                                                                                            >  2:702462/1‑251 (MQ=255)
                                    cccaGACGCAAGTTACTAGCGATCATCCGGCGGCGAGAGGCGTCATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCTGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTtcgc                                                                                                                                                                                                          <  2:173115/250‑3 (MQ=255)
                                     ccaGACGCAAGTTACGCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtca                                                                                                                                                                                                         <  2:653313/251‑1 (MQ=255)
                                        gACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatc                                                                                                                                                                                                       >  1:184898/1‑250 (MQ=255)
                                         aCGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCCAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcctca                                                                                                                                                                                                      >  1:319745/1‑250 (MQ=255)
                                                  tcgccggagaatcCGGCGGCGAGAGGCGACATCTCCACGCCGTGCGCGACGCGAAAAATAACCTTCTTCTACGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCGCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCa                                                                                                                                                                                             <  2:439935/240‑1 (MQ=255)
                                                     ctcGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAATCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTTGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATc                                                                                                                                                                                           <  1:468673/249‑1 (MQ=255)
                                                      tcGATCATCCGGCGGCGCGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCNCAACTCAGTTCTCCTCAAATTCCGCATCttt                                                                                                                                                                                        >  2:331411/1‑250 (MQ=255)
                                                       cGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCtt                                                                                                                                                                                         <  1:336061/249‑1 (MQ=255)
                                                         atcatcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAtt                                                                                                                                                                                     <  2:542668/251‑1 (MQ=255)
                                                         agcatcCGGCGGCGAGAGGCGACATCTCCACNCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAtt                                                                                                                                                                                     <  2:349699/249‑1 (MQ=255)
                                                           catcCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAttt                                                                                                                                                                                    <  1:506801/250‑1 (MQ=255)
                                                               cggcggCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                                                                                                               <  1:302379/251‑1 (MQ=255)
                                                               cggcggCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATc                                                                                                                                                                               <  1:594416/251‑1 (MQ=255)
                                                                      gagaGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACt                                                                                                                                                                        <  1:50440/251‑1 (MQ=255)
                                                                             gACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCg                                                                                                                                                                 >  1:370371/1‑251 (MQ=255)
                                                                             gACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTATCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCcgcg                                                                                                                                                                 >  2:68630/1‑251 (MQ=255)
                                                                              aCATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGt                                                                                                                                                                >  1:428719/1‑251 (MQ=255)
                                                                                  ctcCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTc                                                                                                                                                            >  1:390995/1‑251 (MQ=255)
                                                                                    ccgcGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCGCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTc                                                                                                                                                            <  2:59593/246‑1 (MQ=255)
                                                                                    ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTg                                                                                                                                                          <  1:325975/251‑1 (MQ=255)
                                                                                        gCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTc                                                                                                                                                      <  2:631008/251‑1 (MQ=255)
                                                                                                 gACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAg                                                                                                                                                 >  1:471887/1‑247 (MQ=255)
                                                                                                 gACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGGTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAAGCATGAACTTTCAGCGTATTCTGACTCATAAg                                                                                                                                                 <  2:471887/247‑1 (MQ=255)
                                                                                                    gcgcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTc                                                                                                                                          >  1:461190/1‑251 (MQ=255)
                                                                                                      gcAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCt                                                                                                                                        <  1:242408/251‑1 (MQ=255)
                                                                                                             aaaCTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTg                                                                                                                                 >  1:60730/1‑251 (MQ=255)
                                                                                                                            ccGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAAc                                                                                                                  <  1:31128/251‑1 (MQ=255)
                                                                                                                              gTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCCGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACag                                                                                                                <  2:483542/251‑2 (MQ=255)
                                                                                                                              gTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACAt                                                                                                                >  2:219677/1‑251 (MQ=255)
                                                                                                                                     gttgtgAATAACCAATCGCACCAAGGTAAAGCTGAGTCGCGTACAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGGCAAAAACCTCAACTCCGTTATCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCa                                                                                                         <  2:455659/248‑1 (MQ=255)
                                                                                                                                          gAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGATCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGc                                                                                                    >  1:290308/1‑251 (MQ=255)
                                                                                                                                            aTAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCt                                                                                                  >  1:148876/1‑251 (MQ=255)
                                                                                                                                            aTAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTACTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCt                                                                                                  <  1:158942/251‑1 (MQ=255)
                                                                                                                                                ccAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGACCTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCNCTTGACGGAACATTCAAGCAAAAGCCTGGtt                                                                                              >  2:318390/1‑251 (MQ=255)
                                                                                                                                                 cAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATAACTGGGTTCCAATTCTACTCAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAACTTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTACCTCATAAGGTGTCTCCTACCCGTGATCCCTTGACGGAACCTTCAATCAAAAACCTGGTTc                                                                                             >  2:284803/1‑251 (MQ=255)
                                                                                                                                                  aaTCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTcc                                                                                            >  2:628055/1‑251 (MQ=255)
                                                                                                                                                     ctcACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGATCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTccgcc                                                                                         >  2:468969/1‑251 (MQ=255)
                                                                                                                                                       cACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTccgccg                                                                                        <  1:247584/250‑1 (MQ=255)
                                                                                                                                                       cACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGGGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTccgccg                                                                                        <  1:430005/250‑1 (MQ=255)
                                                                                                                                                       cACCAAGGTAAAGCTGAGTCGCGTCCAAAACACGCTGCGTGGCTCCGTGCGATAACAGTTACTCTTCGGCCAAATCGGTATCCCTGGGTTCCAATTCTACTAAGGCCTTTTCGACAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCCTTTAAATCATGAACTTTCAGCGAATGCTGACTCATAAGGTGGCTCCTACCCGTGAACCATTGACGGAACATTCAAGCAAACGCCTGGTTcctctg                                                                                        >  2:491310/1‑246 (MQ=255)
                                                                                                                                                            aGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTAt                                                                                  >  1:617892/1‑251 (MQ=255)
                                                                                                                                                              gTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCg                                                                                >  2:245208/1‑251 (MQ=255)
                                                                                                                                                                  aGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGc                                                                            <  2:291034/251‑1 (MQ=255)
                                                                                                                                                                     tGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATACGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGc                                                                         >  1:184156/1‑251 (MQ=255)
                                                                                                                                                                      gAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGc                                                                         >  1:290034/1‑250 (MQ=255)
                                                                                                                                                                      gAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGc                                                                         >  1:269379/1‑250 (MQ=255)
                                                                                                                                                                          cgcgTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAAt                                                                     <  1:436882/250‑1 (MQ=255)
                                                                                                                                                                            cgTCGAACTCACGCTGTGTGGCTCCCTGCGATAACTGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTATTCTACTAAGGCCTTTTCGTCAAAACCCTCAACTCCGTTCTCATCAAATTCCGCAGCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCGGCGTGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAAc                                                                  <  1:439843/251‑1 (MQ=255)
                                                                                                                                                                            cgTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAAc                                                                  >  1:83098/1‑251 (MQ=255)
                                                                                                                                                                                   aaaacGCTGTGTGGTTCCCTGCGATAATAGTNCCTCTTCGGCCATATCGTGATCAATGGGTTCCAATTCTCCTAAGGCATTTTCGTAAAAACCCTCAAATCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCAGACTCATAATGTGGCTCCTACCCGTATTCCCTTGACGGAACAATAAAGCTAAGGCCTGGATCCGCCAAAGTATCGCTGCGGCAAATAACGCAGctg                                                           <  2:320948/248‑3 (MQ=255)
                                                                                                                                                                                    acacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCggg                                                          >  1:406555/1‑251 (MQ=255)
                                                                                                                                                                                           gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTAc                                                     <  1:628055/249‑1 (MQ=255)
                                                                                                                                                                                              tgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGAtt                                                >  1:107738/1‑251 (MQ=255)
                                                                                                                                                                                                               aGTTCCTCTTCGGCCAATTCGTTATCACTGGGGCCCAATTCTCATAATCCCTTTGCTGCAAAACTCCCAACACCTTGCCCATCAAAGTCCGCATCTTCATATAAATCATGATCTTTCAGCGTATTCTCACTCAAAAGGTGGCTTCTACCCTTGATCACGCGACAGTACATTCAATCAAAAGCCTTTTTCCGCCGATTTATCGCGGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGGAACGAAtt                               <  2:454834/251‑1 (MQ=255)
                                                                                                                                                                                                                    ctctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAACTCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCaa                           >  1:451854/1‑250 (MQ=255)
                                                                                                                                                                                                                     tctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTACTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCaaa                          >  2:35973/1‑250 (MQ=255)
                                                                                                                                                                                                                      ctTCGGCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGGGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCaaa                          <  1:331411/249‑1 (MQ=255)
                                                                                                                                                                                                                         cGGCCAAAACGTTATCACGGGGTTCCAATTCTACTAAGGCCGTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCTCATCTGCATTTAAATCATGAACTTTCAGCTTATTCTGACTCAAAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGc                      <  2:83098/250‑1 (MQ=255)
                                                                                                                                                                                                                          ggCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCaa                    <  1:335272/251‑1 (MQ=255)
                                                                                                                                                                                                                          ggCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCaa                    <  1:219677/251‑1 (MQ=255)
                                                                                                                                                                                                                           gCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCCGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCa                     >  1:422881/1‑249 (MQ=255)
                                                                                                                                                                                                                           gCCAAATCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAg                   <  1:361437/251‑1 (MQ=255)
                                                                                                                                                                                                                                aTCGTTATCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCgtgt              <  1:352796/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aTCACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACt        <  1:537489/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cACTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTgg      >  1:143798/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aCTGGGTTCCAATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATCTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTGATCGCTGCGGCCAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAACTGCAAGCGTGTTGTACTGGt     >  1:327746/1‑251 (MQ=255)
                                                                                                                                                                                                                                            gggTTCCAATTCTACTACGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTcc  >  2:303490/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                  
AGACCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGGCGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCCTGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGTTGAACTGGTTCC  >  NC_000913/2867211‑2867697

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