Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,970,813 Δ1 bp 5.7% coding (681/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,8130C. 5.7% 441.8 / 18.7 123coding (681/1131 nt)wecBUDP‑N‑acetyl glucosamine‑2‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (65/51);  new base . (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.91e-01

Only 100 of 123 total aligned reads displayed.

GCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTT  >  NC_000913/3970672‑3970959
                                                                                                                                             |                                                                                                                                                   
gCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑Tgg                                                                                                                                             <  1:472082/149‑1 (MQ=255)
  gTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGc                                                                                                                                            <  3:625219/148‑1 (MQ=255)
  gTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCa                                                                                                                                           <  2:844534/149‑1 (MQ=255)
     aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGAc                                                                                                                                        <  1:574638/149‑1 (MQ=255)
     aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGat                                                                                                                                        <  2:1033568/149‑2 (MQ=255)
     aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGa                                                                                                                                         <  4:846628/148‑1 (MQ=255)
        aGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGc                                                                                                                                            >  4:991499/1‑142 (MQ=255)
        aGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGc                                                                                                                                            <  3:991499/142‑1 (MQ=255)
            gatgaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGGAGACATCGcc                                                                                                                                  <  3:601631/148‑1 (MQ=255)
             atgaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGcca                                                                                                                                 >  4:900721/1‑148 (MQ=255)
                 gcagcGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACg                                                                                                                            >  2:450436/1‑149 (MQ=255)
                       aCAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCa                                                                                                                       >  3:223125/1‑148 (MQ=255)
                       aCAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAg                                                                                                                      <  2:431654/149‑1 (MQ=255)
                            cTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGc                                                                                                                                            >  2:506859/1‑122 (MQ=255)
                            cTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGc                                                                                                                                            <  1:506859/122‑1 (MQ=255)
                             tGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATc                                                                                                                <  1:107052/149‑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                                                                                                            <  4:151449/148‑1 (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)
                                     gAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGt                                                                                                        <  2:776391/149‑1 (MQ=255)
                                      aaCTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGt                                                                                                        <  4:130974/148‑1 (MQ=255)
                                          ggcggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTAt                                                                                                    >  3:581125/1‑148 (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)
                                               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‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCaa                                                                                      >  3:531565/1‑148 (MQ=255)
                                                          tttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAcc                                                                                    >  3:40430/1‑148 (MQ=255)
                                                             aTCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGaa                                                                                <  2:49610/149‑1 (MQ=255)
                                                              tCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGaa                                                                                >  4:123539/1‑148 (MQ=255)
                                                               cGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAAc                                                                               >  3:515333/1‑148 (MQ=255)
                                                               cGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAAc                                                                               >  3:412526/1‑148 (MQ=255)
                                                                  ccccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCa                                                                           <  1:170957/149‑1 (MQ=255)
                                                                   cccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCag                                                                          <  1:542862/149‑1 (MQ=255)
                                                                   cccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCa                                                                           >  1:611746/1‑148 (MQ=255)
                                                                   cccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCa                                                                           <  2:611746/148‑1 (MQ=255)
                                                                    ccGATAGAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCaga                                                                         >  2:217436/1‑149 (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:988717/1‑148 (MQ=255)
                                                                                 gatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAAt                                                             >  3:967220/1‑148 (MQ=255)
                                                                                   ttCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                          >  2:260173/1‑149 (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)
                                                                                          gACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTgg                                                   <  2:125287/149‑1 (MQ=255)
                                                                                            ccGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTgggg                                                 >  1:6866/1‑149 (MQ=255)
                                                                                                   caGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgtg                                           >  4:382097/1‑148 (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                                                                                                                 <  3:370990/71‑1 (MQ=255)
                                                                                                          gagaGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACAt                                                                                                                 >  4:370990/1‑71 (MQ=255)
                                                                                                           agagTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTc                                  <  1:969982/149‑1 (MQ=255)
                                                                                                              gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgt                                            <  4:706219/136‑1 (MQ=255)
                                                                                                              gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgt                                            >  3:706219/1‑136 (MQ=255)
                                                                                                                ttCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTc                              >  4:803932/1‑148 (MQ=255)
                                                                                                                ttCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTc                              >  3:670564/1‑148 (MQ=255)
                                                                                                                   ggTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgat                          <  1:299622/149‑1 (MQ=255)
                                                                                                                   ggTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgat                          <  1:716285/149‑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)
                                                                                                                                aagaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAg             >  2:91613/1‑149 (MQ=255)
                                                                                                                                 agaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAg             >  4:95010/1‑148 (MQ=255)
                                                                                                                                  gaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCAtt                                                       <  4:968556/105‑1 (MQ=255)
                                                                                                                                  gaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCAtt                                                       >  3:968556/1‑105 (MQ=255)
                                                                                                                                  gaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTa           >  2:530094/1‑149 (MQ=255)
                                                                                                                                   aaaTCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTa           <  4:531565/148‑1 (MQ=255)
                                                                                                                                     aTCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAttt        >  2:23776/1‑149 (MQ=255)
                                                                                                                                      tCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgat                          <  2:313892/130‑1 (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)
                                                                                                                                           cACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                         <  4:779161/46‑1 (MQ=255)
                                                                                                                                           cACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                         >  3:779161/1‑46 (MQ=255)
                                                                                                                                            aCGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGtt  <  4:515333/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
GCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTT  >  NC_000913/3970672‑3970959

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