Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,642,358 T→A 5.5% intergenic (+27/+22) dtpB → / ← rsmJ dipeptide and tripeptide permease B/16S rRNA m(2)G1516 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,642,3580TA 5.5% 2501.9 / 88.4 687intergenic (+27/+22)dtpB/rsmJdipeptide and tripeptide permease B/16S rRNA m(2)G1516 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (336/313);  new base A (26/12);  total (362/325)
Fisher's exact test for biased strand distribution p-value = 6.47e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

Only 100 of 692 total aligned reads displayed.

ACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGC  >  NC_000913/3642209‑3642506
                                                                                                                                                     |                                                                                                                                                    
aCGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATATCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGt                                                                                                                                                      >  1:3665678/1‑150 (MQ=255)
  gTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTa                                                                                                                                                    >  1:2379115/1‑150 (MQ=255)
  gTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTa                                                                                                                                                    >  1:6593379/1‑150 (MQ=255)
   ttGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTAt                                                                                                                                                   <  1:1346142/150‑1 (MQ=255)
    tGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATc                                                                                                                                                  <  1:2284746/150‑1 (MQ=255)
    tGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATc                                                                                                                                                  <  2:4900519/150‑1 (MQ=255)
         gTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTAc                                                                                                                                             <  1:6393197/150‑1 (MQ=255)
         gTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTAc                                                                                                                                             <  1:106374/150‑1 (MQ=255)
          tCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGATAACGCTAcc                                                                                                                                            <  2:6758172/150‑1 (MQ=255)
              taCCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgt                                                                                                                                        >  2:6153406/1‑150 (MQ=255)
                 cAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtttt                                                                                                                                     >  1:7989582/1‑150 (MQ=255)
                  aaCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgtttta                                                                                                                                    >  1:3365428/1‑150 (MQ=255)
                   aCGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttag                                                                                                                                   <  2:6250640/150‑1 (MQ=255)
                     gtgtTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagtt                                                                                                                                 <  1:6702581/150‑1 (MQ=255)
                      tgtTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagttt                                                                                                                                <  1:7615649/150‑1 (MQ=255)
                        tttGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagtttta                                                                                                                              <  2:2974435/150‑1 (MQ=255)
                        tttGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagtttta                                                                                                                              >  1:5434010/1‑150 (MQ=255)
                        tttGGTAAGATTGGTCAGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTgttttagtttta                                                                                                                              >  2:7140387/1‑150 (MQ=255)
                           ggTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTacac                                                                                                                           >  2:6904426/1‑150 (MQ=255)
                              aaGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTgg                                                                                                                        >  2:3312417/1‑150 (MQ=255)
                               aGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGc                                                                                                                       <  2:2728981/150‑1 (MQ=255)
                                     gtctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCt                                                                                                                 <  1:5873051/150‑1 (MQ=255)
                                     gtctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCt                                                                                                                 <  1:6827024/150‑1 (MQ=255)
                                     gtctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCt                                                                                                                 <  1:6994375/150‑1 (MQ=255)
                                      tctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTg                                                                                                                >  1:5739032/1‑150 (MQ=255)
                                      tctggtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTg                                                                                                                >  1:1242837/1‑150 (MQ=255)
                                          gtcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATa                                                                                                            >  1:1980157/1‑150 (MQ=255)
                                           tcACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAg                                                                                                           >  1:3550901/1‑150 (MQ=255)
                                             aCGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGat                                                                                                         <  1:4530956/150‑1 (MQ=255)
                                                cTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCATAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATc                                                                                                      >  2:8068628/1‑150 (MQ=255)
                                                    gCGTGGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCTAAg                                                                                                  <  2:739605/150‑1 (MQ=255)
                                                         gCAGGAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATg                                                                                             <  2:5336000/150‑1 (MQ=255)
                                                                tGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTa                                                                                      <  1:784278/150‑1 (MQ=255)
                                                                  aTGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTATTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGt                                                                                    >  1:2809579/1‑150 (MQ=255)
                                                                   tGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTg                                                                                   >  2:5984929/1‑150 (MQ=255)
                                                                    gCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAgaga                                                                                  >  2:6532099/1‑150 (MQ=255)
                                                                    gCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGa                                                                                  >  2:1136904/1‑150 (MQ=255)
                                                                    gCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGa                                                                                  <  1:6064930/150‑1 (MQ=255)
                                                                     cTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGAc                                                                                 >  2:8189759/1‑150 (MQ=255)
                                                                      tGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAGTTATTCTTGCAGAAGGCAGGGTAGCGTTATCGCTGCCCTGTTTTTGTTTTACACTGGCGTAGCTGCATAGATATCAAAGCGATGCCCGTGAGTGACt                                                                                >  2:4426846/1‑150 (MQ=255)
                                                                      tGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACt                                                                                >  1:5018037/1‑150 (MQ=255)
                                                                        ttGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTAc                                                                              >  1:83876/1‑150 (MQ=255)
                                                                        ttGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTAc                                                                              >  2:7647842/1‑150 (MQ=255)
                                                                          gATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTGGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACcg                                                                            >  1:3005103/1‑150 (MQ=255)
                                                                          gATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACcg                                                                            <  2:3933642/150‑1 (MQ=255)
                                                                             ggTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGt                                                                         >  2:5436345/1‑150 (MQ=255)
                                                                             ggTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCGGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGt                                                                         >  2:5753836/1‑150 (MQ=255)
                                                                              gTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGtt                                                                        >  2:5034656/1‑150 (MQ=255)
                                                                                  cGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGGTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCGGTAGTGACTACCGCGTTtggc                                                                    >  2:4092033/1‑150 (MQ=255)
                                                                                   gTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcg                                                                   >  2:457212/1‑150 (MQ=255)
                                                                                   gTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcg                                                                   >  2:1996058/1‑150 (MQ=255)
                                                                                      gCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcgtgg                                                                <  2:6485207/150‑1 (MQ=255)
                                                                                      gCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTtggcgtgg                                                                <  1:7768052/150‑1 (MQ=255)
                                                                                      gCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGGTTTAGTTTTACACTGGCGTACCTGCAGAGATATCAAAGCGATGCCCTGTAGTGACTACCGCGGTtggcgtgg                                                                >  2:7852148/1‑150 (MQ=255)
                                                                                             cGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACAtt                                                         >  1:3846952/1‑150 (MQ=255)
                                                                                              gCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTc                                                        <  1:5498741/150‑1 (MQ=255)
                                                                                                 tGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGcc                                                     >  1:2510002/1‑150 (MQ=255)
                                                                                                  gATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCa                                                    <  2:2536051/150‑1 (MQ=255)
                                                                                                  gATAGCGACGCCGGAAAGCCATTAATTATTCGTGCAGCAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCa                                                    <  2:5257007/150‑1 (MQ=255)
                                                                                                    ttGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGc                                                  >  2:7540158/1‑150 (MQ=255)
                                                                                                       cGACGCCGGAAAGCCATTACTTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGt                                               <  2:1272407/150‑1 (MQ=255)
                                                                                                       cGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGt                                               <  1:1721165/150‑1 (MQ=255)
                                                                                                       cGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGt                                               <  1:564457/150‑1 (MQ=255)
                                                                                                        gACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTg                                              <  1:7160184/150‑1 (MQ=255)
                                                                                                        gACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGATTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTg                                              <  1:4700323/150‑1 (MQ=255)
                                                                                                          cGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGgc                                            <  2:1755471/150‑1 (MQ=255)
                                                                                                           gCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGgcg                                           >  2:4267504/1‑150 (MQ=255)
                                                                                                            ccGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGgcgc                                          >  2:7755325/1‑150 (MQ=255)
                                                                                                            ccGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGgcgc                                          <  2:1462786/150‑1 (MQ=255)
                                                                                                                 aaGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGt                                     >  2:240011/1‑150 (MQ=255)
                                                                                                                      attaattaTTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCggg                                >  2:3632809/1‑150 (MQ=255)
                                                                                                                      attaattaTTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCggg                                >  2:6367089/1‑150 (MQ=255)
                                                                                                                        taattaTTCTTGCAGAAAGCAGGGGAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGgcg                              <  1:8116567/150‑1 (MQ=255)
                                                                                                                            tattCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCttt                          <  2:3126920/150‑1 (MQ=255)
                                                                                                                              ttCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTac                        <  1:4442338/150‑1 (MQ=255)
                                                                                                                              ttCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTac                        <  1:1472920/150‑1 (MQ=255)
                                                                                                                              ttCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTac                        <  2:2260323/150‑1 (MQ=255)
                                                                                                                                cTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTacca                      >  2:7276687/1‑150 (MQ=255)
                                                                                                                                 ttGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTaccac                     >  1:6503876/1‑150 (MQ=255)
                                                                                                                                 ttGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTaccac                     >  2:1821068/1‑150 (MQ=255)
                                                                                                                                   gCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTAccacca                   <  1:2617276/150‑1 (MQ=255)
                                                                                                                                    cAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACcaccac                  >  2:284148/1‑150 (MQ=255)
                                                                                                                                      gcaaGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCAcgc                <  1:632792/148‑1 (MQ=255)
                                                                                                                                        aaacAGGGTAGCGATAACGCTACCCTGCTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGt              >  2:1624114/4‑150 (MQ=255)
                                                                                                                                         aacAGGGTAGCGATAACGCTACCCTGCTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGGAGTCCGGGCGCTTTCCCACCAcgcggt             >  1:4946070/3‑148 (MQ=255)
                                                                                                                                            aGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTgg          <  1:7287450/150‑1 (MQ=255)
                                                                                                                                            aGGGTAGCGATAACGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTgg          >  2:3408850/1‑150 (MQ=255)
                                                                                                                                             gggTAGCGATAACGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGt         >  1:1833026/1‑150 (MQ=255)
                                                                                                                                              ggTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTc        <  2:1350386/150‑1 (MQ=255)
                                                                                                                                              ggTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTc        <  2:2164490/150‑1 (MQ=255)
                                                                                                                                              ggTAGCGATAACGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTc        >  2:3506856/1‑150 (MQ=255)
                                                                                                                                              ggTAGCGATAACGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTc        >  2:1059775/1‑150 (MQ=255)
                                                                                                                                                tAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGc      <  2:1893714/150‑1 (MQ=255)
                                                                                                                                                tAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGc      <  1:5173256/150‑1 (MQ=255)
                                                                                                                                                tAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGc      <  1:4260307/150‑1 (MQ=255)
                                                                                                                                                  gCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCa    >  2:5600130/1‑150 (MQ=255)
                                                                                                                                                  gCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCa    <  2:3056136/150‑1 (MQ=255)
                                                                                                                                                  gCGTTATCGCGACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCa    <  2:4010836/150‑1 (MQ=255)
                                                                                                                                                   cGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAg   >  1:3950692/1‑150 (MQ=255)
                                                                                                                                                    gTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGc  >  1:2193351/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACGTTGCCCGTCTATACCAACGTGTTTGGTAAGATTGGTCTGGTCACGCTGGGCGTTGCAGTAGTGATGCTGTTGATGGTGCCGTGGCTGAAACGCATGATTGCGACGCCGGAAAGCCATTAATTATTCTTGCAGAAAGCAGGGTAGCGTTATCGCTACCCTGTTTTAGTTTTACACTGGCGTACCTGCATAGATATCAAAGCGATGCCCTTTAGTGACTACCGCGTTTGGCGTGGCGACATTCGCCAGCGGTGGCGCGTAGTCCGGGCGCTTTACCACCACGCGTTTGGTCGCCAGC  >  NC_000913/3642209‑3642506

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