Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 704,130 A→G I322V (ATC→GTC)  nagE → PTS system, N‑acetylglucosamine‑specific enzyme IIABC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913704,1300AG89.0% 580.2 / 56.3 228I322V (ATC→GTC) nagEPTS system, N‑acetylglucosamine‑specific enzyme IIABC
Reads supporting (aligned to +/- strand):  ref base A (16/9);  new base G (101/102);  total (117/111)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 230 total aligned reads displayed.

TACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTC  >  NC_000913/704040‑704222
                                                                                           |                                                                                            
tACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTc                                                                                            <  1:1773611/93‑1 (MQ=255)
tACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTc                                                                                            <  1:5014078/93‑1 (MQ=255)
 aCCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCg                                                                                           <  1:121567/93‑1 (MQ=255)
  cctcctGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGa                                                                                          >  2:6938742/1‑93 (MQ=255)
   ctcctGCACGCACTTCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGAc                                                                                         >  2:5197307/1‑93 (MQ=255)
    tcctGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACt                                                                                        >  2:2880588/1‑93 (MQ=255)
    tcctGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACt                                                                                        >  2:5665664/1‑93 (MQ=255)
     cctGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTTTCTTTCTCTGCGGGGGCTGTCGCCTa                                                                                       <  1:6971079/93‑1 (MQ=255)
     cctGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTa                                                                                       <  1:5699992/93‑1 (MQ=255)
       tGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTAcg                                                                                     >  1:5099343/1‑93 (MQ=255)
       tGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTAcg                                                                                     <  2:3608066/93‑1 (MQ=255)
       tGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTAcg                                                                                     >  2:203387/1‑93 (MQ=255)
       tGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTAcg                                                                                     >  1:1541605/1‑93 (MQ=255)
          acgcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGt                                                                                  >  2:6697555/1‑93 (MQ=255)
          acgcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGt                                                                                  <  1:2288848/93‑1 (MQ=255)
            gcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTg                                                                                <  1:6301474/93‑1 (MQ=255)
             cacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGa                                                                               <  1:5890021/93‑1 (MQ=255)
             cacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGa                                                                               <  2:4280958/93‑1 (MQ=255)
             cacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGa                                                                               >  2:1776421/1‑93 (MQ=255)
             cacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGa                                                                               <  2:6701373/93‑1 (MQ=255)
               ctgctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATg                                                                             >  1:5334701/1‑93 (MQ=255)
                tgctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGt                                                                            >  1:4600184/1‑93 (MQ=255)
                 gctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGta                                                                           >  1:736935/1‑93 (MQ=255)
                 gctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGta                                                                           <  1:3078188/93‑1 (MQ=255)
                   tgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGtata                                                                         >  2:6084471/1‑93 (MQ=255)
                   tgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGtata                                                                         >  1:4504888/1‑93 (MQ=255)
                   tgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGtata                                                                         >  1:2718686/1‑93 (MQ=255)
                   tgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGtata                                                                         >  2:2804373/1‑93 (MQ=255)
                   tgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGtata                                                                         >  2:6776415/1‑93 (MQ=255)
                    gACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATaa                                                                        >  2:2740570/1‑93 (MQ=255)
                     aCCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAAc                                                                       <  1:2470771/93‑1 (MQ=255)
                     aCCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTAGACTACGCGTTGATGTATAAc                                                                       <  1:503464/93‑1 (MQ=255)
                      ccGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAAcc                                                                      >  1:959878/1‑93 (MQ=255)
                      ccGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAAcc                                                                      >  2:6240179/1‑93 (MQ=255)
                      ccGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAAcc                                                                      >  1:3457670/1‑93 (MQ=255)
                       cGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCt                                                                     <  1:6208258/93‑1 (MQ=255)
                        ggTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTg                                                                    <  1:2774668/93‑1 (MQ=255)
                        ggTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTg                                                                    >  2:145343/1‑93 (MQ=255)
                         gTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTgc                                                                   <  2:1099068/93‑1 (MQ=255)
                         gTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTgc                                                                   <  1:12786/93‑1 (MQ=255)
                         gTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGCGGGCTATCGACTACGCGTTGATGTATAACCTgc                                                                   >  2:1348653/1‑93 (MQ=255)
                          tATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTgcc                                                                  >  1:1027283/1‑93 (MQ=255)
                           aTCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTgccg                                                                 >  2:4757370/1‑93 (MQ=255)
                            tCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGGCTGTCGACTACGCGTTGATGTATAACCTgccg                                                                 >  1:2412850/1‑93 (MQ=255)
                              aGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTgccggcc                                                              >  1:4820919/1‑93 (MQ=255)
                                 cTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGgccgcc                                                           >  1:7026036/1‑93 (MQ=255)
                                    tttGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCgccagc                                                        <  2:5612004/93‑1 (MQ=255)
                                    tttGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCgccagc                                                        >  2:6529417/1‑93 (MQ=255)
                                     ttGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCgccagcc                                                       <  1:1903689/93‑1 (MQ=255)
                                         ggCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGaa                                                   <  1:938631/93‑1 (MQ=255)
                                         ggCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGaa                                                   <  2:3833406/93‑1 (MQ=255)
                                         ggCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGaa                                                   >  2:7024426/1‑93 (MQ=255)
                                          gCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAc                                                  >  2:2314124/1‑93 (MQ=255)
                                             aCGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTc                                               <  1:6089218/93‑1 (MQ=255)
                                                ctgctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTgg                                            <  2:4831193/93‑1 (MQ=255)
                                                 tgctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGa                                           <  2:6458645/93‑1 (MQ=255)
                                                 tgctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGa                                           <  2:4193058/93‑1 (MQ=255)
                                                 tgctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGa                                           <  2:6408339/93‑1 (MQ=255)
                                                  gctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAt                                          >  2:3233002/1‑93 (MQ=255)
                                                   ctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtg                                         >  2:4372141/1‑93 (MQ=255)
                                                     gggTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtgct                                       <  1:3042627/93‑1 (MQ=255)
                                                     gggTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtgct                                       >  2:5849317/1‑93 (MQ=255)
                                                      ggTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtgctg                                      >  2:5346784/1‑93 (MQ=255)
                                                       gTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtgctgc                                     <  1:1193076/93‑1 (MQ=255)
                                                       gTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtgctgc                                     >  2:7395211/1‑93 (MQ=255)
                                                        tATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATgctgct                                    <  1:2758110/93‑1 (MQ=255)
                                                         aTCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCAACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGctgctg                                   <  1:744241/93‑1 (MQ=255)
                                                          tCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTgg                                  >  1:4077509/1‑93 (MQ=255)
                                                          tCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTgg                                  >  1:4913758/1‑93 (MQ=255)
                                                          tCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTgg                                  >  1:438424/1‑93 (MQ=255)
                                                           ccACGCGGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGt                                 <  2:3416814/93‑1 (MQ=255)
                                                                cgGGCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATgg                            <  2:5062345/93‑1 (MQ=255)
                                                                 gggCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATggg                           <  1:1706895/93‑1 (MQ=255)
                                                                 gggCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATggg                           <  1:6674255/93‑1 (MQ=255)
                                                                  ggCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGc                          <  2:6391309/93‑1 (MQ=255)
                                                                  ggCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGc                          >  2:161186/1‑93 (MQ=255)
                                                                   gCTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCg                         >  1:203900/1‑93 (MQ=255)
                                                                    cTTCTCTTTCTCTGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGt                        >  1:3937172/1‑93 (MQ=255)
                                                                     ttctctttctctGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGtt                       >  2:294610/1‑93 (MQ=255)
                                                                      tctctttctctGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTa                      >  2:2132375/1‑93 (MQ=255)
                                                                      tctctttctctGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTa                      >  2:6048831/1‑93 (MQ=255)
                                                                      tctctttctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTa                      >  2:4585524/1‑93 (MQ=255)
                                                                       ctctttctctGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAt                     >  2:333378/1‑93 (MQ=255)
                                                                        tctttctctGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtc                    <  1:825772/93‑1 (MQ=255)
                                                                        tctttctctGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtc                    >  2:7268061/1‑93 (MQ=255)
                                                                        tctttctctGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtc                    >  2:139686/1‑93 (MQ=255)
                                                                         ctttctctGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtct                   >  2:5480700/1‑93 (MQ=255)
                                                                         ctttctctGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtct                   >  2:4077681/1‑93 (MQ=255)
                                                                         ctttctctGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtct                   <  2:3914606/93‑1 (MQ=255)
                                                                             ctctGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATcttctt               <  2:1291469/93‑1 (MQ=255)
                                                                             ctctGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATcttctt               <  2:857061/93‑1 (MQ=255)
                                                                             ctctGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATcttctt               <  1:6872959/93‑1 (MQ=255)
                                                                             ctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATcttctt               <  2:5561533/93‑1 (MQ=255)
                                                                                tGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGc            >  1:4196255/1‑93 (MQ=255)
                                                                                tGCGGGGGCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGc            <  1:1365634/93‑1 (MQ=255)
                                                                                     ggggCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATc        >  1:428628/1‑93 (MQ=255)
                                                                                       ggCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTa      <  1:3792545/93‑1 (MQ=255)
                                                                                       ggCTGTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTa      >  1:1156485/1‑93 (MQ=255)
                                                                                           aTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTc  <  2:2593216/93‑1 (MQ=255)
                                                                                           aTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTc  >  1:3665931/1‑93 (MQ=255)
                                                                                           |                                                                                            
TACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTC  >  NC_000913/704040‑704222

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