Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 704,236 T→C L357S (TTG→TCG)  nagE → fused N‑acetyl glucosamine specific PTS enzyme: IIC, IIB, and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913704,2360TC99.8% 2513.1 / NA 623L357S (TTG→TCG) nagEfused N‑acetyl glucosamine specific PTS enzyme: IIC, IIB, and IIA components
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (345/276);  total (346/277)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 651 total aligned reads displayed.

CGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTTGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTG  >  NC_000913/704138‑704333
                                                                                                  |                                                                                                 
   gTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTa                                                                                              >  2:6702117/1‑101 (MQ=255)
    ttGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTAt                                                                                             >  1:8363666/1‑101 (MQ=255)
      gATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATcc                                                                                           <  2:5874940/101‑1 (MQ=255)
        tGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGc                                                                                         <  1:14044202/101‑1 (MQ=255)
          tataACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCAt                                                                                       >  2:2092866/1‑101 (MQ=255)
          tataACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTGCAGTTCGGTTATCCgnnn                                                                                       >  2:13038416/1‑98 (MQ=255)
           ataACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGcnnn                                                                                      >  2:11897803/1‑98 (MQ=255)
            taACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGt                                                                                     <  1:4822221/101‑1 (MQ=255)
               nnngccggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGttta                                                                                  <  2:12165355/98‑3 (MQ=255)
               ccTGCCGCCCGCCACGCAGCACGTCTGGATGCTGCTGGTCATGGGCATTAACTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCa                                                                                  <  2:4527543/101‑1 (MQ=255)
                cTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCaa                                                                                 <  1:9809791/101‑1 (MQ=255)
                   ccggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCaannn                                                                              >  2:11002002/1‑98 (MQ=255)
                   ccggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCt                                                                              <  1:11167546/101‑1 (MQ=255)
                   ccggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCt                                                                              <  1:2863921/101‑1 (MQ=255)
                   ccggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCCTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCt                                                                              <  1:12571937/101‑1 (MQ=255)
                   ccggccgCCAGCCAGAACGTCTGGATGCTGCTGGGGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGCGTTCAGTTCGGTTATCCGCATGTTCAACCt                                                                              >  1:14375120/1‑101 (MQ=255)
                    cggccggcAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTg                                                                             >  1:14087928/1‑101 (MQ=255)
                    cggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTg                                                                             <  2:802132/101‑1 (MQ=255)
                     ggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGa                                                                            >  1:2324048/1‑101 (MQ=255)
                     ggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGa                                                                            <  1:799271/101‑1 (MQ=255)
                      gccgccAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaa                                                                           <  1:2769606/101‑1 (MQ=255)
                      gccgccAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaa                                                                           <  2:3143243/101‑1 (MQ=255)
                       ccgccAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTgnnn                                                                          >  2:7561279/1‑98 (MQ=255)
                        cgccAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATTTTCAACCTGaaaa                                                                         >  1:6353243/1‑101 (MQ=255)
                        cgccAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaaaa                                                                         <  1:7579838/101‑1 (MQ=255)
                         gccagccaGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAAc                                                                        <  1:5376321/101‑1 (MQ=255)
                         gccagccaGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAAc                                                                        >  2:5150672/1‑101 (MQ=255)
                          ccagccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACg                                                                       <  1:1207882/101‑1 (MQ=255)
                           cagccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaaaannn                                                                      >  2:9155446/1‑98 (MQ=255)
                           cagccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGc                                                                      >  2:2925641/1‑101 (MQ=255)
                              ccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGcnnn                                                                   >  2:7308081/1‑98 (MQ=255)
                              ccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCgg                                                                   <  2:3911725/101‑1 (MQ=255)
                                agAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGt                                                                 >  1:5092696/1‑101 (MQ=255)
                                 gAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTc                                                                <  2:4944652/101‑1 (MQ=255)
                                   nnntCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAGAACGCCGGGTCga                                                              <  2:14265510/98‑2 (MQ=255)
                                   aCGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGt                                                              >  1:4167716/1‑101 (MQ=255)
                                   aCGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGt                                                              >  1:11929742/1‑101 (MQ=255)
                                      tCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGtnnn                                                           >  2:12760972/1‑98 (MQ=255)
                                        tGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGa                                                         >  1:12378461/1‑101 (MQ=255)
                                        tGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGa                                                         >  1:10111631/1‑101 (MQ=255)
                                         ggATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGaannn                                                        >  2:11503454/1‑98 (MQ=255)
                                           aTGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATaa                                                      >  1:9527384/1‑101 (MQ=255)
                                            tgctgcTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGGGAAGAtnnn                                                     >  2:14579085/1‑98 (MQ=255)
                                            tgatgctgGTGATGGGCGTTATCTTCTTCGCTATCTCCTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATaaa                                                     <  1:8619376/98‑1 (MQ=255)
                                                nnnggTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaagaag                                                 <  2:9431902/98‑1 (MQ=255)
                                                gctgGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaagaag                                                 >  2:5298414/1‑101 (MQ=255)
                                                gctgGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaagaag                                                 >  1:8206944/1‑101 (MQ=255)
                                                 ctgGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaagaata                                                >  2:1615713/1‑99 (MQ=255)
                                                 ctgGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAagaaga                                                <  1:11663085/101‑1 (MQ=255)
                                                   nnngATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACg                                              <  2:8278960/98‑1 (MQ=255)
                                                     tGATGGGCGTTATCTTCTTCGCTATCTACGTCGTGGTGTTCAGTTCGGATATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACgag                                            >  1:4050643/1‑101 (MQ=255)
                                                      gATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACgaga                                           >  2:1014137/1‑101 (MQ=255)
                                                       aTGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAagaagaa                                               <  2:11688910/96‑2 (MQ=255)
                                                        tGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATc                                         <  2:629414/101‑1 (MQ=255)
                                                        nnngCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATc                                         <  2:8535719/98‑1 (MQ=255)
                                                          nnngTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGt                                       <  2:13796407/98‑1 (MQ=255)
                                                          ggCGTTATCTTCGTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGt                                       <  1:5854029/101‑1 (MQ=255)
                                                            ggTTATCTTCTTCGCTATCTACTTCGTGGTGTTTAGTTCGGTTATCCACATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTa                                     <  1:1894093/100‑1 (MQ=255)
                                                            cGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGAtggnnn                                     >  2:8909551/1‑96 (MQ=255)
                                                             gTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTAc                                    >  1:5538494/1‑101 (MQ=255)
                                                             gTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTAc                                    >  1:4374145/1‑101 (MQ=255)
                                                                aTCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCTGGTCGTGAAGATAAAGAAGACGAGATCGTTACTga                                 >  1:14530229/1‑101 (MQ=255)
                                                                 tcttctTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAGCGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACtnnn                                >  2:10546965/1‑98 (MQ=255)
                                                                    tcttcGCTATCTCCTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaagaa                             <  1:8932128/101‑1 (MQ=255)
                                                                    nnntcGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaagaa                             <  2:10797062/98‑1 (MQ=255)
                                                                      ttcGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGc                           <  2:3638074/101‑1 (MQ=255)
                                                                      ttcGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGc                           <  1:9251561/101‑1 (MQ=255)
                                                                      nnngCTATCTACTTCGTTGTGTTCAGTTCGGTTATCCGCATGTTCACCCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGc                           <  2:14755298/98‑1 (MQ=255)
                                                                       tcGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGcc                          <  1:11002002/101‑1 (MQ=255)
                                                                        cGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCa                         <  1:1357088/101‑1 (MQ=255)
                                                                         gCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCaa                        >  2:6142725/1‑101 (MQ=255)
                                                                         gCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCaa                        <  1:697305/101‑1 (MQ=255)
                                                                           tATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACa                      >  1:6016104/1‑101 (MQ=255)
                                                                            aTCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAg                     <  2:6948499/101‑1 (MQ=255)
                                                                            aTCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAg                     <  2:1823487/101‑1 (MQ=255)
                                                                             tCTACTTCGTGGTGTTCAGTACGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGc                    <  1:11422758/101‑1 (MQ=255)
                                                                              cTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCa                   >  1:7836316/1‑101 (MQ=255)
                                                                              cTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCa                   >  2:4679643/1‑101 (MQ=255)
                                                                               tACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAgnnn                  >  2:9826748/1‑98 (MQ=255)
                                                                               tACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCaa                  >  1:2100399/1‑101 (MQ=255)
                                                                                nnntCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCatc                 <  2:11421081/98‑3 (MQ=255)
                                                                                   tcttggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACt              <  1:13058515/98‑1 (MQ=255)
                                                                                   tCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACt              <  1:7344006/101‑1 (MQ=255)
                                                                                     gtggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAGGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTga            >  1:13788973/1‑101 (MQ=255)
                                                                                     gtggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAacacnnn            >  2:13548707/1‑98 (MQ=255)
                                                                                     gtggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTga            <  2:6518311/101‑1 (MQ=255)
                                                                                      tggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTAATGAAGAAGCCAACAGCAACACTgaa           <  2:2274306/101‑1 (MQ=255)
                                                                                      tggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGACAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTgaa           >  1:12853242/1‑101 (MQ=255)
                                                                                       nnngtTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTgaag          <  2:12378461/98‑1 (MQ=255)
                                                                                       ggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTgaag          <  1:3692710/101‑1 (MQ=255)
                                                                                         nnntCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTgaagaa        <  2:7836316/98‑1 (MQ=255)
                                                                                          gtTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAACAAGCCAACAGCAACACTGaagaag       >  2:7140089/1‑101 (MQ=255)
                                                                                           ttCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAgg      >  1:4622096/1‑101 (MQ=255)
                                                                                           ttCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGACAACAGCAACACTGAAGAAgg      >  2:2720064/1‑101 (MQ=255)
                                                                                            tCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGt     <  1:9928041/101‑1 (MQ=255)
                                                                                              aGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCt   >  2:4100867/1‑101 (MQ=255)
                                                                                              aGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCt   >  1:10572231/1‑101 (MQ=255)
                                                                                               gTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGCCGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTg  <  1:33025/101‑1 (MQ=255)
                                                                                               gTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTCCTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTg  <  1:6630202/101‑1 (MQ=255)
                                                                                               gTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTg  >  1:1577219/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
CGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTTGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAGAAGGTCTG  >  NC_000913/704138‑704333

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