Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,608,607 C→G 17.0% I141M (ATC→ATG rimI → ribosomal‑protein‑S18‑alanine N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,608,6070CG17.0% 443.8 / 84.7 171I141M (ATC→ATGrimIribosomal‑protein‑S18‑alanine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (74/67);  new base G (14/15);  total (88/83)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

Only 100 of 187 total aligned reads displayed.

CGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCA  >  NC_000913/4608508‑4608706
                                                                                                   |                                                                                                   
cGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatc                                                                                                     >  1:1239359/1‑100 (MQ=255)
cGCTGCCGCCATTCCCCTGTACGAAAGTTTAGGCTTTAGCGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatc                                                                                                     >  2:2429673/1‑100 (MQ=255)
  cTGCCGCCATTGCCCTGTACCAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatcat                                                                                                   >  1:3714755/1‑100 (MQ=255)
  cTGCCGCAATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatcat                                                                                                   <  1:1552541/100‑1 (MQ=255)
    gccgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATgg                                                                                                 >  1:737077/1‑100 (MQ=255)
    gccgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATgg                                                                                                 <  2:509005/100‑1 (MQ=255)
     ccgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGc                                                                                                <  1:466178/100‑1 (MQ=255)
        ccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGtt                                                                                             <  2:3588564/100‑1 (MQ=255)
         cATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtct                                                                                            <  1:998410/100‑3 (MQ=255)
         cATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtct                                                                                            <  1:328987/100‑3 (MQ=255)
         cATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtct                                                                                            <  1:1325893/100‑3 (MQ=255)
         cATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTg                                                                                            <  1:341248/100‑1 (MQ=255)
          aTTGCCCTGTATGAAAGTTTAGGCTTTAACTATGTGTCGATTCGCCGCAATTATGACGCCACCGCGGACGGTCGCGAAGACGCCATGATGTTGGCGtctt                                                                                           <  2:717529/100‑4 (MQ=255)
          aTTGCCCTGTACGAAAGTTTAGGCATTAACGAGGCGACGATTCGCCGCAATTACTACACCCCCACGGACGGTCGCGAAGACGCCATGATGATGGCGtct                                                                                            <  1:1341093/99‑3 (MQ=255)
           ttGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGACATCATCATGGCGTTGcc                                                                                          <  2:2757323/100‑1 (MQ=255)
            tGCCCTGTACGAAATTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                         <  2:3178731/100‑6 (MQ=255)
            tGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                         <  2:2558562/100‑6 (MQ=255)
            tGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                         <  1:2896258/100‑6 (MQ=255)
            tGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                         <  1:3787556/100‑6 (MQ=255)
             gCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCaa                                                                                        >  1:776561/1‑100 (MQ=255)
              cccTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAAt                                                                                       <  1:3526201/100‑1 (MQ=255)
              cccTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAAt                                                                                       >  2:1880977/1‑100 (MQ=255)
                 tGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAg                                                                                    <  1:3843446/100‑1 (MQ=255)
                    aCGAAAGTTTAGGATTTAACGAGGCGACGATACGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTAt                                                                                 >  1:2659822/1‑100 (MQ=255)
                     cGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATg                                                                                >  2:3579388/1‑100 (MQ=255)
                     cGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATg                                                                                <  1:3924144/100‑1 (MQ=255)
                         aGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAAt                                                                            <  1:1245766/100‑1 (MQ=255)
                          gTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATa                                                                           <  2:2967311/100‑1 (MQ=255)
                          gTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATa                                                                           >  2:1049189/1‑100 (MQ=255)
                           tttAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATAc                                                                          >  1:2912099/1‑100 (MQ=255)
                               ggCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAgg                                                                      >  1:3016990/1‑100 (MQ=255)
                               ggCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAgg                                                                      >  1:2832169/1‑100 (MQ=255)
                                 cTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGCAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGt                                                                     >  2:1936132/1‑99 (MQ=255)
                                   ttAACGAGGCGACGATTCGCCGCAATTACTACCCCACCCCGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGa                                                                  <  1:2330988/100‑1 (MQ=255)
                                   ttAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGa                                                                  <  1:832809/100‑1 (MQ=255)
                                    tAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaa                                                                 >  1:3176335/1‑100 (MQ=255)
                                    tAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaa                                                                 >  1:3349632/1‑100 (MQ=255)
                                    tAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaa                                                                 >  1:1780751/1‑100 (MQ=255)
                                     aaCGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaat                                                                >  2:3316567/1‑100 (MQ=255)
                                       cGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaataa                                                              <  2:2642040/100‑1 (MQ=255)
                                       aattgcggcgAATCGTCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaataa                                                              >  1:973092/5‑100 (MQ=255)
                                            cgacgaTTCGCCGCAATTACTACCCCACCAGGGGCGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGAAATAATGAAg                                                         >  2:2588997/1‑100 (MQ=255)
                                            cgacgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAg                                                         >  2:2829615/1‑100 (MQ=255)
                                             gacgaTTCGCCGCAATTATTCCCGCACCACGGACGGTCGCGAAGACGCCATCAGCATGGCGTTGACAATCAGTATGTAATACAAGGTGGAATAATGAAg                                                         <  2:2214293/99‑1 (MQ=255)
                                                  ttCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGAc                                                   <  1:1955260/100‑1 (MQ=255)
                                                  ttCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGAc                                                   >  1:1026395/1‑100 (MQ=255)
                                                    cgccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTg                                                 >  1:826654/1‑100 (MQ=255)
                                                    cgccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTg                                                 >  2:778145/1‑100 (MQ=255)
                                                    cgccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTg                                                 >  1:2454557/1‑100 (MQ=255)
                                                      ccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGa                                               >  2:217027/1‑100 (MQ=255)
                                                          aaTTACTCCCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTGGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGAtttt                                           <  1:220473/100‑1 (MQ=255)
                                                            ttACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGGGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCt                                         <  1:302808/100‑1 (MQ=255)
                                                            ttACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCt                                         <  1:1919616/100‑1 (MQ=255)
                                                             tactacCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCtt                                        <  1:3773269/100‑1 (MQ=255)
                                                              actacCCCACTACGGACGGTCGCGAAGAAGCCATAAACATGAAGTAGCCAGAAAATATTAAATACAAGATGGCAGAATAAAATGGGAATGGAtatcc                                          >  2:1266362/1‑93 (MQ=255)
                                                               ctacCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTg                                      >  1:1544102/1‑100 (MQ=255)
                                                               ctacCCCACCACGGACGGTCGCGAAGACGCCATAATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGCAATAATGAAGTGGGACTGGATTTTCTTTg                                      >  2:2366184/1‑100 (MQ=255)
                                                                  ccccACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATg                                   >  1:653845/1‑100 (MQ=255)
                                                                  ccccACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATg                                   >  2:1919982/1‑100 (MQ=255)
                                                                   cccACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGc                                  <  2:1712739/100‑1 (MQ=255)
                                                                     caccacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCg                                >  1:3250825/1‑100 (MQ=255)
                                                                       ccacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGAt                              <  2:278978/100‑1 (MQ=255)
                                                                        cacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGCTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATg                             <  2:2159724/100‑1 (MQ=255)
                                                                            gacggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAAc                         <  1:2949384/100‑1 (MQ=255)
                                                                              cggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGc                       <  2:2116608/100‑1 (MQ=255)
                                                                              cggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGc                       >  1:1093037/1‑100 (MQ=255)
                                                                              cggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGc                       <  1:3849677/100‑1 (MQ=255)
                                                                               ggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCt                      >  1:2193077/1‑100 (MQ=255)
                                                                               ggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCt                      >  1:2223067/1‑100 (MQ=255)
                                                                                gTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTg                     <  2:2912185/100‑1 (MQ=255)
                                                                                 gattggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGt                    >  1:2136990/9‑100 (MQ=38)
                                                                                  cgcgAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGtt                   <  1:2715931/100‑1 (MQ=255)
                                                                                  cgcgAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGtt                   >  1:2827544/1‑100 (MQ=255)
                                                                                  cgcgAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGt                    >  1:1074650/1‑99 (MQ=255)
                                                                                   gcgAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGCCTGGATTTTCTTTGATGCCGATGAAACGCTGttt                  <  2:3143587/100‑1 (MQ=255)
                                                                                   gcgAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGttt                  >  1:2525719/1‑100 (MQ=255)
                                                                                    tggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTa                 >  2:3706768/6‑100 (MQ=255)
                                                                                    cgAAGACGCCATCATCATGGCGTTGCCAATAAGGATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAGACGCTGTTTa                 >  2:2596877/1‑100 (MQ=255)
                                                                                    cgAAGACGCCATCATCATGGCGTTGCAAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTa                 <  2:3215166/100‑1 (MQ=255)
                                                                                     ggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAc                >  2:3914635/5‑100 (MQ=255)
                                                                                     gAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAc                <  2:1455770/100‑1 (MQ=255)
                                                                                      gcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcc               >  1:556216/4‑100 (MQ=255)
                                                                                      gcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcc               >  1:1540110/4‑100 (MQ=255)
                                                                                      gcaaCGCCATGATGATGGCGTTGCCAACCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcc               >  1:2785389/4‑100 (MQ=255)
                                                                                      aaGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcc               >  2:3120786/1‑100 (MQ=255)
                                                                                      aaGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcc               >  1:2023435/1‑100 (MQ=255)
                                                                                      aaGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcc               <  1:2134828/100‑1 (MQ=255)
                                                                                       caaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCt              >  2:1170905/3‑100 (MQ=255)
                                                                                       caaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCt              >  2:3510725/3‑100 (MQ=255)
                                                                                        aaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCtt             >  1:44300/2‑100 (MQ=255)
                                                                                        aaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCtt             >  1:2641395/2‑100 (MQ=255)
                                                                                          cGCCATGATGATGGCGTTGCCAATCAGTATATAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTg           <  1:1627487/100‑1 (MQ=255)
                                                                                          cGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTg           >  1:768783/1‑100 (MQ=255)
                                                                                          cGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTg           >  1:1861446/1‑100 (MQ=255)
                                                                                            ccATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGGAGTGGGACGGGATCGTCTTTGATGCGGATGAGAAGAGGGTTGCCTTTGAc         >  2:3271179/1‑100 (MQ=255)
                                                                                            ccATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGAc         >  1:1462206/1‑100 (MQ=255)
                                                                                            ccATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGAc         <  2:610950/100‑1 (MQ=255)
                                                                                             catcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACt        <  1:1879753/100‑1 (MQ=255)
                                                                                                catcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCAt     >  2:268500/1‑100 (MQ=255)
                                                                                                 atcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCAtt    >  2:3369787/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCA  >  NC_000913/4608508‑4608706

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