Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,970,813 Δ1 bp 6.1% 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. 6.1% 404.0 / 18.6 115coding (681/1131 nt)wecBUDP‑N‑acetyl glucosamine‑2‑epimerase
Reads supporting (aligned to +/- strand):  ref base C (55/52);  new base . (4/3);  total (59/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

GTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTT  >  NC_000913/3970670‑3970960
                                                                                                                                               |                                                                                                                                                    
gTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTTGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACgcgc                                                                                                                                                  <  1:743213/148‑1 (MQ=255)
      gACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGa                                                                                                                                          <  3:710306/149‑1 (MQ=255)
        ccAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGAc                                                                                                                                         >  2:155149/1‑148 (MQ=255)
             tgatgaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGc                                                                                                                                    >  2:591873/1‑148 (MQ=255)
                 gaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGccacca                                                                                                                               <  3:128281/149‑1 (MQ=255)
                   gcagcGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACg                                                                                                                             >  3:806238/1‑149 (MQ=255)
                     agcGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGc                                                                                                                            >  2:128284/1‑148 (MQ=255)
                       cGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATTATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCAcc                                                                                                                         <  4:451610/149‑1 (MQ=255)
                        gACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCAcc                                                                                                                         <  1:544710/148‑1 (MQ=255)
                        gACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCAcc                                                                                                                         <  1:155149/148‑1 (MQ=255)
                         aCAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCa                                                                                                                        >  1:269289/1‑148 (MQ=255)
                         aCAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAg                                                                                                                       >  4:229372/1‑149 (MQ=255)
                                gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATcc                                                                                                                >  4:792495/1‑149 (MQ=255)
                                gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATc                                                                                                                 >  2:310749/1‑148 (MQ=255)
                                gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCATGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATc                                                                                                                 <  2:670773/148‑1 (MQ=255)
                                  gTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCa                                                                                                               <  1:764651/148‑1 (MQ=255)
                                   ttCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGa                                                                                                             >  4:280141/1‑149 (MQ=255)
                                     cAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGAtt                                                                                                           <  4:450692/149‑1 (MQ=255)
                                       gAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                          <  1:128284/148‑1 (MQ=255)
                                       gAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                          <  2:239699/148‑1 (MQ=255)
                                       gAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                          <  1:222098/148‑1 (MQ=255)
                                       gAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGt                                                                                                         >  4:627098/1‑149 (MQ=255)
                                            ggcggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTAt                                                                                                     <  1:591873/148‑1 (MQ=255)
                                                  aaaTTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGc                                                                                              >  4:72344/1‑149 (MQ=255)
                                                    aTTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCa                                                                                             <  1:274758/148‑1 (MQ=255)
                                                      tACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCAtc                                                                                           <  2:257718/148‑1 (MQ=255)
                                                       aCCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCAtctc                                                                                         >  4:441730/1‑149 (MQ=255)
                                                          cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCaa                                                                                       >  2:798104/1‑148 (MQ=255)
                                                          cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCaa                                                                                       >  1:607626/1‑148 (MQ=255)
                                                          cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAc                                                                                      <  3:751700/149‑1 (MQ=255)
                                                          cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAc                                                                                      <  4:853321/149‑1 (MQ=255)
                                                            tttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAccc                                                                                    >  4:742473/1‑149 (MQ=255)
                                                                tCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGaa                                                                                 <  2:250239/148‑1 (MQ=255)
                                                                   aCCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGt                                                                              <  2:313925/148‑1 (MQ=255)
                                                                    ccccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGat                                                                                                                                         <  4:296941/88‑2 (MQ=255)
                                                                    ccccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCAGCGCGTGGCAGat                                                                                                                                         >  3:296941/1‑87 (MQ=255)
                                                                     cccGATAAAAAGATGATTCTGGTGACCGCTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCag                                                                           >  4:452833/1‑149 (MQ=255)
                                                                       cGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCaga                                                                          >  1:736539/1‑148 (MQ=255)
                                                                        gATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGccacc                                                                                                                                <  2:528852/93‑1 (MQ=255)
                                                                        gATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGccacc                                                                                                                                >  1:528852/1‑93 (MQ=255)
                                                                        gATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCagag                                                                         >  1:409355/1‑148 (MQ=255)
                                                                         aTAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCGCGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGaa                                                                       <  4:549718/149‑1 (MQ=255)
                                                                           aaaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAAc                                                                      <  1:760641/148‑1 (MQ=255)
                                                                                gatgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCa                                                                >  3:49200/1‑149 (MQ=255)
                                                                                gatgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCa                                                                >  3:842565/1‑149 (MQ=255)
                                                                                 atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCaa                                                               <  3:535636/149‑1 (MQ=255)
                                                                                 atgatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCa                                                                <  1:684123/148‑1 (MQ=255)
                                                                                   gatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAAt                                                              >  1:418799/1‑148 (MQ=255)
                                                                                   gatTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATc                                                             <  4:842565/149‑1 (MQ=255)
                                                                                     ttCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                           >  3:384658/1‑149 (MQ=255)
                                                                                      tCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                           >  2:797844/1‑148 (MQ=255)
                                                                                          gTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCt                                                      <  4:806238/149‑1 (MQ=255)
                                                                                             aCCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTggg                                                   <  3:229372/149‑1 (MQ=255)
                                                                                             aCCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTggg                                                   <  3:785330/149‑1 (MQ=255)
                                                                                             aCCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTgg                                                    <  1:215040/148‑1 (MQ=255)
                                                                                               cGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGc                                                 >  3:635430/1‑149 (MQ=255)
                                                                                                ggTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACg                                                                                                                             >  4:467389/1‑72 (MQ=255)
                                                                                                ggTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACg                                                                                                                             <  3:467389/72‑1 (MQ=255)
                                                                                                ggTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGc                                                 >  1:116961/1‑148 (MQ=255)
                                                                                                 gTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAt                                               >  4:322345/1‑149 (MQ=255)
                                                                                                 gTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAt                                               <  3:742473/149‑1 (MQ=255)
                                                                                                  tCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAt                                               >  2:325563/1‑148 (MQ=255)
                                                                                                   cacaGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgt                                             <  4:635430/149‑1 (MQ=255)
                                                                                                   cacaGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgt                                             <  3:792495/149‑1 (MQ=255)
                                                                                                    acaGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGa                                                                                                             <  4:115863/84‑1 (MQ=255)
                                                                                                    acaGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGa                                                                                                             >  3:115863/1‑84 (MQ=255)
                                                                                                     caGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTg                                                     <  1:740765/139‑1 (MQ=255)
                                                                                                     caGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgtg                                            <  1:436003/148‑1 (MQ=255)
                                                                                                      aGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaa                                          >  3:238181/1‑149 (MQ=255)
                                                                                                         cGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaaa                                       <  4:722625/149‑1 (MQ=255)
                                                                                                         cGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaa                                        <  1:797844/148‑1 (MQ=255)
                                                                                                          gTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaaa                                       <  1:722538/148‑1 (MQ=255)
                                                                                                          gTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaaa                                       <  1:245756/148‑1 (MQ=255)
                                                                                                          gTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAAt                                      >  3:371564/1‑149 (MQ=255)
                                                                                                           tGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAAt                                      >  2:772870/1‑148 (MQ=255)
                                                                                                              gagTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTc                                   >  2:739607/1‑148 (MQ=255)
                                                                                                                gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCAtt                                <  3:74567/149‑1 (MQ=255)
                                                                                                                gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCAt                                 >  2:151033/1‑148 (MQ=255)
                                                                                                                       tCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatc                          >  2:616814/1‑148 (MQ=255)
                                                                                                                        cGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatcg                         >  1:864942/1‑148 (MQ=255)
                                                                                                                            gCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATccc                    >  4:125997/1‑149 (MQ=255)
                                                                                                                             cTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATccc                    >  2:372162/1‑148 (MQ=255)
                                                                                                                               ttGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCa                  >  1:586833/1‑148 (MQ=255)
                                                                                                                                 gaagaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa               >  3:336395/1‑149 (MQ=255)
                                                                                                                                   agaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGt             >  4:537655/1‑149 (MQ=255)
                                                                                                                                   agaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGCCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa               >  3:542460/1‑147 (MQ=255)
                                                                                                                                   agaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGCCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa               <  4:542460/147‑1 (MQ=255)
                                                                                                                                    gatgtCTGCCAGCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTa            >  3:487703/5‑149 (MQ=255)
                                                                                                                                    gaaATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTa            >  3:522041/1‑149 (MQ=255)
                                                                                                                                     atgtCTGCCAGCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTa            >  1:149438/4‑148 (MQ=255)
                                                                                                                                      aaTCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAtt          <  4:225998/149‑1 (MQ=255)
                                                                                                                                      aaTCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAtt          <  4:238181/149‑1 (MQ=255)
                                                                                                                                      aaTCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAt           >  2:654764/1‑148 (MQ=255)
                                                                                                                                        tCTGCCAGCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa               >  4:393088/1‑142 (MQ=255)
                                                                                                                                        tCTGCCAGCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa               <  3:393088/142‑1 (MQ=255)
                                                                                                                                         cTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTa        >  1:431649/1‑148 (MQ=255)
                                                                                                                                         cTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTAc       >  3:476727/1‑149 (MQ=255)
                                                                                                                                            ccACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGt    <  4:336395/149‑1 (MQ=255)
                                                                                                                                              aCGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGtt   >  2:75233/1‑148 (MQ=255)
                                                                                                                                               cgcgCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGttt  >  1:881316/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGCTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTT  >  NC_000913/3970670‑3970960

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