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.G99.1% 376.9 / ‑3.9 105coding (273/666 nt)yieHphosphoenolpyruvate and 6‑phosphogluconate phosphatase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (48/56);  total (48/57)
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.54e-01

Only 100 of 119 total aligned reads displayed.

TTTCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGGCT  >  NC_000913/3896908‑3897179
                                                                                                                                           |                                                                                                                                     
tttCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                        <  2:646730/139‑1 (MQ=255)
tttCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                        <  1:841439/139‑1 (MQ=255)
tttCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                        <  2:327566/139‑1 (MQ=255)
 ttCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAg                                                                                                                                      >  1:604422/1‑139 (MQ=255)
  tCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAgg                                                                                                                                     >  2:951515/1‑139 (MQ=255)
         tgtAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAGACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAgggggctgg                                                                                                                              >  2:1339933/1‑134 (MQ=255)
          gtAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGa                                                                                                                              <  2:55330/139‑1 (MQ=255)
           tAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAg                                                                                                                             >  2:562907/1‑139 (MQ=255)
           tAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAg                                                                                                                             >  2:524007/1‑139 (MQ=255)
           tAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAg                                                                                                                             >  2:289263/1‑139 (MQ=255)
             aaaCTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAgcg                                                                                                                           <  2:646168/139‑1 (MQ=255)
                 tGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTcc                                                                                                                       <  1:1470635/139‑1 (MQ=255)
                  gTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGTACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGTGCGCTCCt                                                                                                                      >  2:1089455/1‑139 (MQ=255)
                  gTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCt                                                                                                                      >  2:1108951/1‑139 (MQ=255)
                   tACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTg                                                                                                                     >  1:930140/1‑139 (MQ=255)
                   tACGAAATTATCGATATTGGTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTg                                                                                                                     >  1:1308248/1‑139 (MQ=255)
                        aaTTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGc                                                                                                                >  2:322496/1‑139 (MQ=255)
                          ttATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAg                                                                                                                             <  1:1687389/124‑1 (MQ=255)
                             tCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCa                                                                                                           <  1:322496/139‑1 (MQ=255)
                                  atTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAGACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                      <  1:1339933/139‑1 (MQ=255)
                                  atTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                      <  2:604422/139‑1 (MQ=255)
                                  atTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTgcg                                                                                                      >  2:202591/1‑139 (MQ=255)
                                    tGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                        >  2:1160807/1‑103 (MQ=255)
                                    tGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGcc                                                                                                    <  1:289998/139‑1 (MQ=255)
                                    tGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGcc                                                                                                    <  1:524007/139‑1 (MQ=255)
                                     gTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                                   <  1:1089455/139‑1 (MQ=255)
                                     gTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCa                                                                                                   <  1:289263/139‑1 (MQ=255)
                                       ttCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAt                                                                                                 >  2:1734351/1‑139 (MQ=255)
                                        tCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAgc                                                                                                                            <  1:1561718/111‑1 (MQ=255)
                                        tCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAgc                                                                                                                            >  2:1561718/1‑111 (MQ=255)
                                             tGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtgt                                                                                           >  2:1403306/1‑139 (MQ=255)
                                              gAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtgtg                                                                                          <  1:1108951/139‑1 (MQ=255)
                                              gAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                            <  2:1321924/137‑1 (MQ=255)
                                              gAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                            >  1:1321924/1‑137 (MQ=255)
                                                    ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTg                                                                                                                     <  2:737044/106‑1 (MQ=255)
                                                    ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTg                                                                                                                     >  1:737044/1‑106 (MQ=255)
                                                    ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCt                                                                                    >  1:322597/1‑139 (MQ=255)
                                                    ggTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCt                                                                                    >  2:109009/1‑139 (MQ=255)
                                                         tACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAgggg                                                                                                                                   <  1:988537/86‑1 (MQ=255)
                                                              tAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                        >  1:320256/1‑77 (MQ=255)
                                                              tAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAt                                                                                                 >  2:1325511/1‑116 (MQ=255)
                                                              tAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAt                                                                                                 <  1:1325511/116‑1 (MQ=255)
                                                              tAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCaa                                                                          >  1:611512/1‑139 (MQ=255)
                                                               aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgt                                                                                             <  2:1724255/119‑1 (MQ=255)
                                                               aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCaaa                                                                         >  2:330595/1‑139 (MQ=255)
                                                               aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCaaa                                                                         >  1:1732778/1‑139 (MQ=255)
                                                               aGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCaaa                                                                         >  2:662882/1‑139 (MQ=255)
                                                                gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGaa                                                                                                                                        >  1:1741073/1‑75 (MQ=255)
                                                                gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgt                                                                                               >  2:1510885/1‑116 (MQ=255)
                                                                gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAAt                                                                        >  1:1713816/1‑139 (MQ=255)
                                                                gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAAt                                                                        >  2:367353/1‑139 (MQ=255)
                                                                gCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAAt                                                                        >  2:1747/1‑139 (MQ=255)
                                                                gCGAAAACAGAAGCTGAACACGGTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgt                                                                                               <  1:1510885/116‑1 (MQ=255)
                                                                    aaaCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACa                                                                    <  1:367353/139‑1 (MQ=255)
                                                                    aaaCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACa                                                                    <  1:1734351/139‑1 (MQ=255)
                                                                     aaCAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACaa                                                                   >  1:1107685/1‑139 (MQ=255)
                                                                         gAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATg                                                               <  1:1608028/139‑1 (MQ=255)
                                                                         gAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATg                                                               <  1:1709741/139‑1 (MQ=255)
                                                                          aaGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACt                                                                                                         <  1:361881/96‑1 (MQ=255)
                                                                          aaGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACt                                                                                                         >  2:361881/1‑96 (MQ=255)
                                                                           aGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgca                                                             <  1:662882/139‑1 (MQ=255)
                                                                            gCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcag                                                            <  2:930140/139‑1 (MQ=255)
                                                                             cTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcagc                                                           <  2:322597/139‑1 (MQ=255)
                                                                             cTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcagc                                                           <  2:1089529/139‑1 (MQ=255)
                                                                                        tACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGc                                                <  2:1157339/139‑1 (MQ=255)
                                                                                         aCCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCa                                               <  2:1713816/139‑1 (MQ=255)
                                                                                          ccGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTc                                                                                                                   >  2:536944/1‑70 (MQ=255)
                                                                                          ccGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTc                                                                                                                   <  1:536944/70‑1 (MQ=255)
                                                                                           cGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAg                                             <  2:268345/139‑1 (MQ=255)
                                                                                            gTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAg                                             >  2:1374819/1‑138 (MQ=255)
                                                                                            gTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAg                                             <  1:1374819/138‑1 (MQ=255)
                                                                                             tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAAtgtgtg                                                                                            <  2:1585272/90‑1 (MQ=255)
                                                                                             tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcagc                                                           >  2:1537071/1‑123 (MQ=255)
                                                                                             tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATgcagc                                                           <  1:1537071/123‑1 (MQ=255)
                                                                                             tGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATggg                                                 >  1:1419946/1‑133 (MQ=255)
                                                                                              gCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTcc                                                                                                                       >  1:884896/1‑62 (MQ=255)
                                                                                              gCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTcc                                                                                                                       <  2:884896/62‑1 (MQ=255)
                                                                                               cAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGa                                         <  1:202591/139‑1 (MQ=255)
                                                                                                  aaGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGACCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAata                                      >  1:1154810/1‑139 (MQ=255)
                                                                                                    gTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCAc                                                                                                          <  2:811384/69‑1 (MQ=255)
                                                                                                    gTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCAc                                                                                                          >  1:811384/1‑69 (MQ=255)
                                                                                                    gTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGcc                                                                                                    <  1:593395/75‑1 (MQ=255)
                                                                                                       gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAgc                                                                                                                            <  2:538297/48‑1 (MQ=39)
                                                                                                       gCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAgc                                                                                                                            >  1:538297/1‑48 (MQ=255)
                                                                                                         tCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCa                               >  1:1662711/1‑139 (MQ=255)
                                                                                                            gCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAg                                             <  1:113858/122‑1 (MQ=255)
                                                                                                            gCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAg                                             >  2:113858/1‑122 (MQ=255)
                                                                                                               gTTCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACtt                         <  2:1154810/139‑1 (MQ=255)
                                                                                                                 tCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCact                             <  2:168980/133‑1 (MQ=255)
                                                                                                                 tCGATTCAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCact                             >  1:168980/1‑133 (MQ=255)
                                                                                                                       cAGAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATa                 <  1:562907/139‑1 (MQ=255)
                                                                                                                         gAACTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATaaa               <  2:1662711/139‑1 (MQ=255)
                                                                                                                          aaCTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTa                                                                                   <  1:193920/70‑1 (MQ=255)
                                                                                                                          aaCTGGAAGCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTa                                                                                   >  2:193920/1‑70 (MQ=255)
                                                                                                                                  gCCATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGc      <  2:1732778/139‑1 (MQ=255)
                                                                                                                                    cATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATggg                                                 <  2:10944/94‑1 (MQ=255)
                                                                                                                                    cATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATggg                                                 >  1:10944/1‑94 (MQ=255)
                                                                                                                                    cATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGc      <  1:825516/137‑1 (MQ=255)
                                                                                                                                    cATCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGc      >  2:825516/1‑137 (MQ=255)
                                                                                                                                      tCGAAGGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAGACTGTTCAGCGGCt  >  2:1364186/1‑139 (MQ=255)
                                                                                                                                           |                                                                                                                                     
TTTCAAAGGTGTAAAACTGTACGAAATTATCGATATTGTTTCCCTTGAACATGGTGTTACGTTAGCGAAAACAGAAGCTGAACACGTTTACCGTGCAGAAGTCGCTCGGCTGTTCGATTCAGAACTGGAAGCCATCGAAGGGGCTGGAGCGCTCCTGTCAGCGATCACTGCGCCAATGTGTGTGGTATCTAACGGCCCAAATAACAAAATGCAGCATTCTATGGGCAAGCTGAATATGTTGCACTACTTCCCGGATAAACTGTTCAGCGGCT  >  NC_000913/3896908‑3897179

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