Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,970,817:1 +G 5.7% coding (685/1131 nt) wecB → UDP‑N‑acetyl glucosamine‑2‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,970,8171.G 5.7% 541.3 / 27.2 123coding (685/1131 nt)wecBUDP‑N‑acetyl glucosamine‑2‑epimerase
Reads supporting (aligned to +/- strand):  ref base . (65/51);  new base G (2/5);  total (67/56)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 125 total aligned reads displayed.

GCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGT  >  NC_000913/3970672‑3970962
                                                                                                                                                  |                                                                                                                                                 
gCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑Tgg                                                                                                                                                <  1:472082/149‑1 (MQ=255)
  gTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGc                                                                                                                                               <  3:625219/148‑1 (MQ=255)
  gTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCa                                                                                                                                              <  2:844534/149‑1 (MQ=255)
  gTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCa                                                                                                                                              <  1:548389/149‑1 (MQ=255)
     aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGa                                                                                                                                            >  4:419414/1‑148 (MQ=255)
     aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGAc                                                                                                                                           <  1:574638/149‑1 (MQ=255)
     aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGat                                                                                                                                           <  2:1033568/149‑2 (MQ=255)
     aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGat                                                                                                                                           <  2:505441/149‑2 (MQ=255)
     aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGa                                                                                                                                            <  4:846628/148‑1 (MQ=255)
      ccAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGatt                                                                                                                                          <  2:415460/149‑3 (MQ=255)
        aGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGc                                                                                                                                               >  4:991499/1‑142 (MQ=255)
            gatgaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGGAGACATCGcc                                                                                                                                     <  3:601631/148‑1 (MQ=255)
             atgaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGcca                                                                                                                                    >  4:900721/1‑148 (MQ=255)
                       aCAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCa                                                                                                                          >  3:223125/1‑148 (MQ=255)
                       aCAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAg                                                                                                                         <  2:431654/149‑1 (MQ=255)
                            cTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGc                                                                                                                                               <  1:506859/122‑1 (MQ=255)
                            cTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGc                                                                                                                                               >  2:506859/1‑122 (MQ=255)
                             tGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATc                                                                                                                   <  1:107052/149‑1 (MQ=255)
                              gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATc                                                                                                                   <  3:600558/148‑1 (MQ=255)
                              gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATc                                                                                                                   >  3:205599/1‑148 (MQ=255)
                                gTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAg                                                                                                                >  2:571680/1‑149 (MQ=255)
                                 ttCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAg                                                                                                                >  3:47018/1‑148 (MQ=255)
                                  tCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGa                                                                                                               >  3:258426/1‑148 (MQ=255)
                                  tCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGa                                                                                                               >  4:78317/1‑148 (MQ=255)
                                   cAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGAt                                                                                                              <  3:100935/148‑1 (MQ=255)
                                      aaCTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGt                                                                                                           <  4:130974/148‑1 (MQ=255)
                                          ggcggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATc                                                                                                      <  2:821693/149‑1 (MQ=255)
                                              gcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCgg                                                                                                   >  3:101468/1‑148 (MQ=255)
                                              gcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCgg                                                                                                   >  4:249385/1‑148 (MQ=255)
                                              gcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGt                                                                                                  >  2:1019809/1‑149 (MQ=255)
                                               cAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGt                                                                                                  >  4:702855/1‑148 (MQ=255)
                                                 aaTTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCa                                                                                               >  1:125287/1‑149 (MQ=255)
                                                  aTTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCTCCAGGACATCCAGATTGTCTATCCGGTGCAt                                                                                              >  2:476094/1‑149 (MQ=255)
                                                     aCCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCAtctc                                                                                           >  1:334621/1‑149 (MQ=255)
                                                       ccGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCaa                                                                                         >  2:716285/1‑149 (MQ=255)
                                                        cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCa                                                                                                                          <  4:456256/115‑1 (MQ=255)
                                                        cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCa                                                                                                                          >  3:456256/1‑115 (MQ=255)
                                                        cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCaa                                                                                         >  3:531565/1‑148 (MQ=255)
                                                         gTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAcc                                                                                       <  2:204735/149‑1 (MQ=255)
                                                          tttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAcc                                                                                       >  3:40430/1‑148 (MQ=255)
                                                             aTCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGaa                                                                                   <  2:49610/149‑1 (MQ=255)
                                                             aTCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGa                                                                                    <  3:900721/148‑1 (MQ=255)
                                                              tCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGaa                                                                                   <  4:581125/148‑1 (MQ=255)
                                                              tCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGaa                                                                                   >  4:123539/1‑148 (MQ=255)
                                                               cGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAAc                                                                                  >  3:515333/1‑148 (MQ=255)
                                                                  ccccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCa                                                                              <  1:170957/149‑1 (MQ=255)
                                                                   cccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCa                                                                              <  2:611746/148‑1 (MQ=255)
                                                                   cccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCa                                                                              >  1:611746/1‑148 (MQ=255)
                                                                    ccGATAGAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCaga                                                                            >  2:217436/1‑149 (MQ=255)
                                                                     cGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCaga                                                                            >  4:284783/1‑148 (MQ=255)
                                                                      gATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCagag                                                                           >  3:765725/1‑148 (MQ=255)
                                                                           aaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCgg                                                                     >  2:324346/1‑149 (MQ=255)
                                                                             aGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGt                                                                    <  4:101468/148‑1 (MQ=255)
                                                                               atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCaa                                                                 >  1:969313/1‑149 (MQ=255)
                                                                               atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCa                                                                  >  4:1018058/1‑148 (MQ=255)
                                                                               atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCa                                                                  >  4:988717/1‑148 (MQ=255)
                                                                                 gatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAAt                                                                >  3:967220/1‑148 (MQ=255)
                                                                                    tCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCgg                                                                                                   <  3:390687/110‑1 (MQ=255)
                                                                                    tCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCgg                                                                                                   >  4:390687/1‑110 (MQ=255)
                                                                                    tCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                             <  3:342186/148‑1 (MQ=255)
                                                                                    tCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGa                                                                                                                                            <  3:835963/69‑1 (MQ=38)
                                                                                       ggTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCAtt                                                          >  3:391428/1‑148 (MQ=255)
                                                                                            ccGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTgggg                                                    >  1:6866/1‑149 (MQ=255)
                                                                                               gTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAt                                                 <  2:68675/149‑1 (MQ=255)
                                                                                                      gCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaa                                          >  2:168404/1‑149 (MQ=255)
                                                                                                         tGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATg                                       >  1:772559/1‑149 (MQ=255)
                                                                                                         tGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATg                                       <  2:1045397/149‑1 (MQ=255)
                                                                                                          gagaGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACAt                                                                                                                    >  4:370990/1‑71 (MQ=255)
                                                                                                              gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgt                                               >  3:706219/1‑136 (MQ=255)
                                                                                                                ttCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTc                                 >  4:803932/1‑148 (MQ=255)
                                                                                                                  cGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTg                               <  3:382097/148‑1 (MQ=255)
                                                                                                                   ggTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTga                              <  4:967220/148‑1 (MQ=255)
                                                                                                                      cGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatcg                           <  4:846253/148‑1 (MQ=255)
                                                                                                                        tGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatcgatc                        <  1:244103/149‑1 (MQ=255)
                                                                                                                         ggCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATcc                       >  1:974995/1‑149 (MQ=255)
                                                                                                                            tttGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCGACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCa                    >  2:878035/1‑149 (MQ=255)
                                                                                                                             ttGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATg                                       >  3:281839/1‑129 (MQ=255)
                                                                                                                             ttGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATg                                       <  4:281839/129‑1 (MQ=255)
                                                                                                                             ttGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCa                    <  4:391428/148‑1 (MQ=255)
                                                                                                                             ttGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCa                    >  4:576646/1‑148 (MQ=255)
                                                                                                                             ttGAAGAAATCTGCCACGCGC‑TGGCAGACATCACCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCa                    >  3:669231/1‑148 (MQ=255)
                                                                                                                               gaagaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAgg                  <  3:755369/148‑1 (MQ=255)
                                                                                                                                aagaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAg                >  2:96273/1‑149 (MQ=255)
                                                                                                                                 agaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAg                >  4:95010/1‑148 (MQ=255)
                                                                                                                                  gaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCAtt                                                          >  3:968556/1‑105 (MQ=255)
                                                                                                                                  gaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCAtt                                                          <  4:968556/105‑1 (MQ=255)
                                                                                                                                   atgtCTGCCA‑GCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTa              >  3:744905/4‑148 (MQ=255)
                                                                                                                                   aaaTCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTa              <  4:531565/148‑1 (MQ=255)
                                                                                                                                     aTCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAttt           >  2:23776/1‑149 (MQ=255)
                                                                                                                                     aTCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAttt           >  1:819177/1‑149 (MQ=255)
                                                                                                                                      tCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgat                             <  2:313892/130‑1 (MQ=255)
                                                                                                                                      tCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgat                             >  1:313892/1‑130 (MQ=255)
                                                                                                                                        tGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTAc         <  3:249385/148‑1 (MQ=255)
                                                                                                                                         gCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTAcc        >  3:161994/1‑148 (MQ=255)
                                                                                                                                         gCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTAcc        >  3:334795/1‑148 (MQ=255)
                                                                                                                                          ccACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAt                                                 <  1:893855/106‑1 (MQ=255)
                                                                                                                                          ccACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAt                                                 >  2:893855/1‑106 (MQ=255)
                                                                                                                                           cACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                            <  4:779161/46‑1 (MQ=255)
                                                                                                                                               cgc‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGt  <  3:1018058/148‑1 (MQ=255)
                                                                                                                                               cgc‑TGGCAGACATCGCCACCAAGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGt  <  4:765725/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGT  >  NC_000913/3970672‑3970962

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑