Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,970,817:1 +G 6.1% 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 6.1% 500.6 / 26.7 115coding (685/1131 nt)wecBUDP‑N‑acetyl glucosamine‑2‑epimerase
Reads supporting (aligned to +/- strand):  ref base . (56/52);  new base G (4/3);  total (60/55)
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.90e-01

Only 100 of 119 total aligned reads displayed.

GTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTC  >  NC_000913/3970670‑3970963
                                                                                                                                                    |                                                                                                                                                  
      gACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGa                                                                                                                                             <  3:710306/149‑1 (MQ=255)
        ccAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGAc                                                                                                                                            >  2:155149/1‑148 (MQ=255)
          aGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACAt                                                                                                                                          >  1:250239/1‑148 (MQ=255)
             tgatgaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGc                                                                                                                                       >  2:591873/1‑148 (MQ=255)
                tgaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGccac                                                                                                                                    >  2:245756/1‑148 (MQ=255)
                 gaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGccacca                                                                                                                                  <  3:128281/149‑1 (MQ=255)
                   gcagcGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACg                                                                                                                                >  3:806238/1‑149 (MQ=255)
                     agcGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGc                                                                                                                               >  2:128284/1‑148 (MQ=255)
                        gACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCAcc                                                                                                                            <  1:155149/148‑1 (MQ=255)
                        gACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCAcc                                                                                                                            <  1:544710/148‑1 (MQ=255)
                         aCAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCa                                                                                                                           >  1:269289/1‑148 (MQ=255)
                         aCAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCa                                                                                                                           <  2:495739/148‑1 (MQ=255)
                         aCAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAg                                                                                                                          >  4:229372/1‑149 (MQ=255)
                               tGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATc                                                                                                                    <  4:145821/149‑1 (MQ=255)
                                gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATcc                                                                                                                   >  4:792495/1‑149 (MQ=255)
                                gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATc                                                                                                                    >  2:310749/1‑148 (MQ=255)
                                gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATc                                                                                                                    <  1:722047/148‑1 (MQ=255)
                                  gTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCa                                                                                                                  <  1:764651/148‑1 (MQ=255)
                                   ttCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGa                                                                                                                >  4:280141/1‑149 (MQ=255)
                                     cAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGAtt                                                                                                              <  4:450692/149‑1 (MQ=255)
                                      aGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                             <  4:376746/149‑1 (MQ=255)
                                       gAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                             <  2:239699/148‑1 (MQ=255)
                                       gAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                             <  1:222098/148‑1 (MQ=255)
                                       gAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                             <  1:128284/148‑1 (MQ=255)
                                       gAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGt                                                                                                            >  4:627098/1‑149 (MQ=255)
                                       gAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGt                                                                                                            <  4:722828/149‑1 (MQ=255)
                                               ggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCg                                                                                                     <  1:452509/148‑1 (MQ=255)
                                                    aTTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCa                                                                                                <  1:274758/148‑1 (MQ=255)
                                                      tACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCAtc                                                                                              <  2:257718/148‑1 (MQ=255)
                                                       aCCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCAtctc                                                                                            >  4:441730/1‑149 (MQ=255)
                                                          cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCaa                                                                                          >  1:607626/1‑148 (MQ=255)
                                                          cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAc                                                                                         <  3:751700/149‑1 (MQ=255)
                                                          cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAc                                                                                         <  4:853321/149‑1 (MQ=255)
                                                            tttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAccc                                                                                       >  4:742473/1‑149 (MQ=255)
                                                                   aCCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGt                                                                                 <  2:313925/148‑1 (MQ=255)
                                                                     cccGATAAAAAGATGATTCTGGTGACCGCTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCag                                                                              >  4:452833/1‑149 (MQ=255)
                                                                       cGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCaga                                                                             >  1:736539/1‑148 (MQ=255)
                                                                        gATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGccacc                                                                                                                                   <  2:528852/93‑1 (MQ=255)
                                                                        gATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGccacc                                                                                                                                   >  1:528852/1‑93 (MQ=255)
                                                                        gATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCagag                                                                            >  1:409355/1‑148 (MQ=255)
                                                                         aTAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCGCGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGaa                                                                          <  4:549718/149‑1 (MQ=255)
                                                                           aaaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAAc                                                                         <  1:760641/148‑1 (MQ=255)
                                                                                gatgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCa                                                                   >  3:49200/1‑149 (MQ=255)
                                                                                gatgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCa                                                                   >  3:842565/1‑149 (MQ=255)
                                                                                 atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCaa                                                                  <  3:535636/149‑1 (MQ=255)
                                                                                 atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCa                                                                   <  1:684123/148‑1 (MQ=255)
                                                                                   gatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAAt                                                                 >  1:418799/1‑148 (MQ=255)
                                                                                   gatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATc                                                                <  4:842565/149‑1 (MQ=255)
                                                                                     ttCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                              >  3:384658/1‑149 (MQ=255)
                                                                                      tCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                              >  2:797844/1‑148 (MQ=255)
                                                                                      tCTGGTGACCGCTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCa                                                             <  3:452833/149‑1 (MQ=255)
                                                                                          gTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCt                                                         <  4:806238/149‑1 (MQ=255)
                                                                                             aCCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTggg                                                      <  3:229372/149‑1 (MQ=255)
                                                                                             aCCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTggg                                                      <  3:785330/149‑1 (MQ=255)
                                                                                             aCCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTgg                                                       <  1:215040/148‑1 (MQ=255)
                                                                                               cGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGc                                                    >  3:635430/1‑149 (MQ=255)
                                                                                                ggTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACg                                                                                                                                >  4:467389/1‑72 (MQ=255)
                                                                                                ggTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGc                                                    >  1:116961/1‑148 (MQ=255)
                                                                                                 gTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAt                                                  >  4:322345/1‑149 (MQ=255)
                                                                                                 gTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAt                                                  <  3:742473/149‑1 (MQ=255)
                                                                                                  tCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAt                                                  >  2:325563/1‑148 (MQ=255)
                                                                                                   cacaGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgt                                                <  3:792495/149‑1 (MQ=255)
                                                                                                   cacaGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgt                                                <  4:635430/149‑1 (MQ=255)
                                                                                                    acaGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGa                                                                                                                <  4:115863/84‑1 (MQ=255)
                                                                                                    acaGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGa                                                                                                                >  3:115863/1‑84 (MQ=255)
                                                                                                     caGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTg                                                        >  2:740765/1‑139 (MQ=255)
                                                                                                     caGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTg                                                        <  1:740765/139‑1 (MQ=255)
                                                                                                     caGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgtg                                               <  1:436003/148‑1 (MQ=255)
                                                                                                      aGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaa                                             >  3:238181/1‑149 (MQ=255)
                                                                                                      aGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGa                                              >  1:784950/1‑148 (MQ=255)
                                                                                                         cGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaaa                                          <  4:722625/149‑1 (MQ=255)
                                                                                                         cGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaa                                           <  1:797844/148‑1 (MQ=255)
                                                                                                          gTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaaa                                          <  1:722538/148‑1 (MQ=255)
                                                                                                          gTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAAt                                         >  3:371564/1‑149 (MQ=255)
                                                                                                           tGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAAt                                         >  2:772870/1‑148 (MQ=255)
                                                                                                             agagTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTc                                      >  3:536665/1‑149 (MQ=255)
                                                                                                              gagTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTc                                      >  2:739607/1‑148 (MQ=255)
                                                                                                                gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCAtt                                   <  3:74567/149‑1 (MQ=255)
                                                                                                                gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCAt                                    >  2:151033/1‑148 (MQ=255)
                                                                                                                       tCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatc                             >  2:616814/1‑148 (MQ=255)
                                                                                                                        cGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatcg                            >  1:864942/1‑148 (MQ=255)
                                                                                                                               ttGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCa                     >  1:586833/1‑148 (MQ=255)
                                                                                                                                 gaagaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa                  >  3:336395/1‑149 (MQ=255)
                                                                                                                                 gaagaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa                  <  4:49200/149‑1 (MQ=255)
                                                                                                                                   agaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGt                >  4:537655/1‑149 (MQ=255)
                                                                                                                                    gatgtCTGCCA‑GCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTa               >  3:487703/5‑149 (MQ=255)
                                                                                                                                    gaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTa               >  3:522041/1‑149 (MQ=255)
                                                                                                                                     atgtCTGCCA‑GCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTa               >  1:149438/4‑148 (MQ=255)
                                                                                                                                      aaTCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAtt             <  4:238181/149‑1 (MQ=255)
                                                                                                                                      aaTCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAtt             <  4:225998/149‑1 (MQ=255)
                                                                                                                                      aaTCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAtt             >  3:225998/1‑149 (MQ=255)
                                                                                                                                      aaTCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAt              >  2:654764/1‑148 (MQ=255)
                                                                                                                                        tCTGCCA‑GCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa                  >  4:393088/1‑142 (MQ=255)
                                                                                                                                        tCTGCCA‑GCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa                  <  3:393088/142‑1 (MQ=255)
                                                                                                                                         cTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTAc          >  3:476727/1‑149 (MQ=255)
                                                                                                                                            ccACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGt       <  4:336395/149‑1 (MQ=255)
                                                                                                                                              aCGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGtt      >  2:75233/1‑148 (MQ=255)
                                                                                                                                                 cgc‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTGCCGTTTGTc  <  4:536665/149‑1 (MQ=255)
                                                                                                                                                  gc‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTc  <  2:269289/148‑1 (MQ=255)
                                                                                                                                                  gc‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTc  <  2:489973/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTC  >  NC_000913/3970670‑3970963

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