Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,970,817:1 +G 8.0% 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 8.0% 416.2 / 33.5 100coding (685/1131 nt)wecBUDP‑N‑acetyl glucosamine‑2‑epimerase
Reads supporting (aligned to +/- strand):  ref base . (46/46);  new base G (4/4);  total (50/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

Only 100 of 104 total aligned reads displayed.

GGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGT  >  NC_000913/3970669‑3970962
                                                                                                                                                     |                                                                                                                                                 
ggTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGACACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACgcgc                                                                                                                                                    >  1:284456/1‑149 (MQ=255)
 gTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑t                                                                                                                                                  >  2:90011/1‑149 (MQ=255)
 gTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑t                                                                                                                                                  <  1:215476/149‑1 (MQ=255)
   gCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑Tgg                                                                                                                                                <  2:135058/149‑1 (MQ=255)
    cGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGc                                                                                                                                               <  4:414445/149‑1 (MQ=255)
     gTGACCAGGTGATGGGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCa                                                                                                                                              >  3:45884/1‑149 (MQ=255)
     gTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGATCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCa                                                                                                                                              >  4:468995/1‑149 (MQ=255)
        aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGAc                                                                                                                                           <  2:421718/149‑1 (MQ=255)
        aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGat                                                                                                                                           <  4:249418/149‑2 (MQ=255)
        aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGat                                                                                                                                           <  3:438728/149‑2 (MQ=255)
        aCCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCA‑GCGCGTGGCAGat                                                                                                                                           <  2:4171/149‑2 (MQ=255)
           aGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATc                                                                                                                                        <  1:147591/149‑1 (MQ=255)
              tgatgaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGcc                                                                                                                                     <  1:501967/149‑1 (MQ=255)
                  gaGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGccacca                                                                                                                                 >  3:354311/1‑149 (MQ=255)
                     cagcGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGc                                                                                                                              >  1:287853/1‑149 (MQ=255)
                         gACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCa                                                                                                                          >  2:342667/1‑149 (MQ=255)
                          aCAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAg                                                                                                                         >  4:329153/1‑149 (MQ=255)
                          aCAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAg                                                                                                                         <  3:450388/149‑1 (MQ=255)
                            aaGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGa                                                                                                                       <  4:485597/149‑1 (MQ=255)
                               cTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTTTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACAt                                                                                                                    <  4:545239/149‑1 (MQ=255)
                                tGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATc                                                                                                                   <  3:70144/149‑1 (MQ=255)
                                 gCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGT‑TGGCAGACATCGCCACCACGCACCAGGACATcc                                                                                                                  <  3:391673/149‑1 (MQ=255)
                                   gTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAg                                                                                                                <  4:229750/149‑1 (MQ=255)
                                      cAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGAtt                                                                                                             <  1:363367/149‑1 (MQ=255)
                                       aGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGGCCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTg                                                                                                            >  1:408222/1‑149 (MQ=255)
                                        gAACTGGCGGCAAATTACCCGTTTATCGATCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGt                                                                                                           <  3:468995/149‑1 (MQ=255)
                                        gAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGt                                                                                                           <  3:374473/149‑1 (MQ=255)
                                         aaCTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTc                                                                                                          <  3:444565/149‑1 (MQ=255)
                                         aaCTGGCGGCAAATTACCCGTTTATCGAACCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCGTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTc                                                                                                          >  4:380682/1‑149 (MQ=255)
                                          aCTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCt                                                                                                         <  1:68093/149‑1 (MQ=255)
                                          aCTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCt                                                                                                         <  1:432553/149‑1 (MQ=255)
                                            tGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTAt                                                                                                       <  2:319592/149‑1 (MQ=255)
                                                ggcAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCgg                                                                                                   <  2:287853/149‑1 (MQ=255)
                                                      ttACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCAtc                                                                                             >  3:224339/1‑149 (MQ=255)
                                                      ttACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCAtc                                                                                             >  3:54963/1‑149 (MQ=255)
                                                         cccGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCa                                                                                          >  4:487568/1‑149 (MQ=255)
                                                           cGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAc                                                                                        <  3:44244/149‑1 (MQ=255)
                                                            gTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAcc                                                                                       <  2:146875/149‑1 (MQ=255)
                                                             tttATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAAccc                                                                                      >  3:144904/1‑149 (MQ=255)
                                                                    aCCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTc                                                                               >  1:329292/1‑149 (MQ=255)
                                                                       ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCaga                                                                            <  1:90011/149‑1 (MQ=255)
                                                                       ccGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTGTCGGTCGTGGCGTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGCTTGTCTATCCGGTGCATCTCAACCCGAACGTCaga                                                                            >  2:390310/1‑149 (MQ=255)
                                                                        cGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCagag                                                                           >  3:82695/1‑149 (MQ=255)
                                                                          aTAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGaa                                                                         <  1:571782/149‑1 (MQ=255)
                                                                            aaaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAAcc                                                                       >  2:361862/1‑149 (MQ=255)
                                                                             aaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCg                                                                      <  3:600566/149‑1 (MQ=255)
                                                                             aaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCg                                                                      <  1:15285/149‑1 (MQ=255)
                                                                             aaaaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATACGGTGCATCTCAACCCGAACGTCAGAGAACCg                                                                      <  2:488268/149‑1 (MQ=255)
                                                                               aaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGTCTTTGAAGAAATCTGCCACGCGC‑TGACAGACATCGCCACCACTCACCAGGACATTCAGATTTTCTATCCGGTGCATCTCAACCTGAACGTCAGAGAACCGGt                                                                    <  2:284456/149‑1 (MQ=255)
                                                                               aaGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGt                                                                    <  1:342667/149‑1 (MQ=255)
                                                                                      ttCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                             >  1:195343/1‑149 (MQ=255)
                                                                                      ttCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                             >  4:427237/1‑149 (MQ=255)
                                                                                      ttCTGGTGACCGGTCACAGGCGTGAGAGGTTCCGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCAAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGc                                                             >  4:380562/1‑149 (MQ=255)
                                                                                       tCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCa                                                            <  3:216743/149‑1 (MQ=255)
                                                                                             gACCGGTCACAGGCGTGACAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGACTATCCGGAGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCCCATTCTga                                                      >  2:299437/1‑148 (MQ=255)
                                                                                               ccGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCACATTCTgggg                                                    >  2:488525/1‑149 (MQ=255)
                                                                                                    cacaGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgt                                               >  2:98566/1‑149 (MQ=255)
                                                                                                    cacaGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCAtgt                                               <  2:420216/149‑1 (MQ=255)
                                                                                                     acaGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTgg                                                      <  3:451367/141‑1 (MQ=255)
                                                                                                     acaGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTgg                                                      >  4:451367/1‑141 (MQ=255)
                                                                                                       aGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaa                                            >  3:459138/1‑149 (MQ=255)
                                                                                                       aGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaa                                            >  4:559143/1‑149 (MQ=255)
                                                                                                        ggCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaa                                           >  2:112819/1‑149 (MQ=255)
                                                                                                        ggCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaa                                           >  2:308550/1‑149 (MQ=255)
                                                                                                         gCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGaaaa                                          >  2:483984/1‑149 (MQ=255)
                                                                                                           gTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAAt                                        >  1:123340/1‑149 (MQ=255)
                                                                                                            tGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATg                                       >  3:430814/1‑149 (MQ=255)
                                                                                                            tGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATg                                       >  1:68148/1‑149 (MQ=255)
                                                                                                              agggTTTCGCTCGTGGCTTAGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTc                                     >  2:473857/4‑149 (MQ=255)
                                                                                                                agTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCAt                                   <  4:224339/149‑1 (MQ=255)
                                                                                                                 gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCAtt                                  >  2:472755/1‑149 (MQ=255)
                                                                                                                 gTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCAtt                                  <  1:81466/149‑1 (MQ=255)
                                                                                                                  tttCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCAtt                                  >  3:294335/1‑148 (MQ=255)
                                                                                                                  tttCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCAtt                                  <  4:294335/148‑1 (MQ=255)
                                                                                                                  tttCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTc                                 >  4:452078/1‑149 (MQ=255)
                                                                                                                       gTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatc                            <  1:98566/149‑1 (MQ=255)
                                                                                                                        tCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCa                                    >  1:193258/1‑140 (MQ=255)
                                                                                                                        tCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCa                                    <  2:193258/140‑1 (MQ=255)
                                                                                                                           tGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTgatcgatc                        >  1:209475/1‑149 (MQ=255)
                                                                                                                             gCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATc                                                               >  2:149503/1‑108 (MQ=255)
                                                                                                                             gCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATc                                                               <  1:149503/108‑1 (MQ=255)
                                                                                                                               tttGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACGTCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCa                    <  1:37488/149‑1 (MQ=255)
                                                                                                                                ttGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAg                   >  1:360892/1‑149 (MQ=255)
                                                                                                                                  gaagaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa                 >  3:245293/1‑149 (MQ=255)
                                                                                                                                  gaagaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa                 <  4:354311/149‑1 (MQ=255)
                                                                                                                                  gaagaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGa                 >  4:462688/1‑149 (MQ=255)
                                                                                                                                    agaaATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTATATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGt               >  1:283695/1‑149 (MQ=255)
                                                                                                                                      aaaTCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAt             >  1:291497/1‑149 (MQ=255)
                                                                                                                                      aaaTCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAt             <  1:390310/149‑1 (MQ=255)
                                                                                                                                       tgtCTGCCA‑GCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAtt            >  4:264760/3‑149 (MQ=255)
                                                                                                                                       tgtCTGCCA‑GCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTAtt            >  3:429546/3‑149 (MQ=255)
                                                                                                                                        gtCTGCCA‑GCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTa              >  4:171986/2‑146 (MQ=255)
                                                                                                                                        gtCTGCCA‑GCGCGTGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTa              <  3:171986/145‑1 (MQ=255)
                                                                                                                                         tCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTa          >  2:289407/1‑149 (MQ=255)
                                                                                                                                           tGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTAcc        <  3:487568/149‑1 (MQ=255)
                                                                                                                                             ccACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGt      <  1:158372/149‑1 (MQ=255)
                                                                                                                                              cACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCa                                    <  4:144238/118‑1 (MQ=255)
                                                                                                                                              cACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCa                                    >  3:144238/1‑118 (MQ=255)
                                                                                                                                               aCGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGttt    <  3:216980/149‑1 (MQ=255)
                                                                                                                                                 gcgc‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGt  <  3:462688/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GGTGCGTGACCAGGTGATGAGCAGCGACAAGCTGCGTTCAGAACTGGCGGCAAATTACCCGTTTATCGACCCCGATAAAAAGATGATTCTGGTGACCGGTCACAGGCGTGAGAGTTTCGGTCGTGGCTTTGAAGAAATCTGCCACGCGC‑TGGCAGACATCGCCACCACGCACCAGGACATCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGT  >  NC_000913/3970669‑3970962

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