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,4540TA97.6% 124.0 / NA 43Q33M (CAG→TTG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base A (15/26);  minor base C (1/0);  total (16/27)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009132,867,4550GA97.7% 132.4 / NA 43Q33M (CAG→TTG) rpoSRNA polymerase sigma factor RpoS
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (16/26);  total (16/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑TGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGT  >  NC_000913/2867214‑2867685
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                      
ccACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCGCGATCATCCGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTaa                                                                                                                                                                                                                               <  1:509801/251‑1 (MQ=255)
 cACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCTGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAg                                                                                                                                                                                                                              <  1:475219/251‑1 (MQ=255)
 cACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAg                                                                                                                                                                                                                              <  1:136979/251‑1 (MQ=255)
 cACGATTGCCATAACGGCGGACAAAATTTACCAACACACACACGATAGTCCCGGTCATCCGGAGGGGA‑AGGTCCCCAACTCCCAGCAGAGCGCGACCAGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGGATAACCAATCTCACCAAAGTACAGCAGAGGCGCGTCCAAGCCACGCTGTGGGGAGCCCCGCGATAACAAGTCCTCTTCGGCCAAATCGTTATAAATGGGTTCC‑CATTCTACTaa                                                                                                                                                                                                                               >  2:488080/1‑250 (MQ=255)
     aTTGCCATAACGGCGGGCAATTTTTAACACCAGACGCAAGTTAATCTCGATCATCCGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATATCACCAAGGGAAAGCTGAGTCGCGTCCAACACACGCTGTGGGGCGCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGGTCC‑AATTCTACTAAGGCCt                                                                                                                                                                                                                          >  1:178927/1‑251 (MQ=255)
         ccATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGG‑CGACATCGCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTGCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTc                                                                                                                                                                                                                      <  1:471477/251‑1 (MQ=255)
          cATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCg                                                                                                                                                                                                                     <  1:93460/251‑1 (MQ=255)
                  cggGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCa                                                                                                                                                                                                                  <  1:322200/246‑1 (MQ=255)
                  cggGCAATTTTTACAACCAGACGAAAGTTACTCTCGATCATCCGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCa                                                                                                                                                                                                                  >  2:322200/1‑246 (MQ=255)
                           tttACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCTCGACGCGCAAACTAAACTTCTTGTTCGGCCGTTAACAGAGGTGAATCACCAATCTCACCAAGGTAACGCTGAGACGAGTCCAACACACGCTGTGTGACTCCCTGCGATAACAGTTCCTCTTCGGCCCAATCGTTCTCACTGGTTTCC‑AATTCTACTAAGGCCTTTTCGTCTACAACCTCAACtac                                                                                                                                                                                                    >  1:228986/1‑249 (MQ=255)
                            ttgccaccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTcc                                                                                                                                                                                                    <  1:3088/247‑1 (MQ=255)
                              accaccAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGG‑CGACATCGCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGtt                                                                                                                                                                                                 <  1:21174/251‑1 (MQ=255)
                              accaccAGACGCAAGGTACTCTCGATCATCCGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCATTTTCGTCAAACACCTCAACTCCGtt                                                                                                                                                                                                 >  2:359162/1‑251 (MQ=255)
                                 cccaGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAATGCCTTTTCGTCAAAAACCTCAACTCCGTtctc                                                                                                                                                                                              <  2:253887/250‑1 (MQ=255)
                                    aGACGCAAGTTCCTCTAGATCATCAGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCATCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGACCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCtcatc                                                                                                                                                                                           >  2:62902/1‑251 (MQ=255)
                                            gTTACTCTCGATCATCCGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCg                                                                                                                                                                                   <  2:235985/251‑1 (MQ=255)
                                                               cggcgAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCGGCGATAACAGTTCCTCTTCGGCCAAATCGTTAGCACTTGGTTCC‑AATTCGACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATg                                                                                                                                                                <  1:487122/251‑1 (MQ=255)
                                                                  cgAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCGGTGGGGCTCCCTGCGATAACAGTTCCTCGTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCAGCAAATTCCGCATCTTCATTTAAAGCATGAAc                                                                                                                                                             <  2:27391/251‑1 (MQ=255)
                                                                     gagG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACttt                                                                                                                                                          >  1:283951/1‑251 (MQ=255)
                                                                        g‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAg                                                                                                                                                       <  1:138498/251‑1 (MQ=255)
                                                                                  ccACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTg                                                                                                                                              <  1:73422/251‑1 (MQ=255)
                                                                                                       aaaTAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTGCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCGACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTa                                                                                                                                                   <  1:504537/225‑1 (MQ=255)
                                                                                                       aaaTAAACTTCTTCTTCCGCCGTTAACCGTGGAGACTCACCAATATCACCACGGTAAAGCTGAGTCGAGCCCAACACACGCTGTGTGGCTCACTGCGATAACAGTTCCTCTTAGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTCAATCATGAACTTTCAGCGTa                                                                                                                                                   >  2:504537/1‑225 (MQ=255)
                                                                                                                       tggCCGTTAACAGGGGTGAATAAGCAATCTCCCCAAGGTAAAGCTGAGTAACGTCCACCACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAACTTCCGCATCTTCCTTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGCTATCCCTTGGCgg                                                                                                         >  2:75730/2‑251 (MQ=255)
                                                                                                                             ttAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGCACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACAtt                                                                                                   >  1:191843/1‑251 (MQ=255)
                                                                                                                                  aGTGGTGAATAAACAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCACTGCGATAACAGCTCCTCTTAGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATAATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGc                                                                                              >  2:25160/1‑251 (MQ=255)
                                                                                                                                     ggtgAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCaaa                                                                                           <  1:114613/251‑1 (MQ=255)
                                                                                                                                           tAACCACTCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATCACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTCAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCCTCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTg                                                                                     >  2:74289/1‑251 (MQ=255)
                                                                                                                                              ccAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGtt                                                                                  <  2:121771/251‑1 (MQ=255)
                                                                                                                                                           ggTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTGCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATc                                                                     <  1:138741/251‑1 (MQ=255)
                                                                                                                                                                           gTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACg                                                     <  1:36946/251‑1 (MQ=255)
                                                                                                                                                                            tCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTGCCTCTTCGGCCAAATCGTTATCACTGGGTGCC‑AATTCTACTAAGGCCTTTTCTTCAACAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTGCCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGTTTCCGCCGATTTATCGCTGCGTCATATAACg                                                     <  2:162520/250‑1 (MQ=255)
                                                                                                                                                                                   cacGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCGTCATTTAAATCATGAACTTTCAGCGTATTCTGACGCATAAGGGGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATGTATCGCTGCGGCAAATAACGCAGCGGGt                                             <  2:69045/251‑1 (MQ=255)
                                                                                                                                                                                       atgtgtgGCTCCCGGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCACCTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCCTGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGACTGGTTCCGCCGATTTCTCGCTGCGGCAAATAACGAAGCGGGTTTAc                                         >  2:26683/2‑251 (MQ=255)
                                                                                                                                                                                         gtgtgGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCCGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCCGCGGGTTTACgg                                       >  1:23195/1‑251 (MQ=255)
                                                                                                                                                                                                          aaCAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAAAGCAGCGGGTTTACGGATTTCCCCTTGTAACg                       >  1:349380/1‑250 (MQ=255)
                                                                                                                                                                                                           aCAGTTCCTCTTCGGCCAAATCTTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGaa                     <  1:279798/251‑1 (MQ=255)
                                                                                                                                                                                                            cAGTTCCTCTTCGGCCAAATCGTGATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCAGTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTACCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAAt                    <  2:487636/251‑1 (MQ=255)
                                                                                                                                                                                                                      tCGGCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAAt            <  2:438369/249‑1 (MQ=255)
                                                                                                                                                                                                                      ccGGCCAAATCGTTATCACGGGGTTCCATGT‑CTACTAAGGCCTTTTCGTCAAAAACCTCAACTGCTTTCTCACCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAATGTGGCTCCTTCCCGTGAACCCTTGACGTAACATGCAAGCAAAAGCCTGGTTCCGCCTATTTATCGCTGCGGCAGATAACGGAGCGTGTTTACGGAGTTCCCCTTGTAACGAAGTTCAAAATGc          <  2:203267/250‑1 (MQ=255)
                                                                                                                                                                                                                        ggCCAAATCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCaa        <  1:88747/251‑1 (MQ=255)
                                                                                                                                                                                                                              aTCGTTATCACTGTGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCTCATCTTCATTTAAATCATGAACTTACAGTGTAGTCTGACTCATAAGGTTGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCGGGTTCCGCCGTTTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTATCGAATTTCAAAATGCAAGCgtgt  <  2:262357/251‑1 (MQ=255)
                                                                                                                                                                                                                              aTCGTTATCACTGGGTTCC‑AATTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAAAGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCgtgt  >  2:208236/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                      
CCACGATTGCCATAACGGCGGGCAATTTTTACCACCAGACGCAAGTTACTCTCGATCATCCGGCGGCGAGAGG‑CGACATCTCCACGCAGTGCGCGACGCGCAAAATAAACTTCTTCTTCGGCCGTTAACAGTGGTGAATAACCAATCTCACCAAGGTAAAGCTGAGTCGCGTCCAACACACGCTGTGTGGCTCCCTGCGATAACAGTTCCTCTTCGGCCAAATCGTTATCACTGGGTTCC‑TGTTCTACTAAGGCCTTTTCGTCAAAAACCTCAACTCCGTTCTCATCAAATTCCGCATCTTCATTTAAATCATGAACTTTCAGCGTATTCTGACTCATAAGGTGGCTCCTACCCGTGATCCCTTGACGGAACATTCAAGCAAAAGCCTGGTTCCGCCGATTTATCGCTGCGGCAAATAACGCAGCGGGTTTACGGATTTCCCCTTGTAACGAATTTCAAAATGCAAGCGTGT  >  NC_000913/2867214‑2867685

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