Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,194,106 C→T L21L (CTG→CTA)  mrp ← P‑loop NTPase family protein Mrp

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,194,1060CT97.7% 146.4 / NA 43L21L (CTG→CTA) mrpP‑loop NTPase family protein Mrp
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/20);  minor base G (0/1);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATTCCCGGCTGATTCTTCACGCGTTTCAGAGTGGCTATGTTATGCGACAGTTTCCAGTCGATAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCCAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGT  >  NC_000913/2193859‑2194345
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
tAATTCCCGGCTGATTCTTCACGCGTTTCAGAGTGGCTATGTTATGCGACAGTTTCCAGTCGATAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAgt                                                                                                                                                                                                                                              >  2:448563/1‑251 (MQ=255)
       cGGCTGATTCTTCACGCGTTTCAGAGTGGCTATGTTATGCGACAGTTTCCAGTCGATAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCgg                                                                                                                                                                                                                                       <  2:235243/251‑1 (MQ=255)
               tcttcACGCGTTTCAGATTGGCTATGTTATGCGACAGTTTCCAGTCGATAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTTAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCAt                                                                                                                                                                                                                                <  1:294911/250‑1 (MQ=255)
                        atttCAGAGGGGCTATGTTAGGCGACAGATTCCAGTCGATAGCTTTCGCCCCGGTGATGCGAAGAAGTTCGGCGCTACATTGCTCTTTTAACTCTTAGAACGAGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAAAGTGTCGTCCATCCAGGCGACATGGTGCAACGCATTCAGAGTGGAGAGGTTATGCTTCAGGGTTGGGTGCTGAAAAGTGGCTAGGGTCCCGGGGACAATTGGGCGCa                                                                                                                                                                                                                        >  2:236018/2‑249 (MQ=255)
                               agTGGCTATGTTATGCGACAGTTTCCAGTCGATAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTg                                                                                                                                                                                                               <  1:42085/251‑1 (MQ=255)
                                   gCTATGTTATGCGACAGTTTCCAGTCGATAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTg                                                                                                                                                                                                            >  2:182741/1‑250 (MQ=255)
                                            tGCGACAGTTTACAGGCGATAGCTTTCGCCCCGGCGGTGCATAGGAGTTCGGCGCTACATACCTCTTTTGACTCGTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTCCAGCGTGGTGAGGTTCTGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGGCCGGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTg                                                                                                                                                                                                  >  2:195553/1‑251 (MQ=255)
                                                   gTTTCCAGTCGATAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTgt                                                                                                                                                                                           <  1:239236/251‑1 (MQ=255)
                                                    tttCCAGTCGATAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTgtt                                                                                                                                                                                          >  1:407298/1‑251 (MQ=255)
                                                         cgTNGATAGCTTTCGCCCCGGGTATGCGCAGCAGTTCGGCGCTACATTGCTCGTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGATGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTgttcgttc                                                                                                                                                                                     <  2:411639/250‑1 (MQ=255)
                                                                aGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGa                                                                                                                                                                              <  1:448563/251‑1 (MQ=255)
                                                                      cGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATGCAGGCGACATGGTGCAACGATTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGATGAAAATTGGCTAGTGTCCCGGCGGCCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTccg                                                                                                                                                                         <  1:263810/250‑2 (MQ=255)
                                                                          ccGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCtttttt                                                                                                                                                                    <  1:222435/251‑1 (MQ=255)
                                                                            ggTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGt                                                                                                                                                                  >  2:69702/1‑251 (MQ=255)
                                                                                               gcgcTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAGAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTGCTGTTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTAGAGGTAAATCTTCAGCTg                                                                                                                                               <  2:62979/251‑1 (MQ=255)
                                                                                                      aTTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTcc                                                                                                                                                                          <  1:288433/217‑1 (MQ=255)
                                                                                                      aTTGCTATTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTcc                                                                                                                                                                          >  2:288433/1‑217 (MQ=255)
                                                                                                            ctTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCCTCCAGGAGACATGGTGCAACGCTTTCAGCGTGCTGAGGTTATGCTTCAGGGTTGGGTCCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGCCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTg                                                                                                                                   >  1:55932/1‑250 (MQ=255)
                                                                                                                 aaCTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAat                                                                                                                             <  2:396614/251‑1 (MQ=255)
                                                                                                                     ctTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAatat                                                                                                                           <  1:162952/249‑1 (MQ=255)
                                                                                                                       tCGAACGGGCTAAGCCAGACAAACGGCATAACCAGTTCCACATGCAAAGTGTCGGCCATCCAGGCGACATGGTGCAAAGCTTTCAGCGTGGTGAGGTTATTCTTCAAGGTTGGGTGCTGAAAATTGGCGAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTTGTGATTTGGGCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCAGCAGCTGATGGTTTCACTGAAGAAtata                                                                                                                          <  2:182070/248‑1 (MQ=255)
                                                                                                                         gAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCt                                                                                                                                                <  1:193642/224‑1 (MQ=255)
                                                                                                                         gAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCt                                                                                                                                                >  2:193642/1‑224 (MQ=255)
                                                                                                                               gcTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGtactt                                                                                                               >  1:9773/1‑251 (MQ=255)
                                                                                                                                           aaaCGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGAGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGCTAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACAtttt                                                                                                   >  1:211949/1‑251 (MQ=255)
                                                                                                                                             aCGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTcc                                                                                                 <  1:20827/251‑1 (MQ=255)
                                                                                                                                                   aaaCCAGCTCCCCATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCATTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTg                                                                                           <  1:104322/250‑1 (MQ=255)
                                                                                                                                                       cAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTAc                                                                                       >  1:321993/1‑251 (MQ=255)
                                                                                                                                                        aGTTCCACATGCAACGTGTCGTCCCTCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTAAGCTTCAGCGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCCTTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGCTTGTTCGTTCATCCCGACTCCTTTTTTGCATAGAGAAACCATCATCTGATAGTTTTCCTGAAGAATATAGAGAAGTACATACTTAACATTTTCCCATTTGGTACt                                                                                      <  1:442050/251‑1 (MQ=255)
                                                                                                                                                                aTGCAACGGGTCGTACATCAAGGCGACATGGTGCAACGCTTTAAGCGTGGTGCGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGAAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAAcc                                                                              >  2:45608/1‑251 (MQ=255)
                                                                                                                                                                   cAACGTGTCGTCCATCCAGGCGACCTGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATATATAAACCATCAGCTGATAGTTTACCTGACGAATATAGAGAAGTACTTACTTAACATTTTACCATTTGGTACTATCTAACCCCt                                                                           >  2:275681/1‑251 (MQ=255)
                                                                                                                                                                           cgtcCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTa                                                                   <  1:144330/251‑1 (MQ=255)
                                                                                                                                                                                         aCATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGc                                                     >  1:232285/1‑251 (MQ=255)
                                                                                                                                                                                           aTGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCct                                                   >  1:400309/1‑251 (MQ=255)
                                                                                                                                                                                                gcaccGCTTTCAGCTTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATACATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTTCTTACTTAACATTTTCCCATTTGGTCCTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTAt                                              <  2:399695/247‑1 (MQ=255)
                                                                                                                                                                                                     gCTTTCAGCGTGGTGAGGTTATGCTTCCGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTc                                         >  1:397934/1‑251 (MQ=255)
                                                                                                                                                                                                                      gTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAAtt                        >  1:48420/1‑251 (MQ=255)
                                                                                                                                                                                                                                aGGGTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGt              >  2:58679/1‑251 (MQ=255)
                                                                                                                                                                                                                                 gggTTGGGTGCTGAAAATTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTg             >  2:213618/1‑251 (MQ=255)
                                                                                                                                                                                                                                              aaaaTTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTa                                                                   <  2:318579/184‑1 (MQ=255)
                                                                                                                                                                                                                                              aaaaTTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTa                                                                   >  1:318579/1‑184 (MQ=255)
                                                                                                                                                                                                                                              aaaaTTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGt  >  1:53810/1‑249 (MQ=255)
                                                                                                                                                                                                                                              aaaaTTGGCTAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGt  <  2:53810/249‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               
TAATTCCCGGCTGATTCTTCACGCGTTTCAGAGTGGCTATGTTATGCGACAGTTTCCAGTCGATAGCTTTCGCCCCGGTGATGCGCAGCAGTTCGGCGCTACATTGCTCTTTTAACTCTTCGAACGCGCTATGCCAGACAAACGGCATAACCAGTTCCACATGCAACGTGTCGTCCATCCAGGCGACATGGTGCAACGCTTTCAGCGTGGTGAGGTTATGCTTCAGGGTTGGGTGCTGAAAATTGGCCAGTGTCCCGGCGACCATTGCGCGCAAGGCTTCTGGTGATTTGGCCTGGGATTGTTCGTTCATCCCGACTCCTTTTTTGTATAGATAAACCATCAGCTGATAGTTTACCTGAAGAATATAGAGAAGTACTTACTTAACATTTTCCCATTTGGTACTATCTAACCCCTTTTCACTATTAAGAAGTAATGCCTACTATGACTCAAGTCGCGAAGAAAATTCTGGTGACGTGCGCACTGCCGT  >  NC_000913/2193859‑2194345

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