Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,596,356 Δ1 bp 5.3% coding (382/864 nt) ypfJ ← putative neutral zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,596,3560T. 5.3% 247.3 / 6.4 75coding (382/864 nt)ypfJputative neutral zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base T (33/38);  new base . (2/2);  total (35/40)
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.63e-01

AACAGTTTCTGCACATGATGACCGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGC  >  NC_000913/2596208‑2596503
                                                                                                                                                    |                                                                                                                                                   
aaCAGTTTCTGCACATGATGACCGACTTCATGGGCGATAACGTACCCCTGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCAt                                                                                                                                                     >  3:131266/1‑148 (MQ=255)
   aGTTTCTGCACATGATGACCGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTAt                                                                                                                                                  <  2:220454/149‑1 (MQ=255)
    gTTTCTGCACATGATGACCGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATg                                                                                                                                                 <  7:369848/149‑1 (MQ=255)
    gTTTCTGCACATGATGACCGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATg                                                                                                                                                 <  7:296248/149‑1 (MQ=255)
           cacaTGATGACCGACTTCATGGGCGATAACGTACCCCTTGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGgc                                                                                                                                          <  1:302647/149‑1 (MQ=255)
               tgatgaCCGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGAc                                                                                                                                              >  7:329257/1‑141 (MQ=255)
               tgatgaCCGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGAc                                                                                                                                              <  8:329257/141‑1 (MQ=255)
                gatgaCCGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGgcccgcc                                                                                                                                     >  7:202838/1‑149 (MQ=255)
                 atgaCCGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGgcccgccc                                                                                                                                    <  8:119843/149‑1 (MQ=255)
                  tgaCCGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGTACTGGcccgcccc                                                                                                                                   <  4:124403/149‑1 (MQ=255)
                   gaCCGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGcccgccccg                                                                                                                                  >  6:201073/1‑149 (MQ=255)
                   gaCCGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGcccgccccg                                                                                                                                  >  7:39415/1‑149 (MQ=255)
                   gaCCGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGCGCAATAGAACGGCCCCATTATGGACTGGcccgccccg                                                                                                                                  >  4:105095/1‑149 (MQ=255)
                     ccGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGACCGGCCCCATTATGGACTGGCCCGCCCCGCa                                                                                                                                >  4:236356/1‑149 (MQ=255)
                      cGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAg                                                                                                                               >  3:36352/1‑149 (MQ=255)
                       gACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGc                                                                                                                              >  3:168890/1‑149 (MQ=255)
                       gACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGc                                                                                                                              <  5:273917/149‑1 (MQ=255)
                       gACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATATGGACTGGCCCGCCCCGCAGcc                                                                                                                             >  4:4560/1‑149 (MQ=255)
                          ttCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCgg                                                                                                                           >  7:261987/1‑149 (MQ=255)
                                  cGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGt                                                                                                                    >  1:123436/1‑148 (MQ=255)
                                         gTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCAc                                                                                                            <  3:236356/149‑1 (MQ=255)
                                              ccTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTa                                                                                                        >  4:108379/1‑148 (MQ=255)
                                              ccTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTAc                                                                                                       <  1:260515/149‑1 (MQ=255)
                                                tGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCACATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACAt                                                                                                     >  4:144233/1‑149 (MQ=255)
                                                ctggCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACAt                                                                                                     <  4:131266/147‑1 (MQ=255)
                                                 gggCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGAc                                                                                                                                              <  2:261147/107‑1 (MQ=255)
                                                  tgCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCTGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACAt                                                                                                     <  8:191136/146‑1 (MQ=255)
                                                  ggCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACAt                                                                                                     <  8:261987/147‑1 (MQ=255)
                                                  ggCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACAt                                                                                                     <  8:39415/147‑1 (MQ=255)
                                                   gCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGAc                                                                                                  >  2:414871/1‑149 (MQ=255)
                                                   gCAAAATCGACATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCTGGCAATATAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGa                                                                                                   <  3:105095/148‑1 (MQ=255)
                                                    cAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGAc                                                                                                  >  6:175759/1‑148 (MQ=255)
                                                           gCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCg                                                                                          <  2:123436/149‑1 (MQ=255)
                                                            ccATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGgc                                                                                                                                          <  2:222456/100‑1 (MQ=255)
                                                            ccATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGgc                                                                                                                                          >  1:222456/1‑100 (MQ=255)
                                                            ccATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGt                                                                                        <  3:4560/149‑1 (MQ=255)
                                                              aTCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGtt                                                                                       <  8:202838/149‑1 (MQ=255)
                                                               tCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTg                                                                                      >  8:3904/1‑149 (MQ=255)
                                                                  gcgcCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCACCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTg                                                                                   <  2:112557/149‑1 (MQ=255)
                                                                    gcCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGa                                                                                  >  6:310697/1‑148 (MQ=255)
                                                                    gcCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGAt                                                                                 <  8:18821/149‑1 (MQ=255)
                                                                     ccAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGa                                                                                  >  5:6497/1‑147 (MQ=255)
                                                                      cAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGTGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAg                                                                               <  5:201073/149‑1 (MQ=255)
                                                                         gTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTc                                                                            >  1:13394/1‑149 (MQ=255)
                                                                          tttGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCt                                                                           >  6:55175/1‑149 (MQ=255)
                                                                                tttCATGTAATCATATAAGGAGAGATCGATATAAACCGTGCCATCCTACGGGCAATAGAACGGCCCCATTATGTACTTGACCGCCCCGCAGCCGGTAAGCGTCATTACACGGTACATGACAAGCTTCCGTTGCTGATATGTCTTAcca                                                                      >  4:32390/1‑147 (MQ=255)
                                                                                    aTGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCAtctt                                                                 >  1:377831/1‑149 (MQ=255)
                                                                                    aTGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCAtct                                                                  <  4:168890/148‑1 (MQ=255)
                                                                                       tcatcaTAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGGCTTACCCATCTtctc                                                              <  3:144233/149‑1 (MQ=255)
                                                                                          tcatAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGaa                                                           <  2:377831/149‑1 (MQ=255)
                                                                                            atAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAAc                                                          <  1:414871/148‑1 (MQ=255)
                                                                                             tAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtg                                                        >  6:8290/1‑149 (MQ=255)
                                                                                             tAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtg                                                        <  6:6497/149‑1 (MQ=255)
                                                                                               gAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTAc                                                                        <  6:79698/131‑1 (MQ=255)
                                                                                               gAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTAc                                                                        >  5:79698/1‑131 (MQ=255)
                                                                                               gAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtgt                                                       <  2:13394/148‑1 (MQ=255)
                                                                                                   gagagaTCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTcc                                                  <  7:3904/149‑1 (MQ=255)
                                                                                                         tCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTACACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTTTCCCCAGgt                                            >  2:314546/1‑149 (MQ=255)
                                                                                                               taAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtgt                                                       >  2:328323/1‑132 (MQ=255)
                                                                                                               taAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtgt                                                       <  1:328323/132‑1 (MQ=255)
                                                                                                                aaaCCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTc                                                                            <  6:338634/110‑1 (MQ=255)
                                                                                                                aaaCCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTc                                                                            >  5:338634/1‑110 (MQ=255)
                                                                                                                          aTCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCa                           <  5:175759/149‑1 (MQ=255)
                                                                                                                          aTCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCa                           <  7:200879/149‑1 (MQ=255)
                                                                                                                                  ggCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCAcc                   >  7:25635/1‑149 (MQ=255)
                                                                                                                                  ggCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCAcc                   <  3:108379/149‑1 (MQ=255)
                                                                                                                                  ggCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCAcc                   >  1:23059/1‑149 (MQ=255)
                                                                                                                                     aaTATAACTGCACAATTATGGACTGGCCCTCCCAGCAGCCGGTACGCGTCATTCCACGTTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTATTCCGTGGTTGCCAGAATCACCGa                 <  3:32390/148‑1 (MQ=255)
                                                                                                                                     aaTAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAg                >  2:89670/1‑149 (MQ=255)
                                                                                                                                       tAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGt              >  2:241366/1‑149 (MQ=255)
                                                                                                                                         gAACGGCCCCATATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGaa           >  3:219518/1‑149 (MQ=255)
                                                                                                                                              gCCCCATATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTg      <  4:219518/149‑1 (MQ=255)
                                                                                                                                                cccATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGc     <  5:55175/149‑1 (MQ=255)
                                                                                                                                                cccATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGc     <  5:310697/149‑1 (MQ=255)
                                                                                                                                                   aTTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGc  >  7:66062/1‑149 (MQ=255)
                                                                                                                                                   aTTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGc  >  5:322974/1‑149 (MQ=255)
                                                                                                                                                    ttATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCAtc                                                                   <  5:262556/83‑1 (MQ=255)
                                                                                                                                                    ttATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCAtc                                                                   >  6:262556/1‑83 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACAGTTTCTGCACATGATGACCGACTTCATGGGCGATAACGTACCCCTGGGCAAAATCGCCATCCGCGCCAAGTTTGTCTTTCATGTCATCATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGC  >  NC_000913/2596208‑2596503

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