Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,912,542 T→C D96G (GAT→GGT)  glmU ← N‑acetyl glucosamine‑1‑phosphate uridyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,912,5420TC98.1% 483.7 / ‑0.9 155D96G (GAT→GGT) glmUN‑acetyl glucosamine‑1‑phosphate uridyltransferase
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base C (91/61);  total (94/61)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 162 total aligned reads displayed.

TGCCACCCTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCATCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTT  >  NC_000913/3912453‑3912630
                                                                                         |                                                                                        
tGCCACCCTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTcaccat                                                                                       <  1:2432384/93‑1 (MQ=255)
 gCCACCCTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCaccatc                                                                                      <  2:233151/93‑1 (MQ=255)
 gCCACCCTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCaccatc                                                                                      <  2:424406/93‑1 (MQ=255)
 gCCACCCTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCaccatc                                                                                      <  1:2049662/93‑1 (MQ=255)
  ccACCCTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCg                                                                                     <  2:2673532/93‑1 (MQ=255)
  ccACCCTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCg                                                                                     >  2:4487979/1‑93 (MQ=255)
     cccTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCa                                                                                  >  1:2621150/1‑93 (MQ=255)
       cTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCaaa                                                                                >  1:118109/1‑93 (MQ=255)
         gCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCaaaga                                                                              >  2:1877394/1‑93 (MQ=255)
         gCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCaaaga                                                                              >  2:4535871/1‑93 (MQ=255)
         gCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCaaaga                                                                              <  2:3607683/93‑1 (MQ=255)
            gTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCaaagaaag                                                                           <  1:3246779/93‑1 (MQ=255)
             tttAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAgg                                                                          <  2:2434440/93‑1 (MQ=255)
             tttAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAgg                                                                          <  1:44062/93‑1 (MQ=255)
             tttAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAgg                                                                          <  1:1676636/93‑1 (MQ=255)
              ttAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGt                                                                         <  1:3708783/93‑1 (MQ=255)
              ttAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGt                                                                         >  1:443662/1‑93 (MQ=255)
                aGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGc                                                                       <  1:4252981/93‑1 (MQ=255)
                aGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGc                                                                       >  1:773689/1‑93 (MQ=255)
                  cATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCgg                                                                     <  2:2164166/93‑1 (MQ=255)
                   aTCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGc                                                                    >  2:4270849/1‑93 (MQ=255)
                   aTCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGc                                                                    <  1:257132/93‑1 (MQ=255)
                   aTCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGc                                                                    <  2:4530155/93‑1 (MQ=255)
                     cACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCct                                                                  >  2:1536406/1‑93 (MQ=255)
                     cACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCct                                                                  <  2:1253761/93‑1 (MQ=255)
                      aCGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCctg                                                                 >  1:2108541/1‑93 (MQ=255)
                      aCGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCctg                                                                 >  1:2565332/1‑93 (MQ=255)
                       cGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCctgc                                                                >  1:2227864/1‑93 (MQ=255)
                        gCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCATCATCGGCAAAAAACGGTGCGGCctgct                                                               >  1:3310759/1‑93 (MQ=255)
                        gCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCctgct                                                               >  1:2643457/1‑93 (MQ=255)
                        gCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCctgct                                                               <  1:3078013/93‑1 (MQ=255)
                              gCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCtgcattg                                                         >  2:859252/1‑93 (MQ=255)
                               cTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAGGGTGCGGCCTGCtgcattgc                                                        <  2:4572541/93‑1 (MQ=255)
                               cTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCtgcattgc                                                        >  1:3192658/1‑93 (MQ=255)
                                 ggAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCtgcattgcat                                                      <  1:3601050/93‑1 (MQ=255)
                                 ggAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCtgcattgcat                                                      >  2:2163903/1‑93 (MQ=255)
                                 ggAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCtgcattgcat                                                      >  2:379216/1‑93 (MQ=255)
                                  gAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATg                                                     >  1:2493375/1‑93 (MQ=255)
                                   aGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGa                                                    <  2:3338807/93‑1 (MQ=255)
                                   aGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGa                                                    <  2:2932507/93‑1 (MQ=255)
                                      gtTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGAccc                                                 >  2:1843516/1‑93 (MQ=255)
                                       tttCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCg                                                >  1:1920352/1‑93 (MQ=255)
                                           gACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCATCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTAcc                                            >  1:4018272/1‑93 (MQ=255)
                                            aCAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTAccc                                           <  1:4257486/93‑1 (MQ=255)
                                            aCAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTAccc                                           <  1:599619/93‑1 (MQ=255)
                                            aCAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTAccc                                           <  2:1751869/93‑1 (MQ=255)
                                             cAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCa                                          >  2:3285110/1‑93 (MQ=255)
                                                agaTCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAgct                                       <  1:592397/93‑1 (MQ=255)
                                                   tCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCCgctgct                                    >  2:4741750/1‑93 (MQ=255)
                                                    cAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGctc                                   >  2:3258666/1‑93 (MQ=255)
                                                     aGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGctct                                  <  1:1095543/93‑1 (MQ=255)
                                                     aGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGctct                                  >  1:363718/1‑93 (MQ=255)
                                                     aGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGctct                                  <  2:3969743/93‑1 (MQ=255)
                                                     aGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGGGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGctct                                  >  2:1208075/1‑93 (MQ=255)
                                                      gCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTg                                 >  2:2834039/1‑93 (MQ=255)
                                                      gCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTg                                 >  2:2275804/1‑93 (MQ=255)
                                                      gCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTg                                 >  2:2868258/1‑93 (MQ=255)
                                                      gCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTg                                 >  2:3368860/1‑93 (MQ=255)
                                                        ggCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTCCCCAGCTGCTCTGcc                               >  2:4055132/1‑93 (MQ=255)
                                                        ggCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGcc                               >  1:1190894/1‑93 (MQ=255)
                                                         gCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCt                              <  2:232625/93‑1 (MQ=255)
                                                         gCACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCt                              >  1:869493/1‑93 (MQ=255)
                                                          cACGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTg                             <  1:1874270/93‑1 (MQ=255)
                                                          cACGTCGCCGTAGAGCATTAAAATGTCGTCACCATCGGCAAAGAAAGGTGCGGCGTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTg                             >  1:2424487/1‑93 (MQ=255)
                                                            cGTCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGaa                           <  2:1023652/93‑1 (MQ=255)
                                                              tCGCCGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGc                         >  1:1465857/1‑93 (MQ=255)
                                                                  cGTAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCAccc                     >  1:1723935/1‑93 (MQ=255)
                                                                    tAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAg                   >  2:1341719/1‑93 (MQ=255)
                                                                    tAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAg                   >  2:861841/1‑93 (MQ=255)
                                                                    tAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAg                   >  2:1010758/1‑93 (MQ=255)
                                                                    tAGAGCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAg                   <  1:1496977/93‑1 (MQ=255)
                                                                     agagCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGt                  >  2:4274868/1‑93 (MQ=255)
                                                                     agagCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGt                  <  2:1039630/93‑1 (MQ=255)
                                                                     agagCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGt                  <  1:2196222/93‑1 (MQ=255)
                                                                      gagCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCCGtt                 >  2:981741/1‑93 (MQ=255)
                                                                      gagCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGtt                 >  1:4310192/1‑93 (MQ=255)
                                                                      gagCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGtt                 >  1:1794508/1‑93 (MQ=255)
                                                                       agCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCCTTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACACAGTTa                >  2:2734344/1‑93 (MQ=255)
                                                                        gCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTaa               >  1:722242/1‑93 (MQ=255)
                                                                        gCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTaa               <  1:1513189/93‑1 (MQ=255)
                                                                        gCATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTaa               >  1:1620036/1‑93 (MQ=255)
                                                                         cATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAg              >  2:2535721/1‑93 (MQ=255)
                                                                         cATTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAg              >  2:2400989/1‑93 (MQ=255)
                                                                          aTTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAgg             >  1:4175448/1‑93 (MQ=255)
                                                                          aTTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAgg             >  1:3889252/1‑93 (MQ=255)
                                                                          aTTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAgg             >  1:3570917/1‑93 (MQ=255)
                                                                          aTTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAgg             >  2:4038827/1‑93 (MQ=255)
                                                                          aTTAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAgg             <  2:2170392/93‑1 (MQ=255)
                                                                           ttAAAATGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGt            >  1:1029818/1‑93 (MQ=255)
                                                                             aaaaTGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTg          >  2:671097/1‑93 (MQ=255)
                                                                              aaaTGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTgt         >  2:1503737/1‑93 (MQ=255)
                                                                              aaaTGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTgt         >  1:3040576/1‑93 (MQ=255)
                                                                               aaTGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTgtc        <  1:673305/93‑1 (MQ=255)
                                                                                aTGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTgtcg       <  2:3423413/93‑1 (MQ=255)
                                                                                 tGTCTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTgtcgt      >  1:2931837/1‑93 (MQ=255)
                                                                                    cTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCtt   <  1:3509690/93‑1 (MQ=255)
                                                                                    cTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCtt   <  2:2196666/93‑1 (MQ=255)
                                                                                    cTTCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCtt   <  2:1939318/93‑1 (MQ=255)
                                                                                     ttCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCttt  >  1:858827/1‑93 (MQ=255)
                                                                                     ttCACCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGATTGTCGTCttt  >  1:2410649/1‑93 (MQ=255)
                                                                                         |                                                                                        
TGCCACCCTGCGGTTTAGCATCACGCAGACGCTGGAGTGTTTCGACAGAGATCAGCGGCACGTCGCCGTAGAGCATTAAAATGTCTTCATCATCGGCAAAGAAAGGTGCGGCCTGCTGCATTGCATGACCCGTACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTT  >  NC_000913/3912453‑3912630

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