Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,153,5730TG44.1% 104.6 / 36.9 68Y92D (TAC→GAC) pabC4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base T (1/37);  new base G (30/0);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.73e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAT  >  NC_000913/1153432‑1153722
                                                                                                                                             |                                                                                                                                                      
cGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACAt                                                                                                                                               <  1:423239/151‑1 (MQ=255)
      cAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACa                                                                                                                                         <  2:52165/151‑1 (MQ=255)
       aGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAg                                                                                                                                        <  2:203228/151‑1 (MQ=255)
           tGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCACACTGTAAAGCgg                                                                                                                                     >  2:135059/1‑150 (MQ=255)
            gCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGa                                                                                                                                    <  1:322528/150‑1 (MQ=255)
             cTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGACAGCACATTTaaaa                                                                                                                                         >  2:50277/1‑142 (MQ=255)
               tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGcacaattaaaa                                                                                                                                         >  2:191242/1‑135 (MQ=255)
                  cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCaa                                                                                                                             <  2:592410/151‑1 (MQ=255)
                  cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCa                                                                                                                              <  2:863978/150‑1 (MQ=255)
                   aGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAAc                                                                                                                            <  2:488232/151‑1 (MQ=255)
                       gTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCgg                                                                                                                        <  1:18780/151‑1 (MQ=255)
                             gatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATtct                                                                                                                   <  1:710004/150‑1 (MQ=255)
                              atTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAGCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATtctct                                                                                                                 <  1:623747/151‑1 (MQ=255)
                              atTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATtctct                                                                                                                 <  1:690448/151‑1 (MQ=255)
                               tttCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGAGAGCAAATTGAAAAACGGGGCGGGCAAGCGGGTtctc                                                                                                                  >  1:551521/1‑149 (MQ=255)
                                  ccTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCg                                                                                                              <  2:129600/150‑1 (MQ=255)
                                   cTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGCAAATTGAACAGGGGGCCGGCAACGCGGGGTCTCTcc                                                                                                               >  2:705834/1‑148 (MQ=255)
                                    tgtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAGAACACCTTTAAAAGCTGACCggggaagcgg                                                                                                                        >  1:463041/1‑130 (MQ=255)
                                       gACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACgg                                                                                                       <  2:673156/151‑1 (MQ=255)
                                        aCTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGc                                                                                                      <  1:28162/151‑1 (MQ=255)
                                          ttCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGCAAATTGAAAAACGGGACGGGCAGGCGGGTTCTCTTCCTTTCCGCtt                                                                                                    >  1:886280/1‑151 (MQ=255)
                                          ttCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCACAATGAACAGCGGAGCGGGAACGCGGGGTCTCTCCGTTACgg                                                                                                       >  1:317161/1‑148 (MQ=255)
                                             tGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGAAGCACATTTAAAAGCGGGCCGGCAACGCGGGTTTTCTCCCTTTCGGCtttt                                                                                                  >  2:798926/1‑148 (MQ=255)
                                              ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTa                                                                                                   <  1:461845/148‑1 (MQ=255)
                                               gCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGCACATTTAAAAACGGGCCGGCAACGCGGGGTTTCTCCCTTACGGGTTAt                                                                                                  >  2:843881/1‑148 (MQ=255)
                                                ccTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACCCTTAAGAGGGGGCCGGCGACGCGGGGTCTCTCCCTCACGGCTTAttc                                                                                                >  2:308738/1‑147 (MQ=255)
                                                      cTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGAAAGCACATTTAAAAGCGGGGCGGCAAAGCGGGGTCTCTCCGTTTCGGCTTATCCTGccc                                                                                           >  2:809203/1‑146 (MQ=255)
                                                       ttGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACATTGAACAGCGGACCGGCAACGCGGGTTCTCTCCCTTACGGCCTATCCTGCACAtt                                                                                        >  2:126248/1‑150 (MQ=255)
                                                        tGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTa                                                                                       <  2:361831/150‑1 (MQ=255)
                                                          aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGACAGCACATTTAACAGCGGACCGGCAACGCGGGTTCTCTCCCTTACGGCTTTTCCTGccc                                                                                           >  2:614282/1‑142 (MQ=255)
                                                           aCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAGCACAATGAAAAACGGGGCGGGCAAGCGGGTTCTCTCCCCTACGGGCTATCCTGccc                                                                                           >  1:408277/1‑141 (MQ=255)
                                                             aaGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCg                                                                                 <  1:537207/151‑1 (MQ=255)
                                                              agagATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGt                                                                                <  1:277957/151‑1 (MQ=255)
                                                              agagATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGt                                                                                <  2:293796/151‑1 (MQ=255)
                                                              agagATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGt                                                                                <  2:775100/151‑1 (MQ=255)
                                                                 gaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAAAACACTGAAAAACGGGACGGCGAAGCGGGGTTTCTCCGTTACCGGCTTTTCTTCACCATACGACCGgtt                                                                              >  1:103728/1‑150 (MQ=255)
                                                                  aTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCCCAATTAAAAACGGGCCGGCAACGCGGGGTCTCTCCGTTTCGGCTTATCCTGCACAATACGACCGttt                                                                              >  2:948613/1‑149 (MQ=255)
                                                                    gAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGt                                                                          <  2:883859/151‑1 (MQ=255)
                                                                        aCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACg                                                                      <  2:78559/151‑1 (MQ=255)
                                                                          gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAg                                                                    <  2:74123/151‑1 (MQ=255)
                                                                          gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGACAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAg                                                                    <  1:794110/151‑1 (MQ=255)
                                                                               cagcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGAt                                                               <  2:145779/151‑1 (MQ=255)
                                                                                        aGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACAATGAACAACGGGCCGGCCACGCGGGGTCTCTCCGTTACGGCTTTTTCCGCACACTACAACCGCTTTCGGAAAAAggggggttagttgg                                                       >  1:736424/1‑140 (MQ=255)
                                                                                         gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGgcg                                                     <  1:388028/151‑1 (MQ=255)
                                                                                         gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGgcg                                                     <  2:841254/151‑1 (MQ=255)
                                                                                           aaaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCACATTTAAAAACGGACCGGCAACGCGGGGTCTCTCCCTTACGGGCTTTTCTGCACATTACGACCGCTTTCGTAAAGAGGGGGGTATGTTgggg                                                     >  1:597314/1‑147 (MQ=255)
                                                                                            aaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGAAGCACATTTAACAGCGGACCGGGAACGCGGGTTCTCTCCCTTACGGCTTATTCTGCACATTACGACCCTTTGCGTAAAGAGGGGGTTACGTTGgcg                                                     >  2:632466/1‑148 (MQ=255)
                                                                                              aTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAg                                                <  1:841367/151‑1 (MQ=255)
                                                                                               tGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACATTGAACAGCGGGCCGGCCACGCGGGTTCTCTCCGTTACGGCTTTTCCTGCACATTACGACCCTTTGCGGAACGAGGGGGTTACGCTGGCGCTAAGc                                               >  2:285191/1‑151 (MQ=255)
                                                                                                 gtgtGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGccc                                             <  1:614214/151‑1 (MQ=255)
                                                                                                    tgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACAATTAACAGCGGGCCGGCAAAGCGGGTTTTCTCCGTTACGGCTTATCCTGCACATTAAGAACCGTTGCGGAAAGAGGGGGGTACGTTTGCGCCAAAccc                                             >  2:232524/1‑148 (MQ=255)
                                                                                                     gCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGt                                          <  2:458527/150‑1 (MQ=255)
                                                                                                      cTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGGACATTGTAAAGCGGACCGGCCACGCGGAGTCTCTCCCGTTCGGCGTATCCTGCGCAATACTACCGGTTGCGTAAAGAGGGGGGTACGCTTGCGCTAAGacccgcg                                         >  2:161891/1‑146 (MQ=255)
                                                                                                      cTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAGCACATTTAACAGCGGACCGGCAACGCGGATTCTCTCCGTTTCGGCTTATCCTGCAAATTACGACCCCTTTCGGTACGAGGGGGGTAAGGTGgcggttaggcccggg                                         >  2:269972/1‑138 (MQ=255)
                                                                                                       tGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCACCTTGAACAGCGGACCGGCAACGCGGGGTCTCTCCGTTTCGGGCTATCCTTCACATTAAGACCGCTTTCGGAAAGAGGGGGGTACGTTGGGGCTTAGccc                                             >  2:535620/1‑145 (MQ=255)
                                                                                                            gTCGTGATCAGTCGCGGTAGTGGCGGGCGCGGGGGCAGCACATTGTACAACGGCCCGGCAACGCGGAGTGCGTCCGTTACCGGCTTTTCTGCACATTAAGACCGGTTTCGTCAAGCGGGGGGTTAGTTGGCGCCAAGACCGGGGCgg                                      >  1:100679/1‑147 (MQ=255)
                                                                                                            gTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTgg                                  <  2:398547/151‑1 (MQ=255)
                                                                                                              cGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGGACACTGAAAAAAGGGCCGGCCACGCGGGGTTTCTCCCGTTCGGCGTTTTCTTCACATTACGCCCGGTTTGGTAAAGAGGGGGGTAAGTTGGCGCTAAGCCCCGGGCgggggggg                                >  1:676909/1‑151 (MQ=255)
                                                                                                              cGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAACACATTTAACAGCGGGCCGGCAACGCGGATTCTCTCCCTTACGGCTTTTTCTGCACATTACGACCCTTTGCGTAACGAGGGGGTTAAGTTTGCGCTAAACCCCGGGCGGCTgggg                                >  2:493402/1‑151 (MQ=255)
                                                                                                                 gATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGgcg                              <  2:540934/150‑1 (MQ=255)
                                                                                                                   tCAGTCGCGGTAGTGGCGGGGGGGGGGGCAACACAATGAAAAAAGGGCCGGGAACGCGGGGTCTCTCCCCTACGGCGTTTTCTGCACATTACGACCGCTTGCGTAAAGAGGGGGGTACGTTGGGGCTAAGACCGGGGCGGCGgggggg                              >  1:844875/1‑146 (MQ=255)
                                                                                                                   tCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACATTGTACAGCGGGCCGGCGAAGCGGGGTCTCTCCGTTACGGCGTATCCTGCCCATTTAGAACGGTTGCGTAAAGAGGGGGGTAAGTTGGGGGTAAGCCCCGGGCGGCTgggggg                              >  1:428305/1‑146 (MQ=255)
                                                                                                                    cAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAAt                          <  1:51748/151‑1 (MQ=255)
                                                                                                                                   cgggcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGCGATTACGTTGGCGCTAAGCCCGGTACGGCTGGGGCGCAATCCTCATCTTGCAgct           <  1:272299/151‑3 (MQ=255)
                                                                                                                                    gggcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTa          <  2:214763/151‑1 (MQ=255)
                                                                                                                                       cgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTa       <  1:814979/151‑1 (MQ=255)
                                                                                                                                         aGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTaa      >  1:588975/1‑150 (MQ=255)
                                                                                                                                            gTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAt  <  2:730834/151‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                      
CGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAT  >  NC_000913/1153432‑1153722

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