Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,643,987 A→G 6.9% V99V (GTT→GTC) ‡ rlmN ← dual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,643,9870AG 6.9% 405.3 / 16.5 146V99V (GTT→GTC) ‡rlmNdual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (72/62);  new base G (3/7);  total (76/70)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.98e-01

Only 100 of 147 total aligned reads displayed.

TTTCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATA  >  NC_000913/2643842‑2644136
                                                                                                                                                 |                                                                                                                                                     
tttCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACgg                                                                                                                                                   <  2:624411/150‑1 (MQ=255)
tttCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACgg                                                                                                                                                   <  2:587129/150‑1 (MQ=255)
tttCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACgg                                                                                                                                                   <  2:1017157/150‑1 (MQ=255)
  tCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCg                                                                                                                                                 <  2:1122131/150‑1 (MQ=255)
   ccGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGa                                                                                                                                                >  1:843517/1‑150 (MQ=255)
    cGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGCTCGCCAACGGCgct                                                                                                                                               >  1:300529/1‑148 (MQ=255)
     gACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATc                                                                                                                                              >  1:75260/1‑150 (MQ=255)
     gACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATc                                                                                                                                              <  1:757876/150‑1 (MQ=255)
         cccGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCGTTGGCGATcagcg                                                                                                                                          <  2:1180364/150‑5 (MQ=255)
         cccGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGccc                                                                                                                                          >  2:218956/1‑150 (MQ=255)
          ccGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCTAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCa                                                                                                                                         >  2:1459415/1‑150 (MQ=255)
          ccGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCa                                                                                                                                         <  1:1711370/150‑1 (MQ=255)
           cGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATAAACCGTTTCGACGCGCTGATCGCCGTTGGCGATcagcgcg                                                                                                                                        <  2:1339767/150‑7 (MQ=255)
            gCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCGTTGGCGATcagcgcgt                                                                                                                                       <  1:38792/150‑8 (MQ=255)
            gCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCAtt                                                                                                                                       <  1:218956/150‑1 (MQ=255)
               ggTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGGCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTaa                                                                                                                                    <  1:1627163/150‑1 (MQ=255)
                    ccgtaaaaGCCCTGCAGGGCGGTGGAACACAATTTCCACTCCAGCGCACACCCCACCTGCGACGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTg                                                                                                                               <  2:1610679/145‑1 (MQ=255)
                    cGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGCTCGCCCATTTAATGGTg                                                                                                                               >  2:1461986/1‑150 (MQ=255)
                    cGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTg                                                                                                                               >  2:1027301/1‑150 (MQ=255)
                     ggTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTAAACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGc                                                                                                                              <  2:1289687/150‑1 (MQ=255)
                          aaGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCa                                                                                                                         >  2:1281095/1‑150 (MQ=255)
                            gCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGa                                                                                                                       >  2:1176631/1‑150 (MQ=255)
                            gCCCTGCGGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTCTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGa                                                                                                                       <  1:294289/150‑1 (MQ=255)
                             cccTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGAt                                                                                                                      >  1:359315/1‑150 (MQ=255)
                             cccTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGGCAGAt                                                                                                                      >  2:630461/1‑150 (MQ=255)
                                  ctgGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACg                                                                                                                 >  1:1344492/1‑150 (MQ=255)
                                   ggggCGGTGGAACTGATTTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGc                                                                                                                <  2:1310268/149‑1 (MQ=255)
                                      gCGGTGGAACAGAATTTACACTCCAGCGCACACACCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGt                                                                                                             >  2:891482/1‑150 (MQ=255)
                                        ggTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGttc                                                                                                           >  2:15214/1‑150 (MQ=255)
                                          tGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGttctt                                                                                                         <  1:132173/150‑1 (MQ=255)
                                           ggAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGttcttc                                                                                                        >  1:526545/1‑150 (MQ=255)
                                           ggAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGttcttc                                                                                                        >  1:939573/1‑150 (MQ=255)
                                              aCAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACTCAGAGCGTTGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCACCGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAAc                                                                                                     >  2:1191597/1‑150 (MQ=255)
                                              aCAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAAc                                                                                                     >  2:1191601/1‑150 (MQ=255)
                                               cAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAAcc                                                                                                    <  1:127799/150‑1 (MQ=255)
                                                aGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCa                                                                                                   >  1:226638/1‑150 (MQ=255)
                                                       acacTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCgg                                                                                            <  1:1227161/150‑1 (MQ=255)
                                                         acTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTg                                                                                          <  2:188386/150‑1 (MQ=255)
                                                          cTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTgc                                                                                         >  1:1682014/1‑150 (MQ=255)
                                                          cTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTgc                                                                                         >  1:1400772/1‑150 (MQ=255)
                                                           tCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTgcg                                                                                        >  2:554309/1‑150 (MQ=255)
                                                           tCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTgcg                                                                                        <  2:75260/150‑1 (MQ=255)
                                                                cgcACTCCCCTCCTGCGTAGAGACGCAGAGCGTGGCACGGGCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGAt                                                                                   <  1:1087970/150‑1 (MQ=255)
                                                                   acacCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTc                                                                                >  2:913850/1‑150 (MQ=255)
                                                                    cacCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTcc                                                                               <  2:1400772/150‑1 (MQ=255)
                                                                     ccgcaaCCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCg                                                                              <  1:1306222/145‑1 (MQ=255)
                                                                       cccACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGcc                                                                            >  1:1522021/1‑150 (MQ=255)
                                                                          aCCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCAcc                                                                         >  1:979910/1‑150 (MQ=255)
                                                                          aCCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGAGCGGATTTCCGCCAcc                                                                         >  2:451055/1‑150 (MQ=255)
                                                                             tGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACctct                                                                      <  1:1176631/150‑1 (MQ=255)
                                                                               cGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCttt                                                                    >  2:178603/1‑150 (MQ=255)
                                                                                gAAGAGACGCAGAGCGTGGCACGGGCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTc                                                                   <  1:630461/150‑1 (MQ=255)
                                                                                      aCGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTg                                                             <  1:178603/150‑1 (MQ=255)
                                                                                       cGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGc                                                            <  2:846984/150‑1 (MQ=255)
                                                                                       cGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGc                                                            >  2:150290/1‑150 (MQ=255)
                                                                                       cGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGc                                                            >  1:487448/1‑150 (MQ=255)
                                                                                       cGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGc                                                            >  2:741456/1‑150 (MQ=255)
                                                                                       cGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGc                                                            >  2:28856/1‑150 (MQ=255)
                                                                                         cAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCcg                                                          >  1:500929/1‑150 (MQ=255)
                                                                                         cAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCcg                                                          <  1:1281095/150‑1 (MQ=255)
                                                                                          agagCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCcgc                                                         <  1:1004185/150‑1 (MQ=255)
                                                                                           gagCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCcgcg                                                        <  1:15214/150‑1 (MQ=255)
                                                                                            agCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCgcgc                                                       <  1:1662073/150‑1 (MQ=255)
                                                                                            agCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCgcgc                                                       <  2:293907/150‑1 (MQ=255)
                                                                                              cGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCaa                                                     >  1:1074209/1‑150 (MQ=255)
                                                                                                   cACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACtt                                                >  1:947321/1‑150 (MQ=255)
                                                                                                    aCGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACttt                                               <  2:843517/150‑1 (MQ=255)
                                                                                                    aCGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACttt                                               <  2:1344492/150‑1 (MQ=255)
                                                                                                      ggTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCTCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgt                                             <  1:1339579/150‑1 (MQ=255)
                                                                                                      ggTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgt                                             <  2:526545/150‑1 (MQ=255)
                                                                                                       gtcgtcTTCCGGTCTATTCACCGTTTCGACGCGCTGATCGCCAAAGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgtt                                            <  1:146607/150‑1 (MQ=255)
                                                                                                       gtcgtcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTCATGGTGCCGTCAGATGAACGCTGTGCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgtt                                            <  1:1115027/150‑1 (MQ=255)
                                                                                                       gtcgtcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgtt                                            >  2:1091084/1‑150 (MQ=255)
                                                                                                         cgtcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGa                                          >  1:895332/1‑150 (MQ=255)
                                                                                                           tcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATg                                        <  1:540576/150‑1 (MQ=255)
                                                                                                             ttCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATATc                                      <  1:1191601/150‑1 (MQ=255)
                                                                                                               ccGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCgg                                    <  2:142799/150‑1 (MQ=255)
                                                                                                                  ggATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGtca                                 <  2:1337686/150‑1 (MQ=255)
                                                                                                                       aacacCGTTTCGACGCGCTGATCGCCAACGTCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGtcatctca                            <  1:1276833/149‑1 (MQ=255)
                                                                                                                         cacCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGtcatctcatc                          >  1:1176656/1‑150 (MQ=255)
                                                                                                                          acCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCtcatca                         <  2:1355370/150‑1 (MQ=255)
                                                                                                                                 cGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGt                  >  1:1120929/1‑150 (MQ=255)
                                                                                                                                    cgcgCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgc               <  1:28856/150‑1 (MQ=255)
                                                                                                                                      gggcgATCGCCGTTGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCaa                        >  1:2131/5‑139 (MQ=255)
                                                                                                                                      gggcgATCGCCGTTGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCaa                        <  2:2131/135‑1 (MQ=255)
                                                                                                                                      cgcTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcag             <  2:226638/150‑1 (MQ=255)
                                                                                                                                       ggcgATCGCCGTTGGCGATCGCCCATTTTATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCCCCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcagc            >  2:502651/4‑150 (MQ=255)
                                                                                                                                       ggcgATCGCCGTTGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTATTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTATTTCAGTTTGCCGCGCAACCCTTTGTTGATGTCGGTCATCGCATCAAAGTTGTCgcagc            >  2:1486730/4‑150 (MQ=255)
                                                                                                                                       gcTGATCGCCCACGTCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcagc            <  1:260202/150‑1 (MQ=255)
                                                                                                                                       gcTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcagc            <  1:1285655/150‑1 (MQ=255)
                                                                                                                                       gcTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcagc            <  1:258072/150‑1 (MQ=255)
                                                                                                                                        cTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcagca           >  2:1549650/1‑150 (MQ=255)
                                                                                                                                        cTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcagca           >  1:986516/1‑150 (MQ=255)
                                                                                                                                         tGATCGCCAACGGCGATCGCCCATTTAATGGAGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTATCAGCaa          >  2:437022/1‑150 (MQ=255)
                                                                                                                                          gATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAAt         >  1:134799/1‑150 (MQ=255)
                                                                                                                                          gATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGAGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAAt         >  1:13659/1‑150 (MQ=255)
                                                                                                                                               ccAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGa    <  1:554309/150‑1 (MQ=255)
                                                                                                                                               ccAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTACGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGa    >  1:150108/1‑150 (MQ=255)
                                                                                                                                                cAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTACGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGAt   <  2:300529/150‑1 (MQ=255)
                                                                                                                                                 aaCGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATa  <  1:891482/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TTTCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATA  >  NC_000913/2643842‑2644136

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