Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,557,337 A→T 5.3% P234P (CCT→CCA)  ddpF ← D,D‑dipeptide ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,557,3370AT 5.3% 241.7 / 3.3 76P234P (CCT→CCA) ddpFD,D‑dipeptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (33/39);  new base T (2/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

CCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGC  >  NC_000913/1557192‑1557478
                                                                                                                                                 |                                                                                                                                             
ccAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGt                                                                                                                                            <  3:298574/149‑1 (MQ=255)
 cAGCGGGCAACGTTCGTAGAAAAAACAGTCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGCTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTg                                                                                                                                           >  4:138997/1‑149 (MQ=255)
  aGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGc                                                                                                                                          >  4:428041/1‑149 (MQ=255)
        cAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCAGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAg                                                                                                                                    <  5:345824/149‑1 (MQ=255)
          aCGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCa                                                                                                                                  <  6:131700/149‑1 (MQ=255)
           cGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCAc                                                                                                                                 <  3:197456/149‑1 (MQ=255)
                gTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACctgct                                                                                                                            >  5:460404/1‑149 (MQ=255)
                gTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACctgct                                                                                                                            >  6:48835/1‑149 (MQ=255)
                       aaaCAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTc                                                                                                                     <  1:407717/149‑1 (MQ=255)
                        aaCAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCg                                                                                                                    <  3:188267/149‑1 (MQ=255)
                             ccTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAg                                                                                                               >  1:411604/1‑149 (MQ=255)
                                 gCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGtt                                                                                                             <  3:249866/147‑1 (MQ=255)
                                  cGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCt                                                                                                           <  2:411604/148‑1 (MQ=255)
                                    gCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTAc                                                                                                         >  6:166759/1‑148 (MQ=255)
                                      aaCGGGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAAt                                                                                                      <  4:65094/149‑1 (MQ=255)
                                        cGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCt                                                                                                    <  5:34277/149‑1 (MQ=255)
                                        cGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCAATTCTTCCACCAGCGGTTTATCAATGGCGGGGAGGCAATCCAGCAATAATCGGTTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTTCCCGACGCCCAGTTCTACAATCt                                                                                                    >  2:48262/1‑149 (MQ=255)
                                           gCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGcc                                                                                                 >  1:383606/1‑149 (MQ=255)
                                                ttCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGccc                                                                                                >  1:321540/1‑145 (MQ=255)
                                                 tCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCgaga                                                                                            <  6:418598/148‑1 (MQ=255)
                                                  cccGGCAGATCCGTTTTACGTAATGCCCCTTCTTCCTACAGCGGTGTATCAATGGCGGGTAGGGAATCCAGCAATACTCGGGTGTATTGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGAt                                                                                           <  2:367557/148‑1 (MQ=255)
                                                  cccGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGAt                                                                                           <  3:179850/148‑1 (MQ=255)
                                                  cccGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGAt                                                                                           <  5:48835/148‑1 (MQ=255)
                                                         gATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCAcc                                                                                   <  1:48262/149‑1 (MQ=255)
                                                          aTCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCg                                                                                  >  5:36221/1‑149 (MQ=255)
                                                          aTCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCg                                                                                  >  4:318039/1‑149 (MQ=255)
                                                              gTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAc                                                                              <  1:173304/149‑1 (MQ=255)
                                                              gTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCAc                                                                              >  7:147222/1‑149 (MQ=255)
                                                                      tAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGc                                                                      >  7:506175/1‑149 (MQ=255)
                                                                        aTGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCt                                                                     <  3:428041/148‑1 (MQ=255)
                                                                        aTGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTc                                                                    <  2:236922/149‑1 (MQ=255)
                                                                         tGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCa                                                                   <  3:318039/149‑1 (MQ=255)
                                                                            ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATa                                                                                          <  3:334700/123‑1 (MQ=255)
                                                                            ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATa                                                                                          >  4:334700/1‑123 (MQ=255)
                                                                            ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCatat                                                                >  6:402021/1‑149 (MQ=255)
                                                                               ttcttcCTCCAGCGTTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGAc                                                             <  2:383606/149‑1 (MQ=255)
                                                                               ttcttcCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGAc                                                             <  2:321540/149‑1 (MQ=255)
                                                                                  ttcCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGCTCGCTCATATGACGTa                                                          >  5:90359/1‑149 (MQ=255)
                                                                                     ctccAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCa                                                       >  8:235945/1‑149 (MQ=255)
                                                                                        cAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATAGATCACCGCCACCCGATCGCTCATATGACGTATCAc                                                      <  4:2564/147‑1 (MQ=255)
                                                                                         aGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCTCTCATATGACGTATCACCga                                                   <  6:460404/149‑1 (MQ=255)
                                                                                          gCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCAcc                                                     >  2:72659/1‑146 (MQ=255)
                                                                                                  tCAATGTCGGGGAGGTAATCCAGCAATAATCGGGTGTATGGATGTGCAGCTGGCGTCAGCACCTGCTTCGCGTCGCCCAGTTCTACAATCTGGCCGAGATACATCACCGCCAACCGTTCGCTCATATGACGTATCTCCGAGACATtgt                                           <  8:231235/148‑1 (MQ=255)
                                                                                                    aaTGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgt                                           <  8:414450/146‑1 (MQ=255)
                                                                                                    aaTGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgt                                           >  7:414450/1‑146 (MQ=255)
                                                                                                        gCGGGGAGGGAATCCAGCAATAATCGGGCGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTGCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGaaa                                     <  8:186203/148‑1 (MQ=255)
                                                                                                          ggggAGGGAATCCAGCAATAATCGGGTGTATGTATGTGCAGGGGCCTTCAGCACCTGCTGCGCGTCACCCAGTGCTACAATATGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCa                                  <  6:100782/149‑1 (MQ=255)
                                                                                                          ggggAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACCTCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAAt                                    >  4:299098/1‑147 (MQ=255)
                                                                                                                ggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGcaca                             >  7:282999/1‑148 (MQ=255)
                                                                                                                  aaTCCATCAATAATCGGGGGTATGGATGTGCAGGTGCCGTCAGCTCCTGCTGCGAGTCGCCCAGTTCTACAATCTGCCCGAGATACATCAGCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAg                          <  6:90359/149‑1 (MQ=255)
                                                                                                                  aaTCCAGCAATAATCGGGTGTTTGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTATACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAg                          >  3:41249/1‑149 (MQ=255)
                                                                                                                  aaTCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAg                          >  6:275195/1‑149 (MQ=255)
                                                                                                                       aGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGAc                                                             >  4:391921/1‑109 (MQ=255)
                                                                                                                       aGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGAc                                                             <  3:391921/109‑1 (MQ=255)
                                                                                                                       aGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATc                                   >  6:120074/1‑135 (MQ=255)
                                                                                                                       aGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATc                                   <  5:120074/135‑1 (MQ=255)
                                                                                                                       aGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAg                     >  2:327855/1‑149 (MQ=255)
                                                                                                                       aGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAg                     >  1:56033/1‑149 (MQ=255)
                                                                                                                        gCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGc                    <  6:36221/149‑1 (MQ=255)
                                                                                                                        gCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGc                    >  6:75387/1‑149 (MQ=255)
                                                                                                                          aataatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGccc                  >  7:16089/1‑149 (MQ=255)
                                                                                                                              atCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGa              <  5:166759/149‑1 (MQ=255)
                                                                                                                              atCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGa              <  7:355451/149‑1 (MQ=255)
                                                                                                                              atCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGa              <  1:72659/149‑1 (MQ=255)
                                                                                                                               tCGGGTGTATGGATGTGCTGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGa                                       <  1:437650/123‑1 (MQ=255)
                                                                                                                               tCGGGTGTATGGATGTGCTGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGa                                       <  1:444591/123‑1 (MQ=255)
                                                                                                                               tCGGGTGTATGGATGTGCTGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGa                                       >  2:444591/1‑123 (MQ=255)
                                                                                                                               tCGGGTGTATGGATGTGCTGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGa                                       >  2:437650/1‑123 (MQ=255)
                                                                                                                                   gtgtATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTc         <  5:75387/149‑1 (MQ=255)
                                                                                                                                    tgtATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTcc        <  3:138997/149‑1 (MQ=255)
                                                                                                                                        tGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTg      >  8:256008/1‑147 (MQ=255)
                                                                                                                                        tGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTa    <  1:177213/149‑1 (MQ=255)
                                                                                                                                        tGGATGTGCAGGTGCCGTCACCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTa    <  8:147222/149‑1 (MQ=255)
                                                                                                                                         ggATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAg   >  7:476127/1‑149 (MQ=255)
                                                                                                                                         ggATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAg   >  8:407663/1‑149 (MQ=255)
                                                                                                                                          gATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGc  >  3:178882/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGC  >  NC_000913/1557192‑1557478

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: ‑