Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,897,050 +G coding (277/666 nt) yieH → phosphoenolpyruvate and 6‑phosphogluconate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,897,0461.G97.6% 289.3 / ‑0.9 83coding (273/666 nt)yieHphosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (36/45);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 103 total aligned reads displayed.

TTTCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑G‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGG  >  NC_000913/3896908‑3897177
                                                                                                                                           |                                                                                                                                    
tttCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                       >  2:854502/1‑139 (MQ=255)
  tCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑g‑g                                                                                                                                   >  2:944500/1‑139 (MQ=255)
              aaCTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAgcgc                                                                                                                        <  1:1036727/139‑1 (MQ=255)
                   tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑GTGGGCTGGAGCGCTCCTg                                                                                                                   <  2:855729/139‑1 (MQ=255)
                             tCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         <  1:854502/139‑1 (MQ=255)
                             tCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         <  1:306368/139‑1 (MQ=255)
                             tCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         <  2:819004/139‑1 (MQ=255)
                                  atTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                    >  2:1149064/1‑139 (MQ=255)
                                    tGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGAt                                                                                                           <  1:1118984/130‑1 (MQ=255)
                                    tGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGAt                                                                                                           >  2:1118984/1‑130 (MQ=255)
                                    tGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGcc                                                                                                  <  1:291096/139‑1 (MQ=255)
                                     gTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                                 <  2:967217/139‑1 (MQ=255)
                                     gTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                                 <  2:750179/139‑1 (MQ=255)
                                       ttCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAt                                                                                               <  1:1149064/139‑1 (MQ=255)
                                         cccTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTACGATTCAGAACTGGAAGCCATCGAA‑GCGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgt                                                                                             >  2:829972/1‑139 (MQ=255)
                                          ccTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCTGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtg                                                                                            <  2:1359743/139‑1 (MQ=255)
                                             tGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtgt                                                                                         >  1:975202/1‑139 (MQ=255)
                                              gAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         <  2:1034900/122‑1 (MQ=255)
                                              gAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         >  1:1034900/1‑122 (MQ=255)
                                               aaCATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTgg                                                                                       >  1:1021336/1‑139 (MQ=255)
                                                    ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAGCTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          <  2:784690/131‑1 (MQ=255)
                                                    ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                       >  1:382328/1‑87 (MQ=255)
                                                    ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                       <  1:750410/87‑1 (MQ=255)
                                                    ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                       <  2:382328/87‑1 (MQ=255)
                                                    ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                       >  2:750410/1‑87 (MQ=255)
                                                    ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          >  1:784690/1‑131 (MQ=255)
                                                    ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCt                                                                                  >  2:1124641/1‑139 (MQ=255)
                                                      tgtTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                       <  2:1241679/85‑1 (MQ=255)
                                                      tgtTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                       >  1:1241679/1‑85 (MQ=255)
                                                        ttACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                       <  1:1098769/83‑1 (MQ=255)
                                                        ttACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                       >  2:1098769/1‑83 (MQ=255)
                                                        ttACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACg                                                                              <  2:720489/139‑1 (MQ=255)
                                                         tACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACgg                                                                             >  2:1360272/1‑139 (MQ=255)
                                                          aCGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGc                                                                            <  1:1346557/139‑1 (MQ=255)
                                                            gTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGccc                                                                          <  2:641775/139‑1 (MQ=255)
                                                              tAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgt                                                                                             >  2:366172/1‑118 (MQ=255)
                                                              tAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgt                                                                                             <  1:366172/118‑1 (MQ=255)
                                                               aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                       >  1:903915/1‑76 (MQ=255)
                                                               aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                       <  2:903915/76‑1 (MQ=255)
                                                               aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCaaa                                                                       <  1:944500/139‑1 (MQ=255)
                                                                gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          <  1:439386/119‑1 (MQ=255)
                                                                gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          >  2:439386/1‑119 (MQ=255)
                                                                gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑g‑ggggc                                                                                                                               <  2:1091974/81‑3 (MQ=255)
                                                                gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑g‑ggggc                                                                                                                               >  1:1091974/1‑79 (MQ=255)
                                                                gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑g‑gggg                                                                                                                                <  2:353605/80‑2 (MQ=255)
                                                                gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑g‑gggg                                                                                                                                >  1:353605/1‑79 (MQ=255)
                                                                     aaCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgt                                                                                             >  1:701323/1‑111 (MQ=255)
                                                                     aaCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgt                                                                                             <  2:701323/111‑1 (MQ=255)
                                                                      aCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACaaa                                                                >  2:244215/1‑139 (MQ=255)
                                                                            gCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcag                                                          <  2:975202/139‑1 (MQ=255)
                                                                            gCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcag                                                          <  1:1360272/139‑1 (MQ=255)
                                                                            gCTGAACACGTTTACCGTGCAGAAGCCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          <  2:788418/107‑1 (MQ=255)
                                                                            gCTGAACACGTTTACCGTGCAGAAGCCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                          >  1:788418/1‑107 (MQ=255)
                                                                                       ttACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGc                                                                            >  1:509661/1‑110 (MQ=255)
                                                                                       ttACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGc                                                                            <  2:509661/110‑1 (MQ=255)
                                                                                       ttACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATggg                                               >  1:1148117/1‑139 (MQ=255)
                                                                                          ccGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCaa                                            <  1:1124641/139‑1 (MQ=255)
                                                                                            gTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACaaa                                                                >  2:240994/1‑117 (MQ=255)
                                                                                            gTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACaaa                                                                <  1:240994/117‑1 (MQ=255)
                                                                                             tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                       <  1:206172/46‑1 (MQ=255)
                                                                                             tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                       >  2:206172/1‑46 (MQ=255)
                                                                                             tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                    >  2:49275/1‑80 (MQ=255)
                                                                                             tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                    <  1:49275/80‑1 (MQ=255)
                                                                                             tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACa                                                                  >  2:867793/1‑114 (MQ=255)
                                                                                             tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACa                                                                  <  1:867793/114‑1 (MQ=255)
                                                                                             tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                              >  2:740354/1‑134 (MQ=255)
                                                                                               cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGa                                                                                                                            >  1:983506/1‑54 (MQ=255)
                                                                                               cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         <  2:374713/73‑1 (MQ=255)
                                                                                               cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                         >  1:374713/1‑73 (MQ=255)
                                                                                               cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAt                                                                                               <  1:852904/83‑1 (MQ=255)
                                                                                               cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAt                                                                                               >  2:852904/1‑83 (MQ=255)
                                                                                                  aaGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAata                                    >  2:672121/1‑139 (MQ=255)
                                                                                                    gTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑g                                                                                                                                     >  2:176017/1‑40 (MQ=255)
                                                                                                    gTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑g                                                                                                                                     <  1:176017/40‑1 (MQ=255)
                                                                                                     tCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGt                                                                                      >  1:875412/1‑86 (MQ=255)
                                                                                                     tCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGt                                                                                      <  2:875412/86‑1 (MQ=255)
                                                                                                       gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATg                                  <  1:212597/136‑1 (MQ=255)
                                                                                                       gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATg                                  >  2:212597/1‑136 (MQ=255)
                                                                                                       gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTg                               <  2:490626/139‑1 (MQ=255)
                                                                                                       gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTACGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTg                               >  1:798394/1‑139 (MQ=255)
                                                                                                            gCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgt                                                                                             <  2:1201121/72‑1 (MQ=255)
                                                                                                            gCTGTTCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgt                                                                                             >  1:1201121/1‑72 (MQ=255)
                                                                                                                ttCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACt                                                                                                       >  1:315124/1‑58 (MQ=255)
                                                                                                                ttCGATTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACt                                                                                                       <  2:315124/58‑1 (MQ=255)
                                                                                                                    aTTCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCgg                  <  2:208509/139‑1 (MQ=255)
                                                                                                                     ttCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCt                                                                                                                       <  2:1215638/37‑1 (MQ=255)
                                                                                                                     ttCAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCt                                                                                                                       >  1:1215638/1‑37 (MQ=255)
                                                                                                                       cAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                              >  2:757906/1‑108 (MQ=255)
                                                                                                                       cAGAACTGGAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                              <  1:757906/108‑1 (MQ=255)
                                                                                                                               gAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcagc                                                         >  1:411736/1‑89 (MQ=255)
                                                                                                                               gAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcagc                                                         <  2:411736/89‑1 (MQ=255)
                                                                                                                               gAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcag                                                          <  2:1220808/88‑1 (MQ=255)
                                                                                                                               gAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcag                                                          >  1:1220808/1‑88 (MQ=255)
                                                                                                                               gAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                              >  2:372430/1‑100 (MQ=255)
                                                                                                                               gAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                              <  1:372430/100‑1 (MQ=255)
                                                                                                                               gAAGCCATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCACCAATGTGGGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTATATGGGCAAGCTGAATATGTTGCACTACTTCACGGATAAACTGTTc       >  1:856916/1‑139 (MQ=255)
                                                                                                                                    cATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                                 <  2:1068309/44‑1 (MQ=255)
                                                                                                                                    cATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                                 >  1:1068309/1‑44 (MQ=255)
                                                                                                                                    cATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCgg  <  2:798394/139‑1 (MQ=255)
                                                                                                                                    cATCGAAGG‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCgg  <  1:152621/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                    
TTTCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAA‑G‑GGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGG  >  NC_000913/3896908‑3897177

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