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

CGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAAC  >  NC_000913/1153432‑1153720
                                                                                                                                             |                                                                                                                                                            
cGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAGAGcaaa                                                                                                                                                      >  1:381540/1‑148 (MQ=255)
 gACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACAtt                                                                                                                                                   <  2:202796/151‑1 (MQ=255)
  aCTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTg                                                                                                                                                  <  1:294933/151‑1 (MQ=255)
     aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAAAGcacaat                                                                                                                                                   >  2:388867/1‑145 (MQ=255)
     aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGaa                                                                                                                                                <  1:128707/150‑1 (MQ=255)
     aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAAc                                                                                                                                               <  2:235733/151‑1 (MQ=255)
     aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAAc                                                                                                                                               <  1:230480/151‑1 (MQ=255)
      cAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACa                                                                                                                                              <  2:1008906/151‑1 (MQ=255)
         gATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGAAGcacaattaa                                                                                                                                                >  2:289864/1‑141 (MQ=255)
          aTGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAACACAATGAACa                                                                                                                                              >  1:856080/1‑147 (MQ=255)
            gCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGa                                                                                                                                        <  1:983341/150‑1 (MQ=255)
            gCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGAc                                                                                                                                       <  2:172237/151‑1 (MQ=255)
              ttGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCACATAGaaa                                                                                                                                               >  1:595005/1‑141 (MQ=255)
               tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCg                                                                                                                                     <  2:188594/150‑1 (MQ=255)
                gTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCgg                                                                                                                                    <  1:918558/150‑1 (MQ=255)
                  cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCaa                                                                                                                                 <  2:300414/151‑1 (MQ=255)
                   aGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAGGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGAGAGCACAATGAACAGGGGAGCgg                                                                                                                                    >  2:433741/1‑147 (MQ=255)
                      ggTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGAGAGGACACTGAACAGCGGAcccg                                                                                                                                    >  2:930345/1‑142 (MQ=255)
                        ttGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGAAAACACATTGAACAACGGACCGGGAAcgccg                                                                                                                            >  2:152605/1‑148 (MQ=255)
                           atgatTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATtc                                                                                                                       <  1:186462/151‑1 (MQ=255)
                                 tcatgtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGACAGACGTGATCAGACGCGGCAGTGGCGGCGGCGGGGACAGGGCAATAAACAGCGGGTCCGGCaaaccga                                                                                                                            >  1:1047937/4‑137 (MQ=255)
                                 tCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTcc                                                                                                                 <  1:574057/150‑1 (MQ=255)
                                 tCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTcc                                                                                                                 <  1:417035/150‑1 (MQ=255)
                                        aCTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGc                                                                                                        <  1:577136/151‑1 (MQ=255)
                                                   cAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGca                                                                                             <  2:387344/150‑1 (MQ=255)
                                                   cAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGca                                                                                             <  1:208793/150‑1 (MQ=255)
                                                   cAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGc                                                                                              <  2:147641/149‑1 (MQ=255)
                                                    aGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGcaca                                                                                           <  2:507674/151‑1 (MQ=255)
                                                       ttGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAGCACAATTAACAACGGGCCGGCAACGCGGGTTTTCTCCCTTACGGCCTATTCTTcac                                                                                            >  2:994891/1‑147 (MQ=255)
                                                          aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGGAAGCACATTGAAAAGAGGGGCGGGAACGCGGGTTCTCTCCGTTTCGGGGTATTCTGcac                                                                                            >  2:245545/1‑144 (MQ=255)
                                                          aaCAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACATTGAAAAACGGACCGGCCACGCGGGGTCTCTCCCTTTCGGCTTTTCCTGcc                                                                                             >  1:614817/1‑142 (MQ=255)
                                                              agagATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGt                                                                                <  2:155311/151‑1 (MQ=255)
                                                                  aTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGc                                                                            <  1:501735/151‑1 (MQ=255)
                                                                     aaaaCGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTa                                                                        <  1:116350/151‑1 (MQ=255)
                                                                            tGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGAAGCACATTTAACAACGGACCGGCGACGCCGGGTCTTTCCCCTACGGCTTTTCCCTGCCATTAACACCGCCTGCGTGACGAggg                                                                 >  1:848931/1‑151 (MQ=255)
                                                                                agcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGGAGCACAATGAAAAACGGAGCGGCAACGCGGGGTTTCTCCGTTTCGGGCTTTTCTTCCCATTATCACCCCTTGCGGAACCAgggg                                                                >  1:143401/1‑148 (MQ=255)
                                                                                agcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGAtt                                                             <  1:889561/151‑1 (MQ=255)
                                                                                     aaCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGtt                                                        <  2:120524/151‑1 (MQ=255)
                                                                                           aaaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGGAAGACCATGAACAACGGAGCGGGCACGCGGGGTTTTTCCCTTTACGCCTTTTCCCCCACATTCGAACCCGTTGGTAAAGAGGGGGGTACGtttgggg                                                   >  1:787348/1‑145 (MQ=255)
                                                                                           aaaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCt                                                  <  1:937472/151‑1 (MQ=255)
                                                                                            aaaTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGAGGGGGGGGAGCAAATTCAACAGGGGAGCGGCCACGCGGGGTCTCCCCGTTTCGGCCTTTCCTGCCCATTACAACCGGTTGAGTTAACGGGGGGGTTCCTTTgcg                                                     >  1:502402/1‑146 (MQ=255)
                                                                                                  tgtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGACAACACATTTAAAAGCGGACCGGGAACCCGGACTTTTTCCGCTATCGCCTATTCTGCACAATTCTACCCGTTGCGGTAAGAGGGGGGTATGTTGGCGATCAGccc                                           >  2:696323/1‑150 (MQ=255)
                                                                                                        gAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAGCACATTTAAAAGCGGGCCCGCCACGCGGGGTCTCTCCCCTTCGGCTTATCCTGCACATTACGACCGTTTTCGTGACGAGGGGGTTACGTTTGGGCTAAGCCCCGGGCgg                                    >  2:658081/1‑151 (MQ=255)
                                                                                                        gAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACCTTGAAAAGCGGGCCGGCCACGCGGGTTCTCTCCCTTTCGGCTTATCCTCCACATTACGACCCCTTGCGTAACGAGGGGGTTACGTTTGCGCTAAGCCCGGGGCgg                                    >  2:614729/1‑151 (MQ=255)
                                                                                                          aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGGCAGCACATTGAAAAACGGAACGGGAAAGCGGGTTCTCTTCCGTTTGGGCTTTTCTGCCCACTACGACCGCTTTGGTTACGAGGGGGTTACGTTTgcgcctagacccggggg                                      >  2:485671/1‑135 (MQ=255)
                                                                                                              cGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTgggg                              <  2:442514/151‑1 (MQ=255)
                                                                                                               gTGATCAGTCGCGGTAGTGGCGGGGGGGGGGGCAACACCATTGAAAAAGGAGCGGCCACACGGGGTCTCTCCGTTTCGGGCTTTCCTGCACAATTCGACCGGTTGTGTAAACAGGGGGATGCGTTGGCGGTCAACCCCGGGCGGGCgggg                              >  1:944224/1‑150 (MQ=255)
                                                                                                               gTGATCAGTCGCGGTAGTGGCGGGCGGGGGGACAACACATTGAACAGAGGGCCGGGAACGCGGGTTTTCTCCCTTTCGGGTTATTCCTCACATTAAGAACCGTTTCGTAACGAGGGGGGTACGTTGGCGCTAAAGCCCggggggcgggggg                              >  2:332150/1‑140 (MQ=255)
                                                                                                               gTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACAGCACATTTAACAACGGACCGGCAACGCGGGTTCTCTCCCTTACGGCTTTTCCTGCCCATTACGACCGCTTTCGGAACGAGGGGGGTACGTTTGCGCTCAACCCCGGGGGGCGggggg                             >  2:84231/1‑150 (MQ=255)
                                                                                                               gTGATCAGTCGCGGTAGTGGCGGGCGAGGGGAAAGCACCTTTAACAGCGGGCCGGCAACACGGAGTCTCTCCCGTACGGCCTTTCCCGCACCTTAAGACCGTTTTCGGAACGAGGGGGGTACGTTGGCGCTAAACCCGGCGCGGGTggggg                             >  2:955389/1‑150 (MQ=255)
                                                                                                                     aGTCGCGGTAGTGGCGGGCGGGGGGAAAGAACATTGAACCAAGGGGCGGCCAAGCGGGGTGTCTCCGCTTACGCGTATTCCGCACATTTTCACCGTTTGTGGCACGAGGGGGGTACGTTTGGGCTTAACCCGGGGCGGGCggggggggatc                       >  1:764486/1‑144 (MQ=255)
                                                                                                                        cgcgGTAGTGGCGGGCGGGGGGACAGCACATTGAACAGCGGACCGGGAACGCGGGTTCTCTCCGTTTCGGCGTTTCCTTCACATTACCACCGTTTTCGTAACGAGGGGGATAAGTTGGCGGTAAAGCCCGTGGGGGTGGGgggccaacc                      >  2:990391/1‑144 (MQ=255)
                                                                                                                        cgcgGTAGTGGCGGGCGAGGGTACAGCACTTTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGTTGGGGCGCAATCCt                     <  1:686161/150‑1 (MQ=255)
                                                                                                                        cgcgGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTc                    <  2:345395/151‑1 (MQ=255)
                                                                                                                        cgcgGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTc                    <  1:614819/151‑1 (MQ=255)
                                                                                                                                          gggTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAAc  <  1:388859/151‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                            
CGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTACAGCACATTGAACAGCGGACCGGCAACGCGGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAAAC  >  NC_000913/1153432‑1153720

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