Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,608,607 C→G 15.5% 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,6070CG15.5% 326.7 / 43.6 116I141M (ATC→ATGrimIribosomal‑protein‑S18‑alanine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (51/47);  new base G (13/5);  total (64/52)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 128 total aligned reads displayed.

CTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCA  >  NC_000913/4608510‑4608706
                                                                                                 |                                                                                                   
cTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatcat                                                                                                   <  2:109843/100‑1 (MQ=255)
cTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatcat                                                                                                   >  1:81915/1‑100 (MQ=255)
cTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatcat                                                                                                   >  1:2514806/1‑100 (MQ=255)
cTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatcat                                                                                                   <  2:496166/100‑1 (MQ=255)
  gccgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATgg                                                                                                 <  1:424515/100‑1 (MQ=255)
  gccgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATgg                                                                                                 >  2:2111217/1‑100 (MQ=255)
     gccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGt                                                                                              <  2:1234149/100‑1 (MQ=255)
     gccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGt                                                                                              >  1:1921924/1‑100 (MQ=255)
       cATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTg                                                                                            <  1:1394654/100‑1 (MQ=255)
        aTTGCCCTGTACGAAAGTTTAGGCTTTAACGTGGCGACGATTCGCCGCAATTACTACCCCTCCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGc                                                                                           <  2:1337558/100‑1 (MQ=255)
        aTTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCTTTGc                                                                                           <  2:102255/100‑1 (MQ=255)
         ttGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGcc                                                                                          <  2:2069690/100‑1 (MQ=255)
          tGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                         <  2:1569198/100‑6 (MQ=255)
          tGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                         <  2:1859001/100‑6 (MQ=255)
           gCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCaa                                                                                        >  2:1560389/1‑100 (MQ=255)
            cccTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAAt                                                                                       <  1:538844/100‑1 (MQ=255)
            cccTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAAt                                                                                       >  1:411510/1‑100 (MQ=255)
            cccTGTACGAAAGTTTAGGCTTTAACGAGGCGACGAGTCGCCGCAATTACTACCCCCCCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcgcg                                                                                       <  1:2289278/100‑8 (MQ=255)
              cTGTACTAAAGTTTAGGCTTTAACGCGGCGACGATTCGCCGCAATTACTTCCCCCCCACGGACGGTCGCGAAGACGCCATGATCATGGCGTTGCCAATc                                                                                      <  1:2329669/99‑1 (MQ=255)
               tggaCGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGATGCCATCATCATGGCGTTGCCAATCAg                                                                                    <  1:358604/97‑1 (MQ=255)
               tGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAg                                                                                    >  1:2396612/1‑100 (MQ=255)
                   cGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATg                                                                                <  1:1055498/100‑1 (MQ=255)
                   cGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATg                                                                                >  2:2522845/1‑100 (MQ=255)
                   cGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGACATCATCATGGCGTTGCCAATCAGTATg                                                                                >  1:703817/1‑100 (MQ=255)
                     aaaGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTa                                                                              >  1:1796778/1‑100 (MQ=255)
                       aGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAAt                                                                            >  1:1366820/1‑100 (MQ=255)
                         tttAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATAc                                                                          <  1:1676045/100‑1 (MQ=255)
                             ggCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCGCCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAgg                                                                      <  1:1463133/100‑1 (MQ=255)
                             ggCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAgg                                                                      <  2:272540/100‑1 (MQ=255)
                               cTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTg                                                                    >  2:1084671/1‑100 (MQ=255)
                               cTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTg                                                                    >  2:96199/1‑100 (MQ=255)
                                  tAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaa                                                                 <  2:2140777/100‑1 (MQ=255)
                                    aCGCGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaata                                                               <  2:1714632/100‑1 (MQ=255)
                                    aCGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaata                                                               <  2:352024/100‑1 (MQ=255)
                                    aCGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaata                                                               <  1:2322326/100‑1 (MQ=255)
                                      attgcggcgaATCGTCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaataat                                                             >  1:1896360/8‑100 (MQ=255)
                                       aGGCGACGATTCGCCGCAATTACTACCCCACCACGGACAAGCGCGAAGACGCCAACATCATAGCGTTGCCAATCAGTATTTAATACAAGGTGGAATAATg                                                            <  1:529031/100‑1 (MQ=255)
                                        ggCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGa                                                           >  2:133251/1‑100 (MQ=255)
                                        ggCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGa                                                           >  1:1320249/1‑100 (MQ=255)
                                        ggCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGa                                                           <  1:2169697/100‑1 (MQ=255)
                                          cgacgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAg                                                         >  1:1196213/1‑100 (MQ=255)
                                          cgacgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGATGCCAATCAGTATGTAATACAAGGTGGAATAATGAAg                                                         >  1:1741266/1‑100 (MQ=255)
                                           gacgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGt                                                        >  2:1338946/1‑100 (MQ=255)
                                           gacgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGt                                                        <  2:1037413/100‑1 (MQ=255)
                                             cgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTgg                                                      <  2:2021501/100‑1 (MQ=255)
                                              gaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTggg                                                     <  1:1723079/100‑1 (MQ=255)
                                                ttCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGATGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGAc                                                   >  1:1659089/1‑100 (MQ=255)
                                                 tCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACt                                                  <  2:960014/100‑1 (MQ=255)
                                                  cgccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTg                                                 <  2:1683635/100‑1 (MQ=255)
                                                  cgccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTg                                                 <  1:1956362/100‑1 (MQ=255)
                                                   gccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTgg                                                >  2:2563510/1‑100 (MQ=255)
                                                            actacCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCttt                                       <  2:1015709/100‑1 (MQ=255)
                                                             ctacCCCACCAGGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTg                                      <  1:1501385/100‑1 (MQ=255)
                                                             ctacCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTg                                      <  2:1660686/100‑1 (MQ=255)
                                                              tacCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGa                                     >  1:2216413/1‑100 (MQ=255)
                                                               acCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGAt                                    >  1:2567775/1‑100 (MQ=255)
                                                                ccccACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATg                                   >  2:2330582/1‑100 (MQ=255)
                                                                ccccACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATg                                   >  2:2113174/1‑100 (MQ=255)
                                                                 cccACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGc                                  >  1:801968/1‑100 (MQ=255)
                                                                  ccaccaCGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGcc                                 >  2:1309057/1‑100 (MQ=255)
                                                                  ccaccaCGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGcc                                 >  2:1135434/1‑100 (MQ=255)
                                                                  ccaccaCGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGcc                                 >  1:421432/1‑100 (MQ=255)
                                                                  ccaccaCGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGcc                                 <  2:1659098/100‑1 (MQ=255)
                                                                   caccacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCg                                >  1:469880/1‑100 (MQ=255)
                                                                   caccacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCg                                >  1:1003851/1‑100 (MQ=255)
                                                                    accacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGa                               <  2:632186/100‑1 (MQ=255)
                                                                      cacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATg                             <  1:759639/100‑1 (MQ=255)
                                                                      cacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATg                             >  1:289198/1‑100 (MQ=255)
                                                                      cacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATg                             >  2:2000018/1‑100 (MQ=255)
                                                                              gTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTg                     <  1:1412040/100‑1 (MQ=255)
                                                                              gTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTg                     <  1:2307073/100‑1 (MQ=255)
                                                                              gTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGAGGAAATGGATTTTCTTTGATGCCGATGAAACGCTg                     >  2:1204107/1‑100 (MQ=255)
                                                                               tCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGt                    <  1:1897699/100‑1 (MQ=255)
                                                                               tCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGt                    <  1:2381775/100‑1 (MQ=255)
                                                                               gattggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCt                      >  1:953632/9‑98 (MQ=38)
                                                                               gattggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTg                     >  1:571309/9‑99 (MQ=38)
                                                                                  tggcaaCGCCATGATGATGGCGTTGCCAATGAGAATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGATACGCTGTTTa                 >  1:1697246/6‑100 (MQ=255)
                                                                                   ggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAc                >  1:916145/5‑100 (MQ=255)
                                                                                   ggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAc                >  2:2189141/5‑100 (MQ=255)
                                                                                   gAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAc                >  2:2475023/1‑100 (MQ=255)
                                                                                    gcaaCGCCGTGATGATGGCGGTGCCAATCAGTATGTAAAAGAAGGGGGAATAATGAAGTGGGACGGGCTTTTCTTTGATGCCGGTGAACCGGTGttga                 >  1:2145626/4‑96 (MQ=255)
                                                                                    gcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAGTACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcc               >  1:1108713/4‑100 (MQ=255)
                                                                                    gcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcc               >  2:1480015/4‑100 (MQ=255)
                                                                                    gcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcc               >  2:1236239/4‑100 (MQ=255)
                                                                                    aaGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcc               <  2:2432885/100‑1 (MQ=255)
                                                                                     caaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCt              >  1:1519781/3‑100 (MQ=255)
                                                                                     caaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCt              >  1:1617772/3‑100 (MQ=255)
                                                                                       aCGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCttt            >  2:1972820/1‑100 (MQ=255)
                                                                                        cGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTg           >  2:979157/1‑100 (MQ=255)
                                                                                        cGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTg           >  2:643306/1‑100 (MQ=255)
                                                                                        cGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTg           <  1:643303/100‑1 (MQ=255)
                                                                                        cGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACATTTg           >  2:2043018/1‑100 (MQ=255)
                                                                                          cggtcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGAc         <  2:229265/97‑1 (MQ=255)
                                                                                          ccATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGAc         <  2:2431026/100‑1 (MQ=255)
                                                                                           catcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACt        <  2:779396/100‑1 (MQ=255)
                                                                                           catcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACt        <  2:1016236/100‑1 (MQ=255)
                                                                                            atcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCAATGAAACGCTGTTTACCTTTg           <  1:309960/96‑1 (MQ=255)
                                                                                             tcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCa      >  2:1095839/1‑100 (MQ=255)
                                                                                              catcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCAt     >  2:2160506/1‑100 (MQ=255)
                                                                                                 catGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCa  >  2:1837506/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
CTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCA  >  NC_000913/4608510‑4608706

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