Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,557,194 A→C 5.5% L282R (CTG→CGG)  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,1940AC 5.5% 387.9 / 11.1 110L282R (CTG→CGG) ddpFD,D‑dipeptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (51/53);  new base C (4/2);  total (55/55)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

Only 100 of 111 total aligned reads displayed.

TGAAATCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAG  >  NC_000913/1557046‑1557338
                                                                                                                                                    |                                                                                                                                                
tGAAATCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCa                                                                                                                                                  <  3:73850/149‑1 (MQ=255)
    aTCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCgg                                                                                                                                              <  4:325365/149‑1 (MQ=255)
    aTCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCgg                                                                                                                                              <  3:107354/149‑1 (MQ=255)
    aTCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCgg                                                                                                                                              <  3:520129/149‑1 (MQ=255)
    aTCGTAAAGCAGTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCgg                                                                                                                                              >  3:693847/1‑149 (MQ=255)
       gTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCa                                                                                                                                           >  2:709477/1‑149 (MQ=255)
        tAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCaa                                                                                                                                          >  4:724599/1‑149 (MQ=255)
         aaaGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAAc                                                                                                                                         >  2:562592/1‑149 (MQ=255)
          aaGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCCGCTGGCAAc                                                                                                                                         >  4:377455/1‑148 (MQ=255)
          aaGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCCGCTGGCAAc                                                                                                                                         <  3:377455/148‑1 (MQ=255)
            gCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGtt                                                                                                                                      <  2:443025/149‑1 (MQ=255)
               ttGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGt                                                                                                                                   >  2:658455/1‑149 (MQ=255)
                tGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTa                                                                                                                                  <  3:846142/149‑1 (MQ=255)
                tGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGACCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTa                                                                                                                                  <  4:857948/149‑1 (MQ=255)
                   cgtaTGACGCGTCAGGCGCGTCAATTCCTGATGCGTTGTTAGCATCTCTACGGCGGCCGCCAGCACCGGAGGTCACGTCAGTCCTCCCTTATTGCTAATGATTGCAGGACTTCACATCCGTGAATTGCCATCGGGCAACGTTCGTAGaa                                                                                                                               >  4:548595/4‑149 (MQ=255)
                    ggATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGc                                                                                                                                            <  2:417219/135‑1 (MQ=255)
                    ggATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGc                                                                                                                                            >  1:417219/1‑135 (MQ=255)
                    ggATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGc                                                                                                                                            >  3:612454/1‑135 (MQ=255)
                    ggATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGc                                                                                                                                            <  4:612454/135‑1 (MQ=255)
                       tGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGaaaaaa                                                                                                                           <  2:894450/149‑1 (MQ=255)
                       tGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGaaaaaa                                                                                                                           <  2:557015/149‑1 (MQ=255)
                       tGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGaaaaaa                                                                                                                           <  2:519661/149‑1 (MQ=255)
                           gcgTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGc                                                                                                                       >  1:463555/1‑149 (MQ=255)
                                aGGCGCGTGAATGCCTGATGCGTTGTTAGCACCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGc                                                                                                                  >  1:210834/1‑149 (MQ=255)
                                 ggCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCg                                                                                                                 >  3:579169/1‑149 (MQ=255)
                                   cgcgTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGc                                                                                                               <  1:13044/149‑1 (MQ=255)
                                      gTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAAc                                                                                                            <  3:187682/149‑1 (MQ=255)
                                           tGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCg                                                                                                       <  4:211816/149‑1 (MQ=255)
                                              cTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGtt                                                                                                    <  2:654098/149‑1 (MQ=255)
                                               tGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGttt                                                                                                   >  4:435034/1‑149 (MQ=255)
                                               tGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGttt                                                                                                   >  4:10798/1‑149 (MQ=255)
                                                 aTGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTcc                                                                                                 >  1:167299/1‑149 (MQ=255)
                                                  tGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTcc                                                                                                 <  1:526468/148‑1 (MQ=255)
                                                     gttgttAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGaaa                                                                                                                              <  2:735675/116‑1 (MQ=255)
                                                     gttgttAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGaaa                                                                                                                              >  1:735675/1‑116 (MQ=255)
                                                     gttgttAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCa                                                                                                              <  1:127018/132‑1 (MQ=255)
                                                     gttgttAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGc                                                                                             >  1:289792/1‑149 (MQ=255)
                                                         ttAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGAt                                                                                         >  4:653270/1‑149 (MQ=255)
                                                         ttAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGAt                                                                                         >  4:730707/1‑149 (MQ=255)
                                                            gCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCg                                                                                      >  3:686141/1‑149 (MQ=255)
                                                            gCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCg                                                                                      >  3:402367/1‑149 (MQ=255)
                                                             cATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGt                                                                                     >  3:325560/1‑149 (MQ=255)
                                                              aTCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGt                                                                                     >  2:620135/1‑148 (MQ=255)
                                                               tctcTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTCTTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGttt                                                                                   >  1:203757/1‑149 (MQ=255)
                                                                ctctACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGtttt                                                                                  <  1:413762/149‑1 (MQ=255)
                                                                 tctACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTa                                                                                 >  3:539713/1‑149 (MQ=255)
                                                                  ctACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTAc                                                                                <  3:730707/149‑1 (MQ=255)
                                                                  ctACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTAc                                                                                <  1:381582/149‑1 (MQ=255)
                                                                   tACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACg                                                                               >  3:688275/1‑149 (MQ=255)
                                                                           ccgccAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGccc                                                                       >  3:64280/1‑149 (MQ=255)
                                                                            cgccAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCa                                                                      <  4:402367/149‑1 (MQ=255)
                                                                               gaGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCAttc                                                                   >  2:173342/2‑149 (MQ=255)
                                                                                  cACCGGAGCTCACGTCCGTCGTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCAttcttc                                                                <  3:724599/149‑1 (MQ=255)
                                                                                  cACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCAttcttc                                                                <  1:15246/149‑1 (MQ=255)
                                                                                           tCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCgg                                                       <  3:1203/149‑1 (MQ=255)
                                                                                            cACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGt                                                      <  1:562592/149‑1 (MQ=255)
                                                                                            cACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGt                                                      <  2:167299/149‑1 (MQ=255)
                                                                                            cACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGt                                                      <  2:457872/149‑1 (MQ=255)
                                                                                            cACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGt                                                      <  1:737767/149‑1 (MQ=255)
                                                                                            cACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGt                                                      <  1:832523/149‑1 (MQ=255)
                                                                                             aCGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGtt                                                     <  3:10798/149‑1 (MQ=255)
                                                                                                tccgtccTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGtttgt                                                  <  1:306285/149‑3 (MQ=255)
                                                                                                 ccgtccTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATc                                                 >  1:380609/1‑149 (MQ=255)
                                                                                                   gtccTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTa                                                                                 <  4:815239/115‑1 (MQ=255)
                                                                                                   gtccTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCaa                                               <  4:335623/149‑1 (MQ=255)
                                                                                                   gtccTCCCGTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTa                                                                                 >  3:815239/1‑115 (MQ=255)
                                                                                                    tcctccCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAAt                                              >  4:61303/1‑149 (MQ=255)
                                                                                                    tcctccCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAAt                                              <  1:832296/149‑1 (MQ=255)
                                                                                                       tccCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGc                                           >  4:18509/1‑149 (MQ=255)
                                                                                                         ccTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCgg                                         <  4:693847/149‑1 (MQ=255)
                                                                                                         ccTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCgg                                         >  3:757221/1‑149 (MQ=255)
                                                                                                         ccTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCgg                                         <  4:373373/149‑1 (MQ=255)
                                                                                                           ttATTGCTAATCATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCgggg                                       >  3:514788/1‑149 (MQ=255)
                                                                                                            tATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGggga                                      <  2:692219/149‑1 (MQ=255)
                                                                                                            tATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGggga                                      <  4:688275/149‑1 (MQ=255)
                                                                                                             aTTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGgggag                                     >  3:264594/1‑149 (MQ=255)
                                                                                                              ttGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGgggagg                                    >  1:846378/1‑149 (MQ=255)
                                                                                                              ttGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGgggagg                                    <  1:469708/149‑1 (MQ=255)
                                                                                                               tGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCGCATTCTTCCTCCAGCGGTTTATCAATGGCGgggaggg                                   >  3:574160/1‑149 (MQ=255)
                                                                                                                 cTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGaa                                 >  1:606456/1‑149 (MQ=255)
                                                                                                                   aaTGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATc                               <  4:7754/149‑1 (MQ=255)
                                                                                                                      gATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAg                            >  4:798247/1‑149 (MQ=255)
                                                                                                                      gATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAg                            <  2:463555/149‑1 (MQ=255)
                                                                                                                            cGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCaataa                      >  4:750914/1‑149 (MQ=255)
                                                                                                                              gACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATc                    >  4:612989/1‑149 (MQ=255)
                                                                                                                                cTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCgg                  <  1:249296/149‑1 (MQ=255)
                                                                                                                                cTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCgg                  <  3:659770/149‑1 (MQ=255)
                                                                                                                                 ttCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCggg                 >  3:389289/1‑149 (MQ=255)
                                                                                                                                      aTCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTAt            <  2:289792/149‑1 (MQ=255)
                                                                                                                                      aTCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTAt            <  2:380609/149‑1 (MQ=255)
                                                                                                                                      aTCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTAt            >  4:32684/1‑149 (MQ=255)
                                                                                                                                      aTCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTAt            >  1:738164/1‑149 (MQ=255)
                                                                                                                                      aTCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTAt            >  1:692460/1‑149 (MQ=255)
                                                                                                                                       tCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATg           >  1:73809/1‑149 (MQ=255)
                                                                                                                                        acgaacgTTGCCCGCTGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATgg          >  2:787540/7‑149 (MQ=255)
                                                                                                                                         cgaacgTTGCCCGCTGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGa         >  2:852650/6‑149 (MQ=255)
                                                                                                                                         cGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGa         <  1:686140/149‑1 (MQ=255)
                                                                                                                                          gaacgTTGCCCGCTGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGAt        >  4:192422/5‑149 (MQ=255)
                                                                                                                                          gTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGAt        <  1:321946/149‑1 (MQ=255)
                                                                                                                                                tGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAg  <  2:846378/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TGAAATCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAG  >  NC_000913/1557046‑1557338

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