Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,608,604 C→G 17.1% I140M (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,6040CG17.1% 431.9 / 85.0 171I140M (ATC→ATGrimIribosomal‑protein‑S18‑alanine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (72/68);  new base G (14/15);  total (88/83)
Fisher's exact test for biased strand distribution p-value = 8.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.02e-01

Only 100 of 185 total aligned reads displayed.

AAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCAT  >  NC_000913/4608505‑4608703
                                                                                                   |                                                                                                   
aaaCGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatc                                                                                                     <  2:1847307/100‑1 (MQ=255)
 aaCGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatca                                                                                                    >  2:1504038/1‑100 (MQ=255)
 aaCGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatca                                                                                                    <  2:2479155/100‑1 (MQ=255)
 aaCGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatca                                                                                                    <  2:1638263/100‑1 (MQ=255)
   cGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatc                                                                                                  >  1:1239359/1‑100 (MQ=255)
   cGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatc                                                                                                  >  1:3917235/1‑100 (MQ=255)
   cGCTGCCGCCATTCCCCTGTACGAAAGTTTAGGCTTTAGCGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatc                                                                                                  >  2:2429673/1‑100 (MQ=255)
     cTGCCGCCATTGCCCTGTACCAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatcat                                                                                                >  1:3714755/1‑100 (MQ=255)
       gccgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATgg                                                                                              >  1:737077/1‑100 (MQ=255)
       gccgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATgg                                                                                              <  2:509005/100‑1 (MQ=255)
        ccgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGc                                                                                             <  1:466178/100‑1 (MQ=255)
         cgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCGACCACGGACGGTCGCGAAGACGCCATCATCATGGCg                                                                                            >  1:2770142/1‑100 (MQ=255)
           ccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGtt                                                                                          <  1:2833820/100‑1 (MQ=255)
            cATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtct                                                                                         <  1:1325893/100‑3 (MQ=255)
            cATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtct                                                                                         <  1:998410/100‑3 (MQ=255)
             aTTGCCCTGTATGAAAGTTTAGGCTTTAACTATGTGTCGATTCGCCGCAATTATGACGCCACCGCGGACGGTCGCGAAGACGCCATGATGTTGGCGtctt                                                                                        <  2:717529/100‑4 (MQ=255)
             aTTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCAACCACGGACGGTCGCGAAGACGCCATCATCACGGCGCTGc                                                                                        >  2:3626705/1‑100 (MQ=255)
             aTTGCCCTGTACGAAAGTTTAGGCATTAACGAGGCGACGATTCGCCGCAATTACTACACCCCCACGGACGGTCGCGAAGACGCCATGATGATGGCGtct                                                                                         <  1:1341093/99‑3 (MQ=255)
               tGCCCTGTACGAAATTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                      <  2:3178731/100‑6 (MQ=255)
               tGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                      <  1:3787556/100‑6 (MQ=255)
               tGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                      <  2:2558562/100‑6 (MQ=255)
               tGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                      <  1:2958991/100‑6 (MQ=255)
               tGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCACACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                      <  2:1932788/100‑6 (MQ=255)
                gCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCaa                                                                                     >  1:776561/1‑100 (MQ=255)
                 cccTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAAt                                                                                    >  2:1880977/1‑100 (MQ=255)
                 cccTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAAt                                                                                    <  1:3526201/100‑1 (MQ=255)
                  ccTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATc                                                                                   >  1:1580977/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:1551849/100‑1 (MQ=255)
                        cGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATg                                                                             <  1:3924144/100‑1 (MQ=255)
                          aaaGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCCTGGCGTTGCCAATCAGTATGTa                                                                           <  2:2708130/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:2912203/1‑100 (MQ=255)
                              tttAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATAc                                                                       >  1:3341874/1‑100 (MQ=255)
                                  ggCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAgg                                                                   >  1:3016990/1‑100 (MQ=255)
                                      ttAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGa                                                               <  1:832809/100‑1 (MQ=255)
                                       tAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaa                                                              >  1:3176335/1‑100 (MQ=255)
                                        aaCGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaat                                                             >  1:2211700/1‑100 (MQ=255)
                                        aaCGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaat                                                             >  1:1734642/1‑100 (MQ=255)
                                          cGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGTCGGTCGCGAAGACGCCATCTTCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaataa                                                           <  2:2151749/100‑1 (MQ=255)
                                          aattgcggcgAATCGTCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaataa                                                           >  1:973092/5‑100 (MQ=255)
                                              gCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGaa                                                       >  1:2909057/1‑100 (MQ=255)
                                               cgacgaTTCGCCGCAATTACTACCCCACCAGGGGCGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGAAATAATGAAg                                                      >  2:2588997/1‑100 (MQ=255)
                                               cgacgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAg                                                      >  2:2829615/1‑100 (MQ=255)
                                                gacgaTTCGCCGCAATTATTCCCGCACCACGGACGGTCGCGAAGACGCCATCAGCATGGCGTTGACAATCAGTATGTAATACAAGGTGGAATAATGAAg                                                      <  2:2214293/99‑1 (MQ=255)
                                                  cgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGTAATAATGAAGttg                                                   <  1:2681784/100‑3 (MQ=255)
                                                  cgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTgg                                                   <  1:3418329/100‑1 (MQ=255)
                                                    aTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGa                                                 <  1:77149/100‑1 (MQ=255)
                                                     ttCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGAc                                                >  1:1026395/1‑100 (MQ=255)
                                                     ttCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGAc                                                <  1:1895658/100‑1 (MQ=255)
                                                     ttCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGAc                                                <  1:1955260/100‑1 (MQ=255)
                                                       cgccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTg                                              >  2:778145/1‑100 (MQ=255)
                                                       cgccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTg                                              >  1:826654/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)
                                                                tactacCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCtt                                     <  1:3773269/100‑1 (MQ=255)
                                                                 actacCCCACTACGGACGGTCGCGAAGAAGCCATAAACATGAAGTAGCCAGAAAATATTAAATACAAGATGGCAGAATAAAATGGGAATGGAtatcc                                       >  2:1266362/1‑93 (MQ=255)
                                                                  ctacCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTg                                   >  1:1544089/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                               <  1:3401943/100‑1 (MQ=255)
                                                                        caccacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCg                             >  1:3250825/1‑100 (MQ=255)
                                                                          ccacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGAt                           <  2:278978/100‑1 (MQ=255)
                                                                          ccacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGAt                           >  2:3333728/1‑100 (MQ=255)
                                                                             cggacggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATATAAGGTGGAAAAATGACGTGGGACTGGATTTTCTTTGATGCCGATGaa                        >  2:2661394/1‑100 (MQ=255)
                                                                                 cggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGc                    >  1:1093037/1‑100 (MQ=255)
                                                                                  ggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCt                   >  1:2193077/1‑100 (MQ=255)
                                                                                   gTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTg                  <  2:2912185/100‑1 (MQ=255)
                                                                                   gTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTg                  >  1:2292552/1‑100 (MQ=255)
                                                                                    tCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGAGGACTGGATTTTCTTTGATGCCGATGAAACGCTGtt                <  2:3302468/100‑1 (MQ=255)
                                                                                    gattggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGt                 >  1:2136990/9‑100 (MQ=38)
                                                                                     cgcgAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGtt                >  1:2827544/1‑100 (MQ=255)
                                                                                     cgcgAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGt                 >  1:1074650/1‑99 (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)
                                                                                         gcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcc            >  1:556216/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:3510725/3‑100 (MQ=255)
                                                                                          caaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCt           >  2:1170905/3‑100 (MQ=255)
                                                                                           aaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCtt          >  1:2641395/2‑100 (MQ=255)
                                                                                           aaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCtt          >  1:44300/2‑100 (MQ=255)
                                                                                             cGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTg        <  2:3176313/100‑1 (MQ=255)
                                                                                             cGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTg        >  1:1861446/1‑100 (MQ=255)
                                                                                               ccATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGAc      <  2:610950/100‑1 (MQ=255)
                                                                                               ccATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGAc      >  1:1462206/1‑100 (MQ=255)
                                                                                                catcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACt     <  1:1879753/100‑1 (MQ=255)
                                                                                                   catcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCAt  <  2:2062159/100‑1 (MQ=255)
                                                                                                   catcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCAt  <  2:2782839/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCAT  >  NC_000913/4608505‑4608703

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: