Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,557,347 A→G 5.9% L231P (CTG→CCG)  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,3470AG 5.9% 205.6 / 3.3 68L231P (CTG→CCG) ddpFD,D‑dipeptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (27/36);  new base G (2/2);  total (30/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

AACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATT  >  NC_000913/1557201‑1557492
                                                                                                                                                  |                                                                                                                                                 
aaCGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGc                                                                                                                                                 >  7:254405/1‑149 (MQ=255)
 aCGTGCGTAGACAAAACAGCCTTGCGGCAACGTGCGGTTGCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCACTGGCGGTGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCa                                                                                                                                                <  2:309326/149‑1 (MQ=255)
   gTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCAcc                                                                                                                                              <  1:403208/149‑1 (MQ=255)
      cGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACctgc                                                                                                                                           <  2:338163/149‑1 (MQ=255)
       gTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACctgct                                                                                                                                          >  2:59825/1‑149 (MQ=255)
       gTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACctgc                                                                                                                                           <  3:80217/148‑1 (MQ=255)
           aaaaaaCAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCTGACGGCACctgcacatc                                                                                                                                      <  3:377001/149‑6 (MQ=255)
               aaCAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCGGCCCCTGCTGCGCGTCg                                                                                                                                  >  6:278209/1‑149 (MQ=255)
                   gCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCa                                                                                                                              <  1:161578/149‑1 (MQ=255)
                    ccTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAg                                                                                                                             >  6:383065/1‑149 (MQ=255)
                        gCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTCATGCCCATTCTTCCTCCAGCGGGTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTc                                                                                                                          <  4:114332/148‑1 (MQ=255)
                         cGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTa                                                                                                                        <  3:332812/149‑1 (MQ=255)
                           gCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACa                                                                                                                      <  6:43571/149‑1 (MQ=255)
                             aaCGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAAt                                                                                                                    <  5:17942/149‑1 (MQ=255)
                             aaCGTGCGGTTTCCCGGCAGAGCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGCTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGGGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAAt                                                                                                                    <  6:330016/149‑1 (MQ=255)
                              aCGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATc                                                                                                                   >  1:106536/1‑149 (MQ=255)
                                 ttcGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCTGGGAGGGAATCCAGCAATAATCTGGTGTATGGATGTGCAGGTGCCGTCAACACCTGCTGCGCGTCGCCCAGTTCTACAATCTGc                                                                                                                >  7:194939/3‑149 (MQ=255)
                                 tGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGc                                                                                                                <  1:653/149‑1 (MQ=255)
                                 tGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGc                                                                                                                >  2:389162/1‑149 (MQ=255)
                                 tGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGc                                                                                                                <  5:383065/149‑1 (MQ=255)
                                        tCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGGGGTTTATCAATGGCGGGGAGGGAATCCAGCATTAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCTCGTCGCCCAGTTCTACAATCTGCCCgaga                                                                                                          <  6:65554/148‑1 (MQ=255)
                                        tCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCgaga                                                                                                          <  3:309691/148‑1 (MQ=255)
                                        tCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCgaga                                                                                                          <  3:329280/148‑1 (MQ=255)
                                        tCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCgaga                                                                                                          <  1:389162/148‑1 (MQ=255)
                                         cccGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGCTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGGGTATGGATGTGCAGGGGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGAt                                                                                                         <  6:279702/148‑1 (MQ=255)
                                          ccGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATa                                                                                                        <  5:100437/148‑1 (MQ=255)
                                          ccGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATa                                                                                                        >  3:212054/1‑148 (MQ=255)
                                            ggCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACAt                                                                                                     <  7:71153/149‑1 (MQ=255)
                                            ggCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACAt                                                                                                     >  2:64230/1‑149 (MQ=255)
                                               ataTCCGTTTTACGTACTGCCCATTCTTCATCCAGCTGTTTATCAATGGCGTGTAGGGAATCCAGCAATAATCGGGTTTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCAc                                                                                                  <  5:213661/147‑1 (MQ=255)
                                                 ctCCGTTTTACGTAATGCCCATTCTTCCTCCAGAGGGTTATCAATGGCGGGGAGGGAATCCCGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCCGCACCTGCTGCGCGTCGCCCAGTTCGACAATCTGCACGCGATAAATCTCCg                                                                                                >  6:258503/2‑149 (MQ=255)
                                                    cGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCa                                                                                             <  7:372608/149‑1 (MQ=255)
                                                                   ccATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCatat                                                                              >  3:264438/1‑149 (MQ=255)
                                                                     aTTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGa                                                                            <  8:70676/149‑1 (MQ=255)
                                                                      ttcttcCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGAc                                                                           >  7:59078/1‑149 (MQ=255)
                                                                      ttcttcCTCCAGCGGTTAATCACTGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGAc                                                                           <  2:85892/149‑1 (MQ=255)
                                                                               cAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCCCCGCCACCCGATCGCTCATATGACGTATCAcc                                                                   >  6:181786/1‑148 (MQ=255)
                                                                                     tttATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACAt                                                            <  8:59078/149‑1 (MQ=255)
                                                                                         tCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATtgt                                                         <  1:386577/148‑1 (MQ=255)
                                                                                                 ggggAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCa                                                <  8:254405/149‑1 (MQ=255)
                                                                                                  gggagggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATAc                                                                                                       >  4:191542/1‑93 (MQ=255)
                                                                                                    gagggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGca                                             <  7:112798/149‑1 (MQ=255)
                                                                                                       ggaATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGcaca                                           >  3:270962/1‑148 (MQ=255)
                                                                                                             cAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACctgct                                                                                                                                          <  8:59486/47‑1 (MQ=255)
                                                                                                             cAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACctgct                                                                                                                                          >  7:59486/1‑47 (MQ=255)
                                                                                                             cAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGt                                      >  2:194877/1‑147 (MQ=255)
                                                                                                                cAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGcc                                 >  7:27890/1‑149 (MQ=255)
                                                                                                                cAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGcc                                 >  3:165347/1‑149 (MQ=255)
                                                                                                                  ataatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCg                               >  2:43496/1‑149 (MQ=255)
                                                                                                                   taatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCa                                    >  7:109772/1‑143 (MQ=255)
                                                                                                                   taatCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCa                                    <  8:109772/143‑1 (MQ=255)
                                                                                                                           tgtATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTcc                      <  1:194877/149‑1 (MQ=255)
                                                                                                                             tATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCt                     <  4:264438/148‑1 (MQ=255)
                                                                                                                             tATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTg                    >  5:242404/1‑149 (MQ=255)
                                                                                                                               tGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGt                   <  1:59825/148‑1 (MQ=255)
                                                                                                                                ggATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAg                 >  3:174270/1‑149 (MQ=255)
                                                                                                                                  aTGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCg               >  7:374910/1‑149 (MQ=255)
                                                                                                                                     tgCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGtt             <  1:64230/148‑1 (MQ=255)
                                                                                                                                     tgCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTa            <  4:165347/149‑1 (MQ=255)
                                                                                                                                       cAGGTGCTGACGGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTa                  >  7:146459/1‑141 (MQ=255)
                                                                                                                                       cAGGTGCTGACGGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTa                  <  8:146459/141‑1 (MQ=255)
                                                                                                                                       cAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTAc           <  4:212054/148‑1 (MQ=255)
                                                                                                                                        aGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTAc           >  6:229977/1‑147 (MQ=255)
                                                                                                                                         ggTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCAACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACCTATCACAGAGACATTGTGTGAAATCAGCAAATAGCTCAGCCCGTGATTTACCAGTAGCGTTACCAg        >  7:58722/1‑149 (MQ=255)
                                                                                                                                              cGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAAt   >  2:375747/1‑149 (MQ=255)
                                                                                                                                                tCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGAc                                                                           <  5:182758/75‑1 (MQ=255)
                                                                                                                                                tCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGAc                                                                           >  6:182758/1‑75 (MQ=255)
                                                                                                                                                tCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAAtt  >  2:281671/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGTCAGCACCTGCTGCGCGTCGCCCAGTTCTACAATCTGCCCGAGATACATCACCGCCACCCGATCGCTCATATGACGTATCACCGAGACATTGTGTGAAATCAGCACATAGGTCAGCCCGTGATTTTCCTGTAGCGTTACCAGTAAATT  >  NC_000913/1557201‑1557492

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