Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,643,988 A→T 6.1% V99D (GTT→GAT) ‡ 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,9880AT 6.1% 369.1 / 8.9 146V99D (GTT→GAT) ‡rlmNdual specificity 23S rRNA m(2)A2503, tRNA m(2)A37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (70/62);  new base T (3/7);  total (76/70)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

Only 100 of 148 total aligned reads displayed.

TTTCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATAC  >  NC_000913/2643842‑2644137
                                                                                                                                                  |                                                                                                                                                     
tttCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACgg                                                                                                                                                    <  2:624411/150‑1 (MQ=255)
tttCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACgg                                                                                                                                                    <  2:587129/150‑1 (MQ=255)
  tCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCg                                                                                                                                                  <  2:1122131/150‑1 (MQ=255)
   ccGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGa                                                                                                                                                 >  2:1004185/1‑150 (MQ=255)
     gACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATc                                                                                                                                               <  1:757876/150‑1 (MQ=255)
         cccGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCGTTGGCGATcagcg                                                                                                                                           <  2:1180364/150‑5 (MQ=255)
          ccGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCTAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCa                                                                                                                                          >  2:1459415/1‑150 (MQ=255)
          ccGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCa                                                                                                                                          <  1:1711370/150‑1 (MQ=255)
          ccGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGGCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCGTTGGCGATcagcgc                                                                                                                                          <  1:547642/150‑6 (MQ=255)
            gCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCAtt                                                                                                                                        <  1:218956/150‑1 (MQ=255)
            gCAGGTTGCGGTTAAAGCCCTGCTGGGCGCTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTATCGACGCGCTGATCGCCAACGGCGATCGCCCAtt                                                                                                                                        >  2:1609340/1‑150 (MQ=255)
               ggTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGGCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTaa                                                                                                                                     <  1:1627163/150‑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                                                                                                                          <  1:1671514/150‑1 (MQ=255)
                            gCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGa                                                                                                                        >  1:293907/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)
                              ccTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATg                                                                                                                      >  2:1339579/1‑150 (MQ=255)
                                  ctgGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACg                                                                                                                  >  1:1344492/1‑150 (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)
                                              aCAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAAc                                                                                                      >  2:1191601/1‑150 (MQ=255)
                                               cAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAAcc                                                                                                     <  1:127799/150‑1 (MQ=255)
                                                aGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCa                                                                                                    >  1:226638/1‑150 (MQ=255)
                                                 gAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCAc                                                                                                   >  2:21692/1‑150 (MQ=255)
                                                  aaTTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGGTTCGACGCGCTGATCGCACCCGGCGATCGCCCATTTAATGGGGCCGTCAGCTGAACGCTGTTCGTCAACCACt                                                                                                  >  1:1518769/1‑150 (MQ=255)
                                                       acacTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCgg                                                                                             <  1:1227161/150‑1 (MQ=255)
                                                        cacTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGt                                                                                            >  1:1355869/1‑150 (MQ=255)
                                                          cTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACCCCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACATCCGGTgc                                                                                          >  1:96304/1‑150 (MQ=255)
                                                          cTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTgc                                                                                          >  1:1682014/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)
                                                                    cacCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTcc                                                                                >  1:735644/1‑150 (MQ=255)
                                                                     ccgcaaCCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCg                                                                               <  1:1306222/145‑1 (MQ=255)
                                                                     acCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGTACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCg                                                                               <  1:894715/150‑1 (MQ=255)
                                                                          aCCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGAGCGGATTTCCGCCAcc                                                                          >  2:451055/1‑150 (MQ=255)
                                                                           ccTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACct                                                                         >  2:183985/1‑150 (MQ=255)
                                                                             tGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACctct                                                                       <  1:1176631/150‑1 (MQ=255)
                                                                               cGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCttt                                                                     >  2:1182156/1‑150 (MQ=255)
                                                                               cGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCttt                                                                     >  2:178603/1‑150 (MQ=255)
                                                                               cGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTAGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCttt                                                                     >  1:508910/1‑150 (MQ=255)
                                                                                gAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTc                                                                    <  1:1459415/150‑1 (MQ=255)
                                                                                gAAGAGACGCAGAGCGTGGCACGGGCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTc                                                                    <  1:630461/150‑1 (MQ=255)
                                                                                      aCGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTg                                                              <  1:178603/150‑1 (MQ=255)
                                                                                       cGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGc                                                             >  2:741456/1‑150 (MQ=255)
                                                                                       cGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGc                                                             <  2:846984/150‑1 (MQ=255)
                                                                                       cGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGc                                                             >  1:401440/1‑150 (MQ=255)
                                                                                       cGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGc                                                             >  2:150290/1‑150 (MQ=255)
                                                                                       cGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGc                                                             >  1:487448/1‑150 (MQ=255)
                                                                                         cAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCcg                                                           >  2:258072/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                                                        <  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)
                                                                                                      ggTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgt                                              >  1:347243/1‑150 (MQ=255)
                                                                                                       gtcgtcTTCCGGTCTATTCACCGTTTCGACGCGCTGATCGCCAAAGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgtt                                             <  1:146607/150‑1 (MQ=255)
                                                                                                       gtcgtcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTCATGGTGCCGTCAGATGAACGCTGTGCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTttgtt                                             <  1:1115027/150‑1 (MQ=255)
                                                                                                         cgtcTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACGTTGTTGa                                           >  1:153944/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)
                                                                                                                cGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGt                                    <  2:1074209/150‑1 (MQ=255)
                                                                                                                cGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGt                                    <  2:1300702/150‑1 (MQ=255)
                                                                                                                  ggATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGtca                                  <  2:1337686/150‑1 (MQ=255)
                                                                                                                       taCACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGtcatctca                             >  1:827701/1‑150 (MQ=255)
                                                                                                                       aacacCGTTTCGACGCGCTGATCGCCAACGTCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGtcatctca                             <  1:1276833/149‑1 (MQ=255)
                                                                                                                        acacCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGtcatctcat                            >  2:1232277/1‑150 (MQ=255)
                                                                                                                         cacCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGtcatctcatc                           >  1:1176656/1‑150 (MQ=255)
                                                                                                                          acCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCtcatca                          <  2:1355370/150‑1 (MQ=255)
                                                                                                                                tCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTg                    <  2:1135055/150‑1 (MQ=255)
                                                                                                                                 cGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGt                   >  1:1120929/1‑150 (MQ=255)
                                                                                                                                      gggcgATCGCCGTTGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCaa                         <  2:2131/135‑1 (MQ=255)
                                                                                                                                      gggcgATCGCCGTTGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCaa                         >  1:2131/5‑139 (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            >  1:986516/1‑150 (MQ=255)
                                                                                                                                        cTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCgcagca            >  2:1549650/1‑150 (MQ=255)
                                                                                                                                         tGATCGCCAACGGCGATCGCCCATTTAATGGAGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTATCAGCaa           >  2:437022/1‑150 (MQ=255)
                                                                                                                                          gATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGAGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAAt          >  1:13659/1‑150 (MQ=255)
                                                                                                                                                cAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTACGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGAt    <  2:300529/150‑1 (MQ=255)
                                                                                                                                                 aaCGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATa   <  1:891482/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TTTCCGACACCCGCAGGTTGCGGTTAAAGCCCTGCTGGGCGGTGGAACAGAATTTACACTCCAGCGCACACCCCACCTGCGAAGAGACGCAGAGCGTGGCACGGTCGTCTTCCGGGATATACACCGTTTCGACGCGCTGATCGCCAACGGCGATCGCCCATTTAATGGTGCCGTCAGATGAACGCTGTTCTTCAACCACTTCCGGTGCGCGGATTTCCGCCACCTCTTTCAGTTTGCCGCGCAACACTTTGTTGATGTCGGTCATCTCATCAAAGTTGTCGCAGCAATAGTGATAC  >  NC_000913/2643842‑2644137

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