Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,608,604 C→G 16.0% 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,6040CG16.0% 310.1 / 43.4 116I140M (ATC→ATGrimIribosomal‑protein‑S18‑alanine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (51/46);  new base G (13/6);  total (64/52)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 130 total aligned reads displayed.

AAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCAT  >  NC_000913/4608505‑4608703
                                                                                                   |                                                                                                   
 aaCGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatca                                                                                                    >  1:1890986/1‑100 (MQ=255)
 aaCGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatca                                                                                                    >  1:801361/1‑100 (MQ=255)
     cTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatcat                                                                                                >  1:2514806/1‑100 (MQ=255)
     cTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatcat                                                                                                >  1:81915/1‑100 (MQ=255)
       gccgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATgg                                                                                              <  1:424515/100‑1 (MQ=255)
          gccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGt                                                                                           >  1:1921924/1‑100 (MQ=255)
            cATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTg                                                                                         <  1:1394654/100‑1 (MQ=255)
             aTTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtctt                                                                                        <  2:157792/100‑4 (MQ=255)
             aTTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCTTTGc                                                                                        <  2:102255/100‑1 (MQ=255)
              ttGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGcc                                                                                       <  2:2069690/100‑1 (MQ=255)
              ttGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGcc                                                                                       <  2:1051264/100‑1 (MQ=255)
               tGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                      <  2:1859001/100‑6 (MQ=255)
               tGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                      <  2:1569198/100‑6 (MQ=255)
                gCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGAGGATTCGCGGCAATTACTACCCGACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtct                                                                                         <  1:160680/96‑3 (MQ=255)
                 cccTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAAt                                                                                    <  1:538844/100‑1 (MQ=255)
                 cccTGTACGAAAGTTTAGGCTTTAACGAGGCGACGAGTCGCCGCAATTACTACCCCCCCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcgcg                                                                                    <  1:2289278/100‑8 (MQ=255)
                  ccTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATc                                                                                   >  2:2097979/1‑100 (MQ=255)
                   cTGTACTAAAGTTTAGGCTTTAACGCGGCGACGATTCGCCGCAATTACTTCCCCCCCACGGACGGTCGCGAAGACGCCATGATCATGGCGTTGCCAATc                                                                                   <  1:2329669/99‑1 (MQ=255)
                   cTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATc                                                                                   <  1:397264/99‑1 (MQ=255)
                    tggaCGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGATGCCATCATCATGGCGTTGCCAATCAg                                                                                 <  1:358604/97‑1 (MQ=255)
                    tGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAg                                                                                 >  1:2396612/1‑100 (MQ=255)
                      tACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTa                                                                               <  2:352989/100‑1 (MQ=255)
                        cGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATg                                                                             >  1:805696/1‑100 (MQ=255)
                        cGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATg                                                                             >  2:2522845/1‑100 (MQ=255)
                        cGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATg                                                                             <  1:1055498/100‑1 (MQ=255)
                        cGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGACATCATCATGGCGTTGCCAATCAGTATg                                                                             >  1:703817/1‑100 (MQ=255)
                          aaaGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTa                                                                           >  1:1799295/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)
                                  ggCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAgg                                                                   >  2:1977149/1‑100 (MQ=255)
                                    cTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTg                                                                 >  2:96199/1‑100 (MQ=255)
                                    cTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTg                                                                 >  2:1084671/1‑100 (MQ=255)
                                       tAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaa                                                              <  2:2140777/100‑1 (MQ=255)
                                         aCGCGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaata                                                            <  2:1714632/100‑1 (MQ=255)
                                         aCGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaata                                                            <  1:2322326/100‑1 (MQ=255)
                                         aCGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaata                                                            <  2:352024/100‑1 (MQ=255)
                                          cGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaataa                                                           >  1:2376834/1‑100 (MQ=255)
                                           gAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaataat                                                          >  2:1765520/1‑100 (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)
                                               cgacgaTTCGCCGCAATTACTACCCCACCACGGACTGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTTGAATAATGAAg                                                      <  1:266566/100‑1 (MQ=255)
                                               cgacgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAg                                                      >  1:1196213/1‑100 (MQ=255)
                                               cgacgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAg                                                      >  1:633849/1‑100 (MQ=255)
                                               cgacgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGATGCCAATCAGTATGTAATACAAGGTGGAATAATGAAg                                                      >  1:1741266/1‑100 (MQ=255)
                                                gacgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGt                                                     >  2:1338946/1‑100 (MQ=255)
                                                  cgaTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTgg                                                   <  2:1306236/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                                              <  1:1956362/100‑1 (MQ=255)
                                                       cgccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTg                                              <  2:1683635/100‑1 (MQ=255)
                                                        gccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTgg                                             >  2:2563510/1‑100 (MQ=255)
                                                        gccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTgg                                             >  2:1459117/1‑100 (MQ=255)
                                                          cgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGAt                                           <  2:1154174/100‑1 (MQ=255)
                                                                 actacCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCttt                                    <  2:1015709/100‑1 (MQ=255)
                                                                  ctacCCCACCAGGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTg                                   <  1:1501385/100‑1 (MQ=255)
                                                                   tacCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGa                                  >  1:2216413/1‑100 (MQ=255)
                                                                   tacCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGa                                  >  2:217485/1‑100 (MQ=255)
                                                                    acCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGAt                                 >  1:2567775/1‑100 (MQ=255)
                                                                     ccccACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATg                                >  2:2113174/1‑100 (MQ=255)
                                                                     ccccACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATg                                >  2:2330582/1‑100 (MQ=255)
                                                                       ccaccaCGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGcc                              <  2:1659098/100‑1 (MQ=255)
                                                                       ccaccaCGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGcc                              >  2:1135434/1‑100 (MQ=255)
                                                                       ccaccaCGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGcc                              >  2:1309057/1‑100 (MQ=255)
                                                                       ccaccaCGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGcc                              >  1:421432/1‑100 (MQ=255)
                                                                        caccacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCg                             >  1:469880/1‑100 (MQ=255)
                                                                        caccacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCg                             >  1:1003851/1‑100 (MQ=255)
                                                                         accacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGa                            <  2:632186/100‑1 (MQ=255)
                                                                           cacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATg                          >  2:2000018/1‑100 (MQ=255)
                                                                           cacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATg                          <  1:759639/100‑1 (MQ=255)
                                                                                   gTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGAGGAAATGGATTTTCTTTGATGCCGATGAAACGCTg                  >  2:1204107/1‑100 (MQ=255)
                                                                                   gTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAATATGTAATACAAGGTGGAAGAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGc                    <  1:1414039/98‑1 (MQ=255)
                                                                                   gTCGCGAAGACGCCATCATCATGGCGTTGCCAAGCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTg                  <  1:2235085/100‑1 (MQ=255)
                                                                                    gattggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTg                  >  1:571309/9‑99 (MQ=38)
                                                                                        ggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAc             >  2:2368755/5‑100 (MQ=255)
                                                                                        ggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAc             >  1:916145/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:1617772/3‑100 (MQ=255)
                                                                                            aCGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCttt         >  2:1972820/1‑100 (MQ=255)
                                                                                             cGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTg        >  2:643306/1‑100 (MQ=255)
                                                                                             cGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTg        >  2:979157/1‑100 (MQ=255)
                                                                                             cGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTg        >  1:2011027/1‑100 (MQ=255)
                                                                                             cGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACATTTg        >  2:2043018/1‑100 (MQ=255)
                                                                                               cggtcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGAc      <  2:229265/97‑1 (MQ=255)
                                                                                               ccATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGAc      <  1:979154/100‑1 (MQ=255)
                                                                                                catcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACt     <  2:1016236/100‑1 (MQ=255)
                                                                                                catcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACt     <  2:779396/100‑1 (MQ=255)
                                                                                                 atcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTc    >  1:1709194/1‑100 (MQ=255)
                                                                                                 atcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCAATGAAACGCTGTTTACCTTTg        <  1:309960/96‑1 (MQ=255)
                                                                                                  tcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCa   >  2:1095839/1‑100 (MQ=255)
                                                                                                   catcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCAt  >  2:2160506/1‑100 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: