Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,921,059 T→G P279P (CCT→CCG)  rsmF → 16S rRNA m(5)C1407 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,921,0590TG98.2% 177.9 / NA 54P279P (CCT→CCG) rsmF16S rRNA m(5)C1407 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (30/23);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGTAAAGATCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATG  >  NC_000913/1920830‑1921279
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                             
tgCGTAAAGATCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGc                                                                                                                                                                                                                          <  2:409463/235‑1 (MQ=255)
          tCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAATAGACTTACCCCGACGCAGTAGAGATTTTACCACTTGGCGATCTCTTCCCGGGTGcaaccaa                                                                                                                                                                                                                    >  1:200928/1‑231 (MQ=255)
          tCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTTCCGGGTGcaaacaa                                                                                                                                                                                                                    >  1:76392/1‑231 (MQ=255)
          tCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGcaaacaa                                                                                                                                                                                                                    >  1:474642/1‑231 (MQ=255)
             cGATGCGCTAAAAAGCTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGa                                                                                                                                                                                                            <  1:700698/236‑1 (MQ=255)
             cGATGCGATAAAACACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACAGAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGGCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAgcg                                                                                                                                                                                                                <  2:314297/232‑1 (MQ=255)
              gATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGAc                                                                                                                                                                                                           <  2:474642/236‑1 (MQ=255)
               aTGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGAcc                                                                                                                                                                                                          >  1:264976/1‑236 (MQ=255)
                 gcgcTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGa                                                                                                                                                                                                            >  2:191223/1‑232 (MQ=255)
                 gcgcTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGa                                                                                                                                                                                                            <  1:191223/232‑1 (MQ=255)
                   gcTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGCTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGTTGACCgaag                                                                                                                                                                                                      >  1:156700/1‑236 (MQ=255)
                       gaataCTGGTCACCAGAAATCAATCAGTGAATCGCATCTACACAACGGGATCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAgg                                                                                                                                                                                                  <  2:155128/232‑1 (MQ=255)
                               gTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTCCCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGACCTGACCGAAGAAGGCttt                                                                                                                                                                                              >  1:612985/1‑232 (MQ=255)
                                    cAGAAAGCAATCAGGAACTCGCAGATACCCAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGt                                                                                                                                                                                      <  2:15107/235‑1 (MQ=255)
                                         aGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGCAGAAGGATTCTTGCCGGGTTTccc                                                                                                                                                                                >  2:751155/1‑236 (MQ=255)
                                              tCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTATGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTcccc                                                                                                                                                                               >  2:700910/1‑231 (MQ=255)
                                               cAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCca‑c                                                                                                                                                                             <  1:700910/232‑1 (MQ=255)
                                                 ggAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCCAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAAttt                                                                                                                                                                       >  2:642701/1‑236 (MQ=255)
                                                       cGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGa                                                                                                                                                                   <  2:594259/234‑1 (MQ=255)
                                                          aGCTGCCCAACGTGAGCTTATCGACACAGCCTTTCATGCATTCCGGCCTGTTGGTACGCTGGTTTACTCGACCTGTACCCTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTGCCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCg                                                                                                                                                              <  2:200928/236‑1 (MQ=255)
                                                                 caACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGctt                                                                                                                                                       <  1:280702/236‑1 (MQ=255)
                                                                    cGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGcttctt                                                                                                                                                    >  1:17321/1‑236 (MQ=255)
                                                                      ggAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCttcttc                                                                                                                                                   >  1:308082/1‑235 (MQ=255)
                                                                                aaaaggacttttCATTCATGACGTCCTGGTGATACGCTGGTTTACTCGACCTGGACCTTAAACCAGGAAGAAAACGAATCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTGACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCGCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCt                                                                                                                                        <  2:76392/227‑1 (MQ=255)
                                                                                   aGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCca‑c                                                                                                                                                                             <  1:678561/196‑1 (MQ=255)
                                                                                   aGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCca‑c                                                                                                                                                                             >  2:678561/1‑196 (MQ=255)
                                                                                    gcgcCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGt                                                                                                                                    >  2:614616/1‑236 (MQ=255)
                                                                                     cgcCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGt                                                                                                                                    >  1:646219/1‑235 (MQ=255)
                                                                                         tttCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTa                                                                                                                                   <  1:642701/232‑1 (MQ=255)
                                                                                                  ttACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGc                                                                                                                         >  2:155126/1‑233 (MQ=255)
                                                                                                   tACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGcaaac                                                                                                                                                                                                                      <  1:307590/140‑1 (MQ=255)
                                                                                                   tACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGcaaac                                                                                                                                                                                                                      >  2:307590/1‑140 (MQ=255)
                                                                                                    ccGTCCTGGTGGTACGCGGTTTGACTCGACCTGTACCTTAACCCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTCCCCCGACGCAGTAGAGTTTTTCCCACGTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑AGAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGAtt                                                                                                                     <  2:706740/234‑1 (MQ=255)
                                                                                                          tggtggTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTTAGGCGATTTCCGcc                                                                                                               >  2:357837/1‑235 (MQ=255)
                                                                                                                 aCGCTGTTTTACTGAACCTGTACCTTAAACCAGGAAGAAACAGAAGCCGTTT‑GCCTTGGGCTGAAAGAGACTTCCCCCGACGCAGTAGATTTTTTACCACTTGGCGATCTCTTCCCGGGTACAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGCCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCg                                                                                                        <  2:308082/235‑1 (MQ=255)
                                                                                                                 aCGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCg                                                                                                        >  1:40823/1‑235 (MQ=255)
                                                                                                                    cTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCACGCCTTACCCGcccc                                                                                                    >  1:240472/1‑236 (MQ=255)
                                                                                                                     gtgTTTACTCGA‑CTGTACCCTAAACCAAGAAAAAAACCGAGCCCTTTGGCCTGTGGCTAAAGGAGCCTTCCCCCGACACATTAGAGTTTTTACCACTTGGCTATCTCTTCCGGGGTGAAAACAAAGCGCTGACCGAAGAAGGCTTTTTTCATTTTTTCCCA‑AAAATTTACGCCTGCAAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGcccc                                                                                                    <  2:637278/233‑1 (MQ=255)
                                                                                                                       gTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGTCCCCa                                                                                                  >  1:637278/1‑235 (MQ=255)
                                                                                                                          tACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATa                                                                                              >  1:476577/1‑236 (MQ=255)
                                                                                                                           actcaccttGTCCCTAAACCCGGGAAGAAACCGAGCCGGTTC‑GCCGGTGCCTGAAGTGGCTTTCCCCCGACGAAGTGGGGTTTTTCCCACTTGGCGATCTCTCCCGGGTTGCAAAAAAAGCGCTGACGAAAGAAGGCTTTTTGCATGTTTTCCCA‑AAATTTTCCGACTGCAAGGGTTTCTTCGTTGCTCGTCTGCGAAAACCTCAGGCGATTCCCGCCTTACCCGCCCCCAAATAc                                                                                             <  2:17321/228‑1 (MQ=255)
                                                                                                                                 ccTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTc                                                                                       >  1:590795/1‑236 (MQ=255)
                                                                                                                                           aaaCCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTATCCGCCCCCAAATCCAAAGTCGGTAATTTTc                                                                             >  1:192894/1‑236 (MQ=255)
                                                                                                                                                      aaacCGAAGCCTTTT‑GCCTGGGGCTAAAAGGGCCTAACCCCGCCGCAGTGAGGTTTTCACCATTTGGCGATCTCTTCCGGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCTTGTTTTCCCA‑AAAATTTCCGCCTGCGAGGGCTTCTTCGTTGCTCGTCGGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGcc                                                                    <  2:192894/234‑1 (MQ=255)
                                                                                                                                                             aGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAg                                                           >  1:393142/1‑236 (MQ=255)
                                                                                                                                                                           tgGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGCAAGATCGCGAAGCCgg                                              >  1:390689/1‑236 (MQ=255)
                                                                                                                                                                                         gtACCACGACGCAGTGGAGTTTTTACCACTTGGCTATCTCTTCCCGGGTGAAAACAAGCCGTTGACCGAAGAAGGTTTTTTGCATTTTTTCCCA‑CAAATTTCCGACTGCGAAGCCTTCTTCGTTGCTCGTCTGCGAAAAACTCAGGCGATCCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGt                                    <  2:556924/231‑1 (MQ=255)
                                                                                                                                                                                          tACCACGACTCATTAAATTTTATCCCACTTTCCTATCTCTTCCCTGCTCCAACCAAACCCATCACAAAACAATCCTTTTTGCAT‑TTTTCCCACCAAATTTACGACTCCCAATCTTTCTTCTTTCATCCTCTTCTTAAAACTCATCCATTTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTTAAAGATCGCGAAGCTGGACAAATTCGTCAggc                               <  2:156700/236‑1 (MQ=255)
                                                                                                                                                                                                       tAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGgtg                    >  2:734957/1‑234 (MQ=255)
                                                                                                                                                                                                                 aCCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAAca        >  1:707433/1‑235 (MQ=255)
                                                                                                                                                                                                                    aCTTGGCTATCTCTCCCCGGTTGCAAACAAAGCCCTGACCGAAGAAGGCTTTTTCCTTTTTTCCCCA‑AAAATTACCGATGGCAAGGGTTTCTTCGTTCCTCGTTTGCTAAAACCTCGGGCGTTTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTggg     <  2:393142/236‑1 (MQ=255)
                                                                                                                                                                                                                      ttGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGAt   >  2:727468/1‑236 (MQ=255)
                                                                                                                                                                                                                        ggCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATg  >  1:219180/1‑235 (MQ=255)
                                                                                                                                                                                                                             tctcTTCCCGGGTCCAAAAAAGGCGCTGACCGAAGAAGGCTTTTTGCATGTTTCCCCA‑AAAATTTACGCCTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTAAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGt                                                                                        <  2:45254/144‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                             
TGCGTAAAGATCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTT‑GCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCA‑CAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATG  >  NC_000913/1920830‑1921279

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