Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 704,236 T→C L357S (TTG→TCG)  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,2360TC99.6% 826.3 / ‑5.2 239L357S (TTG→TCG) nagEPTS system, N‑acetylglucosamine‑specific enzyme IIABC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (140/98);  minor base A (0/1);  total (140/99)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 247 total aligned reads displayed.

TGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTTGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAG  >  NC_000913/704146‑704325
                                                                                          |                                                                                         
tGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCgg                                                                                         <  1:1854788/93‑1 (MQ=255)
 gTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGt                                                                                        <  1:5093570/93‑1 (MQ=255)
   ataACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTa                                                                                      >  1:8526932/1‑93 (MQ=255)
   ataACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTa                                                                                      >  2:1886538/1‑93 (MQ=255)
    taACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTAt                                                                                     <  1:769867/93‑1 (MQ=255)
    taACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTAt                                                                                     >  2:7668342/1‑93 (MQ=255)
     aaCCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATc                                                                                    <  1:78218/93‑1 (MQ=255)
         tGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCa                                                                                <  1:4778865/93‑1 (MQ=255)
          gccggcgtCCAGCCAGAACGTCTCGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCAt                                                                               <  1:4810263/93‑1 (MQ=255)
          gccggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCAt                                                                               >  2:2737294/1‑93 (MQ=255)
            cggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGt                                                                             <  2:2443922/93‑1 (MQ=255)
              gccgccAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTc                                                                           >  2:5786759/1‑93 (MQ=255)
               ccgccAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCa                                                                          >  1:6118988/1‑93 (MQ=255)
                cgccAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCaa                                                                         <  2:2910372/93‑1 (MQ=255)
                 gccagccaGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAAc                                                                        <  2:6188170/93‑1 (MQ=255)
                 gccagccaGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAAc                                                                        >  2:4893590/1‑93 (MQ=255)
                  ccagccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAAcc                                                                       <  2:260854/93‑1 (MQ=255)
                    agccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTg                                                                     >  1:408176/1‑93 (MQ=255)
                      ccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaa                                                                   >  2:4777700/1‑93 (MQ=255)
                      ccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaa                                                                   >  2:6485537/1‑93 (MQ=255)
                      ccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaa                                                                   <  1:6010067/93‑1 (MQ=255)
                       cagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaaa                                                                  >  1:2881592/1‑93 (MQ=255)
                        agtaCGTCTGGATGCTGCTGGTGATGTGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGcaaa                                                                 >  2:4959467/4‑93 (MQ=255)
                        agAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCCACCTGcaac                                                                 >  1:4456505/1‑92 (MQ=255)
                        agAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaaaa                                                                 <  1:921935/93‑1 (MQ=255)
                        agAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaaaa                                                                 <  2:5685013/93‑1 (MQ=255)
                        agAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaaaa                                                                 <  2:655914/93‑1 (MQ=255)
                          aaCGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACg                                                               <  2:8013348/93‑1 (MQ=255)
                           aCGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGc                                                              >  2:6898732/1‑93 (MQ=255)
                           aCGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGc                                                              <  1:7636006/93‑1 (MQ=255)
                             gTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCg                                                            >  1:3797262/1‑93 (MQ=255)
                             gTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCg                                                            >  1:127381/1‑93 (MQ=255)
                              tCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCgg                                                           >  2:3328688/1‑93 (MQ=255)
                              tCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCgg                                                           >  1:6457578/1‑93 (MQ=255)
                                  gATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCg                                                       >  2:7446638/1‑93 (MQ=255)
                                   aTGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGt                                                      <  2:5134622/93‑1 (MQ=255)
                                   aTGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGt                                                      <  2:3943214/93‑1 (MQ=255)
                                     gctgctGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGa                                                    <  1:4359423/93‑1 (MQ=255)
                                      ctgctgGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGaa                                                   <  1:1950880/93‑1 (MQ=255)
                                       tgctgGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAg                                                  >  2:5311426/1‑93 (MQ=255)
                                        gctgGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAAGCTGAAAACGCCGGGTCGTGAAGa                                                 >  2:4755728/1‑93 (MQ=255)
                                        gctgGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGa                                                 <  2:3507480/93‑1 (MQ=255)
                                         ctgGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGAt                                                >  2:4540369/1‑93 (MQ=255)
                                          tgGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATa                                               <  2:1696556/93‑1 (MQ=255)
                                              gATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaaga                                           >  1:323802/1‑93 (MQ=255)
                                              gATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaaga                                           >  1:163346/1‑93 (MQ=255)
                                               aTGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaagaa                                          <  1:3073255/93‑1 (MQ=255)
                                                tGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaagaag                                         >  1:7180201/1‑93 (MQ=255)
                                                  ggCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGAc                                       >  2:3447609/1‑93 (MQ=255)
                                                  ggCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGAc                                       >  1:7299050/1‑93 (MQ=255)
                                                  ggCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGAc                                       <  2:5596649/93‑1 (MQ=255)
                                                    ggTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACga                                     <  1:1522465/92‑1 (MQ=255)
                                                    cGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACga                                     <  1:2738666/93‑1 (MQ=255)
                                                     gTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACgag                                    >  2:2157695/1‑93 (MQ=255)
                                                      ttATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACgaga                                   >  1:8217605/1‑93 (MQ=255)
                                                        aTCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATc                                 >  2:3797832/1‑93 (MQ=255)
                                                        aTCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATc                                 <  1:6990627/93‑1 (MQ=255)
                                                         tcttctTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCg                                >  1:8697664/1‑93 (MQ=255)
                                                         tcttctTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCg                                <  1:7367565/93‑1 (MQ=255)
                                                          cttcttCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGt                               <  1:4089091/93‑1 (MQ=255)
                                                           ttcttcGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGtt                              <  2:1232797/93‑1 (MQ=255)
                                                           ttcttcGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGtt                              >  1:2612559/1‑93 (MQ=255)
                                                           ttcttcGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGtt                              >  1:3712965/1‑93 (MQ=255)
                                                            tcttcGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTa                             >  1:3286363/1‑93 (MQ=255)
                                                             cttcGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTAc                            >  1:6184012/1‑93 (MQ=255)
                                                                 gCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaa                        <  1:5791000/93‑1 (MQ=255)
                                                                  cTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaag                       >  2:1462214/1‑93 (MQ=255)
                                                                   tATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaaga                      >  2:6761304/1‑93 (MQ=255)
                                                                   tATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaaga                      <  1:1142929/93‑1 (MQ=255)
                                                                   tATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaaga                      >  1:2168683/1‑93 (MQ=255)
                                                                    aTCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaagaa                     >  2:8084829/1‑93 (MQ=255)
                                                                    aTCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaagaa                     >  2:5978599/1‑93 (MQ=255)
                                                                      cttcTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGc                   <  1:1694058/90‑1 (MQ=255)
                                                                      cTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGCTCGTCACCGCAGAAGc                   >  1:8014952/1‑93 (MQ=255)
                                                                      cTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGc                   >  1:6928554/1‑93 (MQ=255)
                                                                       tACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGcc                  <  2:5508718/93‑1 (MQ=255)
                                                                       tACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGcc                  >  2:8671187/1‑93 (MQ=255)
                                                                        aCTTCGTGGTGTTCAGTTCGGTTATCCGCATCTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCa                 <  1:8114323/93‑1 (MQ=255)
                                                                          ttCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCagc               <  2:8390968/93‑3 (MQ=255)
                                                                          ttCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCagc               <  2:1452987/93‑3 (MQ=255)
                                                                             gtggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCcagcaga            >  1:5480442/1‑92 (MQ=255)
                                                                             gtggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCagcagc            >  1:1107658/1‑93 (MQ=255)
                                                                             gtggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCagcagc            >  1:519908/1‑93 (MQ=255)
                                                                             gtggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCagcagc            >  2:6820298/1‑93 (MQ=255)
                                                                              tggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAgcagca           >  1:1027235/1‑93 (MQ=255)
                                                                              tggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAgcagca           <  1:3344263/93‑1 (MQ=255)
                                                                              cggtgTTAAGTTCGGTTATCCGCATGTTCACCCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAgcagca           <  2:7257535/92‑1 (MQ=255)
                                                                               ggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAGCAGCaa          <  1:5852061/93‑1 (MQ=255)
                                                                               ggtgTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAGCAGCaa          >  2:164472/1‑93 (MQ=255)
                                                                                gtgtTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAGCAGCAac         >  2:7419018/1‑93 (MQ=255)
                                                                                  tgtCAGTTGTGTTATCCGCATCTCCAACCTGAAAACGCCTGGTCGTGAAGATAAGGAAGACGAGATCGTTACTGAAGAAGCCAGCAGCAacac       <  1:30594/91‑1 (MQ=255)
                                                                                  gtTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAGCAGCAacac       >  2:5734077/1‑93 (MQ=255)
                                                                                  gtTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAGCAGCAacac       <  2:1789511/93‑1 (MQ=255)
                                                                                  gtTCAGTTAGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAGCAGCAacac       <  2:4475531/93‑1 (MQ=255)
                                                                                   ttCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGGGATCGTTACTGAAGAAGCCAGCAGCAACACt      >  1:5936019/1‑93 (MQ=255)
                                                                                   ttCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAGCAGCAACACt      >  1:5863705/1‑93 (MQ=255)
                                                                                     cAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAGCAGCAACACTga    >  2:7698384/1‑93 (MQ=255)
                                                                                      aGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAGCAGCAACACTgaa   >  1:6906281/1‑93 (MQ=255)
                                                                                       gTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAGCAGCAACACTgaag  >  2:1314291/1‑93 (MQ=255)
                                                                                       gTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAGCAGCAACACTgaag  >  2:2356795/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTTGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACTGAAG  >  NC_000913/704146‑704325

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