Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,557,197 G→T 5.7% P281Q (CCG→CAG)  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,1970GT 5.7% 372.6 / 9.4 106P281Q (CCG→CAG) ddpFD,D‑dipeptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (50/50);  new base T (4/2);  total (54/52)
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.42e-01

Only 100 of 113 total aligned reads displayed.

ATCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGT  >  NC_000913/1557050‑1557345
                                                                                                                                                   |                                                                                                                                                    
aTCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCgg                                                                                                                                                     <  4:325365/149‑1 (MQ=255)
aTCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCgg                                                                                                                                                     <  3:520129/149‑1 (MQ=255)
aTCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCgg                                                                                                                                                     <  3:107354/149‑1 (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)
       aGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGt                                                                                                                                              >  4:368244/1‑149 (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                                                                                                                                                   >  1:417219/1‑135 (MQ=255)
                ggATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGc                                                                                                                                                   <  2:417219/135‑1 (MQ=255)
                ggATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGc                                                                                                                                                   <  4:612454/135‑1 (MQ=255)
                ggATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGc                                                                                                                                                   >  3:612454/1‑135 (MQ=255)
                   tGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGaaaaaa                                                                                                                                  <  2:519661/149‑1 (MQ=255)
                   tGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGaaaaaa                                                                                                                                  <  1:35679/149‑1 (MQ=255)
                   tGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGaaaaaa                                                                                                                                  <  2:894450/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)
                                       tGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCg                                                                                                              <  4:211816/149‑1 (MQ=255)
                                          cTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGtt                                                                                                           <  2:654098/149‑1 (MQ=255)
                                           tGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGttt                                                                                                          >  4:10798/1‑149 (MQ=255)
                                           tGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGttt                                                                                                          >  4:435034/1‑149 (MQ=255)
                                             aTGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTcc                                                                                                        >  1:167299/1‑149 (MQ=255)
                                              tGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTcc                                                                                                        <  1:526468/148‑1 (MQ=255)
                                              tGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTcc                                                                                                        >  2:526468/1‑148 (MQ=255)
                                                 gttgttATCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTTGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCa                                                                                                                     >  2:127018/1‑132 (MQ=255)
                                                 gttgttAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGaaa                                                                                                                                     >  1:735675/1‑116 (MQ=255)
                                                 gttgttAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGaaa                                                                                                                                     <  2:735675/116‑1 (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:402367/1‑149 (MQ=255)
                                                        gCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCg                                                                                             >  3:686141/1‑149 (MQ=255)
                                                         cATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGt                                                                                            >  3:325560/1‑149 (MQ=255)
                                                          aTCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGt                                                                                            >  2:620135/1‑148 (MQ=255)
                                                          aTCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGt                                                                                            <  1:620135/148‑1 (MQ=255)
                                                           tctcTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTCTTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGttt                                                                                          >  1:203757/1‑149 (MQ=255)
                                                           tctcTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGttt                                                                                          >  3:373373/1‑149 (MQ=255)
                                                            ctctACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGtttt                                                                                         <  1:413762/149‑1 (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)
                                                                           gaGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCAttc                                                                          >  2:173342/2‑149 (MQ=255)
                                                                              cACCGGAGCTCACGTCCGTCGTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCAttcttc                                                                       <  3:724599/149‑1 (MQ=255)
                                                                              cACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCAttcttc                                                                       <  3:615693/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                                                             <  1:832523/149‑1 (MQ=255)
                                                                                        cACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGt                                                             <  2:457872/149‑1 (MQ=255)
                                                                                        cACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGt                                                             <  1:737767/149‑1 (MQ=255)
                                                                                        cACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGt                                                             <  2:167299/149‑1 (MQ=255)
                                                                                         aCGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGtt                                                            <  3:10798/149‑1 (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)
                                                                                                     ccTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCgg                                                <  4:373373/149‑1 (MQ=255)
                                                                                                     ccTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCgg                                                >  3:757221/1‑149 (MQ=255)
                                                                                                       ttATTGCTAATCATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCgggg                                              >  3:514788/1‑149 (MQ=255)
                                                                                                        tATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGggga                                             <  4:688275/149‑1 (MQ=255)
                                                                                                        tATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGggga                                             <  2:692219/149‑1 (MQ=255)
                                                                                                         aTTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGgggag                                            >  3:264594/1‑149 (MQ=255)
                                                                                                          ttGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGgggagg                                           <  1:469708/149‑1 (MQ=255)
                                                                                                          ttGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGgggagg                                           <  3:776905/149‑1 (MQ=255)
                                                                                                          ttGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGgggagg                                           >  1:846378/1‑149 (MQ=255)
                                                                                                           tGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCGCATTCTTCCTCCAGCGGTTTATCAATGGCGgggaggg                                          >  3:574160/1‑149 (MQ=255)
                                                                                                             cTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGaa                                        >  1:606456/1‑149 (MQ=255)
                                                                                                               aaTGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATc                                      <  4:7754/149‑1 (MQ=255)
                                                                                                                  gATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAg                                   <  2:463555/149‑1 (MQ=255)
                                                                                                                  gATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAg                                   >  4:798247/1‑149 (MQ=255)
                                                                                                                        cGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCaataa                             >  4:750914/1‑149 (MQ=255)
                                                                                                                          gACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATc                           >  4:612989/1‑149 (MQ=255)
                                                                                                                            cTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCgg                         <  3:659770/149‑1 (MQ=255)
                                                                                                                            cTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCgg                         <  1:249296/149‑1 (MQ=255)
                                                                                                                             ttCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCggg                        >  3:389289/1‑149 (MQ=255)
                                                                                                                                  aTCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTAt                   >  1:738164/1‑149 (MQ=255)
                                                                                                                                  aTCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTAt                   <  2:380609/149‑1 (MQ=255)
                                                                                                                                  aTCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTAt                   >  1:692460/1‑149 (MQ=255)
                                                                                                                                  aTCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTAt                   >  4:32684/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)
                                                                                                                                          gTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGc           >  4:656926/1‑149 (MQ=255)
                                                                                                                                            tGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAg         <  2:846378/149‑1 (MQ=255)
                                                                                                                                                   gggCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCgg                                                                                                             >  2:491157/1‑42 (MQ=255)
                                                                                                                                                   gggCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGt  <  3:147159/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCGTAAAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTATGGATGTGCAGGTGCCGT  >  NC_000913/1557050‑1557345

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