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

CGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAT  >  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)
                                       gACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑Cgg                                                                                                       <  2:673156/151‑1 (MQ=255)
                                        aCTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGc                                                                                                      <  1:28162/151‑1 (MQ=255)
                                          ttCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAGCAAATTGAAAAACGGGACGGGCAGGCGGGTTCTCTTCCTTT‑CCGCtt                                                                                                    >  1:886280/1‑151 (MQ=255)
                                          ttCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCACAATGAACAGCGGAGCGGGAACGCGGGGTCTCTCCGTTA‑Cgg                                                                                                       >  1:317161/1‑148 (MQ=255)
                                             tGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGAAGCACATTTAAAAGCGGGCCGGCAACGCGGGTTTTCTCCCTTT‑CGGCtttt                                                                                                  >  2:798926/1‑148 (MQ=255)
                                              ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTa                                                                                                   <  1:461845/148‑1 (MQ=255)
                                               gCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGCACATTTAAAAACGGGCCGGCAACGCGGGGTTTCTCCCTTA‑CGGGTTAt                                                                                                  >  2:843881/1‑148 (MQ=255)
                                                ccTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACCCTTAAGAGGGGGCCGGCGACGCGGGGTCTCTCCCTCA‑CGGCTTAttc                                                                                                >  2:308738/1‑147 (MQ=255)
                                                      cTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGAAAGCACATTTAAAAGCGGGGCGGCAAAGCGGGGTCTCTCCGTTT‑CGGCTTATCCTGccc                                                                                           >  2:809203/1‑146 (MQ=255)
                                                       ttGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACATTGAACAGCGGACCGGCAACGCGGGTTCTCTCCCTTA‑CGGCCTATCCTGCACAtt                                                                                        >  2:126248/1‑150 (MQ=255)
                                                        tGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTa                                                                                       <  2:361831/150‑1 (MQ=255)
                                                          aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGACAGCACATTTAACAGCGGACCGGCAACGCGGGTTCTCTCCCTTA‑CGGCTTTTCCTGccc                                                                                           >  2:614282/1‑142 (MQ=255)
                                                           aCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAGCACAATGAAAAACGGGGCGGGCAAGCGGGTTCTCTCCCCTA‑CGGGCTATCCTGccc                                                                                           >  1:408277/1‑141 (MQ=255)
                                                             aaGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCg                                                                                 <  1:537207/151‑1 (MQ=255)
                                                              agagATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGt                                                                                <  1:277957/151‑1 (MQ=255)
                                                              agagATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGt                                                                                <  2:293796/151‑1 (MQ=255)
                                                              agagATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGt                                                                                <  2:775100/151‑1 (MQ=255)
                                                                 gaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAAAACACTGAAAAACGGGACGGCGAAGCGGGGTTTCTCCGTTACCGGCTT‑TTCTTCACCATACGACCGgtt                                                                              >  1:103728/1‑150 (MQ=255)
                                                                  aTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCCCAATTAAAAACGGGCCGGCAACGCGGGGTCTCTCCGTTT‑CGGCTTATCCTGCACAATACGACCGttt                                                                              >  2:948613/1‑149 (MQ=255)
                                                                    gAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGt                                                                          <  2:883859/151‑1 (MQ=255)
                                                                        aCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACg                                                                      <  2:78559/151‑1 (MQ=255)
                                                                          gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAg                                                                    <  2:74123/151‑1 (MQ=255)
                                                                          gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGACAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAg                                                                    <  1:794110/151‑1 (MQ=255)
                                                                               cagcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGAt                                                               <  2:145779/151‑1 (MQ=255)
                                                                                        aGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACAATGAACAACGGGCCGGCCACGCGGGGTCTCTCCGTTA‑CGGCTTTTTCCGCACACTACAACCGCTTTCGGAAAAAggggggttagttgg                                                       >  1:736424/1‑140 (MQ=255)
                                                                                         gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGgcg                                                     <  1:388028/151‑1 (MQ=255)
                                                                                         gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGgcg                                                     <  2:841254/151‑1 (MQ=255)
                                                                                           aaaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCACATTTAAAAACGGACCGGCAACGCGGGGTCTCTCCCTTA‑CGGGCTTTTCTGCACATTACGACCGCTTTCGTAAAGAGGGGGGTATGTTgggg                                                     >  1:597314/1‑147 (MQ=255)
                                                                                            aaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGAAGCACATTTAACAGCGGACCGGGAACGCGGGTTCTCTCCCTTA‑CGGCTTATTCTGCACATTACGACCCTTTGCGTAAAGAGGGGGTTACGTTGgcg                                                     >  2:632466/1‑148 (MQ=255)
                                                                                              aTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAg                                                <  1:841367/151‑1 (MQ=255)
                                                                                               tGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACATTGAACAGCGGGCCGGCCACGCGGGTTCTCTCCGTTA‑CGGCTTTTCCTGCACATTACGACCCTTTGCGGAACGAGGGGGTTACGCTGGCGCTAAGc                                               >  2:285191/1‑151 (MQ=255)
                                                                                                 gtgtGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGccc                                             <  1:614214/151‑1 (MQ=255)
                                                                                                    tgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACAATTAACAGCGGGCCGGCAAAGCGGGTTTTCTCCGTTA‑CGGCTTATCCTGCACATTAAGAACCGTTGCGGAAAGAGGGGGGTACGTTTGCGCCAAAccc                                             >  2:232524/1‑148 (MQ=255)
                                                                                                     gCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGt                                          <  2:458527/150‑1 (MQ=255)
                                                                                                      cTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGGACATTGTAAAGCGGACCGGCCACGCGGAGTCTCTCCCGTT‑CGGCGTATCCTGCGCAATACTACCGGTTGCGTAAAGAGGGGGGTACGCTTGCGCTAAGacccgcg                                         >  2:161891/1‑146 (MQ=255)
                                                                                                      cTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAGCACATTTAACAGCGGACCGGCAACGCGGATTCTCTCCGTTT‑CGGCTTATCCTGCAAATTACGACCCCTTTCGGTACGAGGGGGGTAAGGTGgcggttaggcccggg                                         >  2:269972/1‑138 (MQ=255)
                                                                                                       tGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCACCTTGAACAGCGGACCGGCAACGCGGGGTCTCTCCGTTT‑CGGGCTATCCTTCACATTAAGACCGCTTTCGGAAAGAGGGGGGTACGTTGGGGCTTAGccc                                             >  2:535620/1‑145 (MQ=255)
                                                                                                            gTCGTGATCAGTCGCGGTAGTGGCGGGCGCGGGGGCAGCACATTGTACAACGGCCCGGCAACGCGGAGTGCGTCCGTTACCGGCTT‑TTCTGCACATTAAGACCGGTTTCGTCAAGCGGGGGGTTAGTTGGCGCCAAGACCGGGGCgg                                      >  1:100679/1‑147 (MQ=255)
                                                                                                            gTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTgg                                  <  2:398547/151‑1 (MQ=255)
                                                                                                              cGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGGACACTGAAAAAAGGGCCGGCCACGCGGGGTTTCTCCCGTT‑CGGCGTTTTCTTCACATTACGCCCGGTTTGGTAAAGAGGGGGGTAAGTTGGCGCTAAGCCCCGGGCgggggggg                                >  1:676909/1‑151 (MQ=255)
                                                                                                              cGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAACACATTTAACAGCGGGCCGGCAACGCGGATTCTCTCCCTTA‑CGGCTTTTTCTGCACATTACGACCCTTTGCGTAACGAGGGGGTTAAGTTTGCGCTAAACCCCGGGCGGCTgggg                                >  2:493402/1‑151 (MQ=255)
                                                                                                                 gATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGgcg                              <  2:540934/150‑1 (MQ=255)
                                                                                                                   tCAGTCGCGGTAGTGGCGGGGGGGGGGGCAACACAATGAAAAAAGGGCCGGGAACGCGGGGTCTCTCCCCTA‑CGGCGTTTTCTGCACATTACGACCGCTTGCGTAAAGAGGGGGGTACGTTGGGGCTAAGACCGGGGCGGCGgggggg                              >  1:844875/1‑146 (MQ=255)
                                                                                                                   tCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACATTGTACAGCGGGCCGGCGAAGCGGGGTCTCTCCGTTA‑CGGCGTATCCTGCCCATTTAGAACGGTTGCGTAAAGAGGGGGGTAAGTTGGGGGTAAGCCCCGGGCGGCTgggggg                              >  1:428305/1‑146 (MQ=255)
                                                                                                                    cAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAAt                          <  1:51748/151‑1 (MQ=255)
                                                                                                                                   cgggcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGCGATTACGTTGGCGCTAAGCCCGGTACGGCTGGGGCGCAATCCTCATCTTGCAgct           <  1:272299/151‑3 (MQ=255)
                                                                                                                                    gggcgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTa          <  2:214763/151‑1 (MQ=255)
                                                                                                                                       cgAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTa       <  1:814979/151‑1 (MQ=255)
                                                                                                                                         aGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTaa      >  1:588975/1‑150 (MQ=255)
                                                                                                                                            gTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAt  <  2:730834/151‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                      
CGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTA‑CGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAACAT  >  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: ‑