Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,608,604 C→G 5.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,6040CG 5.1% 365.5 / 9.8 118I140M (ATC→ATGrimIribosomal‑protein‑S18‑alanine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (60/51);  new base G (2/4);  total (63/55)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

Only 100 of 124 total aligned reads displayed.

AACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTACCGCT  >  NC_000913/4608456‑4608753
                                                                                                                                                    |                                                                                                                                                     
aaCTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatca                                                                                                                                                      <  2:235254/150‑1 (MQ=255)
  cTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatc                                                                                                                                                    >  1:525967/1‑150 (MQ=255)
    ggAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatcat                                                                                                                                                  <  1:1077764/150‑1 (MQ=255)
    ggAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCcatcatcat                                                                                                                                                  >  2:164938/1‑150 (MQ=255)
      aaaaaCGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATgg                                                                                                                                                >  2:631308/1‑150 (MQ=255)
       aaaaCGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGc                                                                                                                                               >  1:646262/1‑150 (MQ=255)
           cgcgGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtct                                                                                                                                           <  1:416893/150‑3 (MQ=255)
            gcggcgTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACTAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGc                                                                                                                                          <  1:1429610/150‑1 (MQ=255)
            gcggcgTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGc                                                                                                                                          <  1:1651776/150‑1 (MQ=255)
             cggcgTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGcc                                                                                                                                         >  2:1412871/1‑150 (MQ=255)
              ggcgTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcg                                                                                                                                        <  2:1158652/150‑6 (MQ=255)
                cgTGGCGACACTCTGGCTGGACGGCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGGTTAGGCTTTAACGAGGCGACGATGCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGtcttcgcg                                                                                                                                      <  1:84029/150‑8 (MQ=255)
                 gTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCCCCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATc                                                                                                                                     >  1:1084645/1‑150 (MQ=255)
                     cGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTa                                                                                                                                 >  1:68645/1‑150 (MQ=255)
                        cacTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGt                                                                                                                              >  1:1647231/1‑150 (MQ=255)
                         acTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTa                                                                                                                             <  2:650305/150‑1 (MQ=255)
                             tGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGTCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATAc                                                                                                                         >  1:70607/1‑150 (MQ=255)
                                 tGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAgg                                                                                                                     <  1:997510/150‑1 (MQ=255)
                                 tGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAgg                                                                                                                     >  1:1483952/1‑150 (MQ=255)
                                 tGGAAGTCCGTGCTTCAAACGCTGCCGACATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAgg                                                                                                                     >  1:1029316/1‑150 (MQ=255)
                                  ggAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGt                                                                                                                    <  1:1542946/150‑1 (MQ=255)
                                    aaGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTgg                                                                                                                  <  1:42333/150‑1 (MQ=255)
                                       tCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGTaat                                                                                                               <  2:1137727/150‑1 (MQ=255)
                                       tCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaat                                                                                                               >  1:758724/1‑150 (MQ=255)
                                           tGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATg                                                                                                           <  2:95229/150‑1 (MQ=255)
                                            gCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGa                                                                                                          <  2:1771458/150‑1 (MQ=255)
                                             cTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGCGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGCCGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAAGAATGaa                                                                                                         >  2:1012358/1‑150 (MQ=255)
                                             cTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGCGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGaa                                                                                                         >  2:150298/1‑150 (MQ=255)
                                              ttCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACTATTCGACGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAg                                                                                                        <  1:75362/150‑1 (MQ=255)
                                              ttCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAg                                                                                                        <  1:131423/150‑1 (MQ=255)
                                               tCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGt                                                                                                       <  1:415757/150‑1 (MQ=255)
                                               tCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGt                                                                                                       >  1:755660/1‑150 (MQ=255)
                                                cAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTg                                                                                                      <  1:156940/150‑1 (MQ=255)
                                                  aaCGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTggg                                                                                                    >  1:1208143/1‑150 (MQ=255)
                                                     gCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGCCGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGAcg                                                                                                 >  1:763679/1‑149 (MQ=255)
                                                      cTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCACAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTg                                                                                                >  2:441238/1‑150 (MQ=255)
                                                        gccgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGa                                                                                              >  1:1372979/1‑150 (MQ=255)
                                                                    cTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATg                                                                                  >  1:777204/1‑150 (MQ=255)
                                                                    cTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATg                                                                                  <  2:654334/150‑1 (MQ=255)
                                                                    cTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATg                                                                                  >  1:1217370/1‑150 (MQ=255)
                                                                    cTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATg                                                                                  <  1:1693661/150‑1 (MQ=255)
                                                                    cTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATg                                                                                  <  1:1686712/150‑1 (MQ=255)
                                                                       tACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCg                                                                               <  1:1051794/150‑1 (MQ=255)
                                                                        aCGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACTGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGa                                                                              <  2:579055/150‑1 (MQ=255)
                                                                         cGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTTGATTTTCTTTGATTCCGAt                                                                             >  2:1490868/1‑150 (MQ=255)
                                                                         cGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGAt                                                                             >  1:1270625/1‑150 (MQ=255)
                                                                          gAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACAACGGACGGTCGCGATGACGCCATCATCATGGCGTTGACAATCAGTATGTAATACAAGGTGGAATAATGAAGTCGGACTGGATTTTCTTTTCTGCCGATg                                                                            >  1:925081/1‑150 (MQ=255)
                                                                           aaaGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGGATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGa                                                                           >  1:950553/1‑150 (MQ=255)
                                                                              gTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTTATGCCGATGAAAc                                                                        >  2:1640058/1‑150 (MQ=255)
                                                                                ttAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGGCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGc                                                                      <  1:508304/150‑1 (MQ=255)
                                                                                 tAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGTAGCCGCCATCATCATGGCGTTGCCAATCAGGATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGATACGCt                                                                     >  2:59569/1‑150 (MQ=255)
                                                                                  aGGCTTTAACGAGGCGACTATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTg                                                                    <  2:1396809/150‑1 (MQ=255)
                                                                                  aGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTg                                                                    >  1:1697309/1‑150 (MQ=255)
                                                                                  aGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTg                                                                    >  2:1744140/1‑150 (MQ=255)
                                                                                    gCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGGTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGtt                                                                  >  1:899368/1‑150 (MQ=255)
                                                                                     cTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGttt                                                                 >  2:1575250/1‑150 (MQ=255)
                                                                                      tttAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTa                                                                <  2:182311/150‑1 (MQ=255)
                                                                                         aaCGAGGCGACGATTCGCCGCAATTACAACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACAGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcgt                                                             >  1:217004/1‑148 (MQ=255)
                                                                                          aCGAGGCGACGATTCGCCGCAATTACTACCCCCCCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCtt                                                            >  1:1319632/1‑150 (MQ=255)
                                                                                          aCGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCtt                                                            >  1:454656/1‑150 (MQ=255)
                                                                                           cGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCttt                                                           <  2:821675/150‑1 (MQ=255)
                                                                                            gAGGCGACGATTCGCCGCTAGTTCTAGCCCACCACGGACGGTCGCGAAGACGCCATCATGATGGCGTTGCCAATCAGTATGTAATAGAAGGTGGAATAATGTAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTg                                                          <  1:270850/150‑1 (MQ=255)
                                                                                            gAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTg                                                          <  1:1193468/150‑1 (MQ=255)
                                                                                                     aTTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCa                                                 >  1:679659/1‑150 (MQ=255)
                                                                                                     aTTCGCCACAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCa                                                 <  1:441238/150‑1 (MQ=255)
                                                                                                         gccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCgg                                             >  2:740099/1‑150 (MQ=255)
                                                                                                         gccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGAAGAAACGCTGTTTACCTTTGACTCATTCACCgg                                             >  1:703983/1‑150 (MQ=255)
                                                                                                          ccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTCCCTTTGACTCATTCACCGGc                                            >  2:732877/1‑150 (MQ=255)
                                                                                                          ccgcAATTACTACCCCACCACGGACGGTAGCGAAGACGCCATCATCATGGCGGTGCCAATCAGTTTGTTATCCAAGGTGGAATAATGATGAGGGACTTGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACACATTCACCGGc                                            >  1:1618636/1‑150 (MQ=255)
                                                                                                          ccgcAATTACTACCCCACCACGGACGGGCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGGGGAATACTGAAGTGGGACTGGATGTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGc                                            <  1:1226500/150‑1 (MQ=255)
                                                                                                           cgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGcc                                           >  2:1412274/1‑150 (MQ=255)
                                                                                                           cgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGcc                                           >  2:962122/1‑150 (MQ=255)
                                                                                                                 tactacCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGc                                     <  2:1461828/150‑1 (MQ=255)
                                                                                                                    tacCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACATGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGa                                  >  1:647248/1‑150 (MQ=255)
                                                                                                                    tacCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGa                                  <  2:68645/150‑1 (MQ=255)
                                                                                                                     acCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGAt                                 <  2:1647231/150‑1 (MQ=255)
                                                                                                                     acCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGAt                                 >  1:1722355/1‑150 (MQ=255)
                                                                                                                      ccccACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATg                                <  2:1208143/150‑1 (MQ=255)
                                                                                                                        ccaccaCGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGtt                              <  2:70607/150‑1 (MQ=255)
                                                                                                                        ccaccaCGGACGGTCCCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGATGTTTACCTTTGACTCATACACCGGCCTGCAGCGGATGtt                              >  2:434106/1‑150 (MQ=255)
                                                                                                                         caccacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGttt                             >  2:1247001/1‑150 (MQ=255)
                                                                                                                           ccacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCt                           >  2:250910/1‑150 (MQ=255)
                                                                                                                            cacGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCtt                          <  2:1084645/150‑1 (MQ=255)
                                                                                                                             acggacggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTg                         <  1:1435468/150‑1 (MQ=255)
                                                                                                                               ggacggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGAt                       <  1:1408541/150‑1 (MQ=255)
                                                                                                                                gacggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGAtt                      <  2:355602/150‑1 (MQ=255)
                                                                                                                                gacggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGAtt                      <  2:1245759/150‑1 (MQ=255)
                                                                                                                                   ggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACa                   <  1:1389916/150‑1 (MQ=255)
                                                                                                                                    gTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAg                  >  2:851334/1‑150 (MQ=255)
                                                                                                                                    gTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCAGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAg                  <  1:478343/150‑1 (MQ=255)
                                                                                                                                      cgcgAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCCCCGGCCTGCAGCGGATGTTTCTTGATTACAGCg                <  1:1066161/150‑1 (MQ=255)
                                                                                                                                      cgcgAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCg                >  2:340909/1‑150 (MQ=255)
                                                                                                                                      cgcgAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCg                >  1:1487695/1‑150 (MQ=255)
                                                                                                                                        tggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACTGCGTc              >  1:49968/6‑150 (MQ=255)
                                                                                                                                         taaGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCa             <  1:557832/149‑1 (MQ=255)
                                                                                                                                         ggcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCa             >  2:492644/5‑150 (MQ=255)
                                                                                                                                                ccATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTAc      <  2:758724/150‑1 (MQ=255)
                                                                                                                                                ccATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGAAGTTTCTTGATTACAGCGTCACCTTTAc      <  2:1483952/150‑1 (MQ=255)
                                                                                                                                                  atcatcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTACCg    <  2:975086/150‑1 (MQ=255)
                                                                                                                                                    catcatGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTACCGCt  <  2:561412/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCACCTTTACCGCT  >  NC_000913/4608456‑4608753

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