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→CCGrsmF16S 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

TGCGTAAAGATCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATG  >  NC_000913/1920830‑1921279
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                             
tgCGTAAAGATCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGc                                                                                                                                                                                                                          <  2:409463/235‑1 (MQ=255)
          tCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAATAGACTTACCCCGACGCAGTAGAGATTTTACCACTTGGCGATCTCTTCCCGGGTGcaaccaa                                                                                                                                                                                                                    >  1:200928/1‑231 (MQ=255)
          tCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTTCCGGGTGcaaacaa                                                                                                                                                                                                                    >  1:76392/1‑231 (MQ=255)
          tCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGcaaacaa                                                                                                                                                                                                                    >  1:474642/1‑231 (MQ=255)
             cGATGCGCTAAAAAGCTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGa                                                                                                                                                                                                            <  1:700698/236‑1 (MQ=255)
             cGATGCGATAAAACACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACAGAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGGCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAgcg                                                                                                                                                                                                                <  2:314297/232‑1 (MQ=255)
              gATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGAc                                                                                                                                                                                                           <  2:474642/236‑1 (MQ=255)
               aTGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGAcc                                                                                                                                                                                                          >  1:264976/1‑236 (MQ=255)
                 gcgcTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGa                                                                                                                                                                                                            >  2:191223/1‑232 (MQ=255)
                 gcgcTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGa                                                                                                                                                                                                            <  1:191223/232‑1 (MQ=255)
                   gcTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGCTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGTTGACCgaag                                                                                                                                                                                                      >  1:156700/1‑236 (MQ=255)
                       gaataCTGGTCACCAGAAATCAATCAGTGAATCGCATCTACACAACGGGATCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAgg                                                                                                                                                                                                  <  2:155128/232‑1 (MQ=255)
                               gTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTCCCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGACCTGACCGAAGAAGGCttt                                                                                                                                                                                              >  1:612985/1‑232 (MQ=255)
                                    cAGAAAGCAATCAGGAACTCGCAGATACCCAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGt                                                                                                                                                                                      <  2:15107/235‑1 (MQ=255)
                                         aGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGCAGAAGGATTCTTGCCGGGTTTccc                                                                                                                                                                                >  2:751155/1‑236 (MQ=255)
                                              tCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTATGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTcccc                                                                                                                                                                               >  2:700910/1‑231 (MQ=255)
                                               cAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCcac                                                                                                                                                                             <  1:700910/232‑1 (MQ=255)
                                                 ggAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCCAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAAttt                                                                                                                                                                       >  2:642701/1‑236 (MQ=255)
                                                       cGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGa                                                                                                                                                                   <  2:594259/234‑1 (MQ=255)
                                                          aGCTGCCCAACGTGAGCTTATCGACACAGCCTTTCATGCATTCCGGCCTGTTGGTACGCTGGTTTACTCGACCTGTACCCTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTGCCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCg                                                                                                                                                              <  2:200928/236‑1 (MQ=255)
                                                                 caACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGctt                                                                                                                                                       <  1:280702/236‑1 (MQ=255)
                                                                    cGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGcttctt                                                                                                                                                    >  1:17321/1‑236 (MQ=255)
                                                                      ggAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCttcttc                                                                                                                                                   >  1:308082/1‑235 (MQ=255)
                                                                                aaaaggacttttCATTCATGACGTCCTGGTGATACGCTGGTTTACTCGACCTGGACCTTAAACCAGGAAGAAAACGAATCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTGACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCGCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCt                                                                                                                                        <  2:76392/227‑1 (MQ=255)
                                                                                   aGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCcac                                                                                                                                                                             <  1:678561/196‑1 (MQ=255)
                                                                                   aGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCcac                                                                                                                                                                             >  2:678561/1‑196 (MQ=255)
                                                                                    gcgcCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGt                                                                                                                                    >  2:614616/1‑236 (MQ=255)
                                                                                     cgcCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGt                                                                                                                                    >  1:646219/1‑235 (MQ=255)
                                                                                         tttCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTa                                                                                                                                   <  1:642701/232‑1 (MQ=255)
                                                                                                  ttACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGc                                                                                                                         >  2:155126/1‑233 (MQ=255)
                                                                                                   tACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGcaaac                                                                                                                                                                                                                      <  1:307590/140‑1 (MQ=255)
                                                                                                   tACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGcaaac                                                                                                                                                                                                                      >  2:307590/1‑140 (MQ=255)
                                                                                                    ccGTCCTGGTGGTACGCGGTTTGACTCGACCTGTACCTTAACCCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTCCCCCGACGCAGTAGAGTTTTTCCCACGTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCAAGAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGAtt                                                                                                                     <  2:706740/234‑1 (MQ=255)
                                                                                                          tggtggTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTTAGGCGATTTCCGcc                                                                                                               >  2:357837/1‑235 (MQ=255)
                                                                                                                 aCGCTGTTTTACTGAACCTGTACCTTAAACCAGGAAGAAACAGAAGCCGTTTGCCTTGGGCTGAAAGAGACTTCCCCCGACGCAGTAGATTTTTTACCACTTGGCGATCTCTTCCCGGGTACAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGCCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCg                                                                                                        <  2:308082/235‑1 (MQ=255)
                                                                                                                 aCGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCg                                                                                                        >  1:40823/1‑235 (MQ=255)
                                                                                                                    cTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCACGCCTTACCCGcccc                                                                                                    >  1:240472/1‑236 (MQ=255)
                                                                                                                     gtgTTTACTCGACTGTACCCTAAACCAAGAAAAAAACCGAGCCCTTTGGCCTGTGGCTAAAGGAGCCTTCCCCCGACACATTAGAGTTTTTACCACTTGGCTATCTCTTCCGGGGTGAAAACAAAGCGCTGACCGAAGAAGGCTTTTTTCATTTTTTCCCAAAAATTTACGCCTGCAAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGcccc                                                                                                    <  2:637278/233‑1 (MQ=255)
                                                                                                                       gTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGTCCCCa                                                                                                  >  1:637278/1‑235 (MQ=255)
                                                                                                                          tACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATa                                                                                              >  1:476577/1‑236 (MQ=255)
                                                                                                                           actcaccttGTCCCTAAACCCGGGAAGAAACCGAGCCGGTTCGCCGGTGCCTGAAGTGGCTTTCCCCCGACGAAGTGGGGTTTTTCCCACTTGGCGATCTCTCCCGGGTTGCAAAAAAAGCGCTGACGAAAGAAGGCTTTTTGCATGTTTTCCCAAAATTTTCCGACTGCAAGGGTTTCTTCGTTGCTCGTCTGCGAAAACCTCAGGCGATTCCCGCCTTACCCGCCCCCAAATAc                                                                                             <  2:17321/228‑1 (MQ=255)
                                                                                                                                 ccTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTc                                                                                       >  1:590795/1‑236 (MQ=255)
                                                                                                                                           aaaCCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTATCCGCCCCCAAATCCAAAGTCGGTAATTTTc                                                                             >  1:192894/1‑236 (MQ=255)
                                                                                                                                                      aaacCGAAGCCTTTTGCCTGGGGCTAAAAGGGCCTAACCCCGCCGCAGTGAGGTTTTCACCATTTGGCGATCTCTTCCGGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCTTGTTTTCCCAAAAATTTCCGCCTGCGAGGGCTTCTTCGTTGCTCGTCGGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGcc                                                                    <  2:192894/234‑1 (MQ=255)
                                                                                                                                                             aGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAg                                                           >  1:393142/1‑236 (MQ=255)
                                                                                                                                                                           tgGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGCAAGATCGCGAAGCCgg                                              >  1:390689/1‑236 (MQ=255)
                                                                                                                                                                                         gtACCACGACGCAGTGGAGTTTTTACCACTTGGCTATCTCTTCCCGGGTGAAAACAAGCCGTTGACCGAAGAAGGTTTTTTGCATTTTTTCCCACAAATTTCCGACTGCGAAGCCTTCTTCGTTGCTCGTCTGCGAAAAACTCAGGCGATCCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGt                                    <  2:556924/231‑1 (MQ=255)
                                                                                                                                                                                          tACCACGACTCATTAAATTTTATCCCACTTTCCTATCTCTTCCCTGCTCCAACCAAACCCATCACAAAACAATCCTTTTTGCATTTTTCCCACCAAATTTACGACTCCCAATCTTTCTTCTTTCATCCTCTTCTTAAAACTCATCCATTTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTTAAAGATCGCGAAGCTGGACAAATTCGTCAggc                               <  2:156700/236‑1 (MQ=255)
                                                                                                                                                                                                       tAGAGTTTTTACCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGgtg                    >  2:734957/1‑234 (MQ=255)
                                                                                                                                                                                                                 aCCACTTGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAAca        >  1:707433/1‑235 (MQ=255)
                                                                                                                                                                                                                    aCTTGGCTATCTCTCCCCGGTTGCAAACAAAGCCCTGACCGAAGAAGGCTTTTTCCTTTTTTCCCCAAAAATTACCGATGGCAAGGGTTTCTTCGTTCCTCGTTTGCTAAAACCTCGGGCGTTTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTggg     <  2:393142/236‑1 (MQ=255)
                                                                                                                                                                                                                      ttGGCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGAt   >  2:727468/1‑236 (MQ=255)
                                                                                                                                                                                                                        ggCGATCTCTTCCCGGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATg  >  1:219180/1‑235 (MQ=255)
                                                                                                                                                                                                                             tctcTTCCCGGGTCCAAAAAAGGCGCTGACCGAAGAAGGCTTTTTGCATGTTTCCCCAAAAATTTACGCCTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTAAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGt                                                                                        <  2:45254/144‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                             
TGCGTAAAGATCCCGATGCGCTAAAAAACTGGTCACCAGAAAGCAATCAGGAAATCGCAGCTACACAACGGGAGCTTATCGACAGCGCCTTTCATGCATTACGTCCTGGTGGTACGCTGGTTTACTCGACCTGTACCTTAAACCAGGAAGAAAACGAAGCCGTTTGCCTGTGGCTGAAAGAGACTTACCCCGACGCAGTAGAGTTTTTACCACTTGGCGATCTCTTCCCTGGTGCAAACAAAGCGCTGACCGAAGAAGGCTTTTTGCATGTTTTCCCACAAATTTACGACTGCGAAGGCTTCTTCGTTGCTCGTCTGCGTAAAACTCAGGCGATTCCCGCCTTACCCGCCCCCAAATACAAAGTCGGTAATTTTCCGTTCAGCCCGGTGAAAGATCGCGAAGCTGGACAAATTCGTCAGGCGGCTACAGGTGTTGGCTTAAACTGGGATG  >  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: