Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,775,4670AC27.9% 244.9 / 13.6 111intergenic (+23/‑120)aroD/ydiF3‑dehydroquinate dehydratase/putative acetate‑CoA transferase
Reads supporting (aligned to +/- strand):  ref base A (60/20);  new base C (1/30);  total (61/50)
Fisher's exact test for biased strand distribution p-value = 9.70e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 130 total aligned reads displayed.

GTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATA  >  NC_000913/1775418‑1775516
                                                 |                                                 
gTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGGaca                                                   >  1:4676216/1‑50 (MQ=255)
 tATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAacac                                                  >  1:107351/1‑50 (MQ=255)
 tATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAacac                                                  >  1:8145740/1‑50 (MQ=255)
 tATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAacac                                                  >  1:3401685/1‑50 (MQ=255)
 tATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAacac                                                  <  1:6804198/50‑1 (MQ=255)
 tATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAacac                                                  >  1:2084442/1‑50 (MQ=255)
  aTTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACAcc                                                 >  1:5858243/1‑50 (MQ=255)
   ttAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACAccc                                                <  1:1225150/50‑1 (MQ=255)
   ttAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACAccc                                                >  1:798985/1‑50 (MQ=255)
   ttAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACAccc                                                >  1:6834325/1‑50 (MQ=255)
     aaCTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACAccccc                                              >  1:864901/1‑48 (MQ=255)
     aaCTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACAccccc                                              >  1:5642944/1‑48 (MQ=255)
      aCTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGCACACCCTcc                                             >  1:6311017/1‑50 (MQ=255)
      aCTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTcc                                             >  1:4416479/1‑50 (MQ=255)
      aCTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTcc                                             >  1:11975/1‑50 (MQ=255)
      aCTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTcc                                             >  1:1366161/1‑50 (MQ=255)
       cTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACcccccc                                            >  1:542176/1‑50 (MQ=255)
       cTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTccc                                            <  1:6207532/50‑1 (MQ=255)
        tATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACCCCCTcccc                                           >  1:5898082/1‑50 (MQ=255)
        tATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTcccc                                           >  1:5219757/1‑50 (MQ=255)
         aTTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCccg                                          >  1:2214210/1‑50 (MQ=255)
         aTTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCccg                                          >  1:7153740/1‑50 (MQ=255)
         aTTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAAAACCCACCccg                                          >  1:7605624/1‑50 (MQ=255)
          ttttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTcccc                                           <  1:3118278/48‑1 (MQ=255)
          ttttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTcccc                                           <  1:4573448/48‑1 (MQ=255)
          ttttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCccgc                                         >  1:7253370/1‑50 (MQ=255)
           tttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACcccccccccc                                        >  1:3558010/1‑50 (MQ=255)
            ttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCccgccg                                       >  1:533733/1‑50 (MQ=255)
            ttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCccgcc                                        <  1:3341431/49‑1 (MQ=255)
            ttACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCCCCccgccg                                       >  1:5554918/1‑50 (MQ=255)
              acacCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGaa                                     >  1:4919835/1‑50 (MQ=255)
              acacCAGGCATAAGCAATAATATTTCGGCGGGAACACCCCCCCCGCCGaa                                     >  1:7309322/1‑50 (MQ=255)
               cacCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAAc                                    <  1:7649304/50‑1 (MQ=255)
               cacCAGGCATAAGCAATAATATTTCGGCGGGAACACCCCCCCCGCCGAAc                                    >  1:3600161/1‑50 (MQ=255)
                acCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAAc                                    >  1:6604389/1‑49 (MQ=255)
                 ccAGGCATAAGCAATAATATTTCGGCGGGAACACCCCCCCCGCCGAACt                                   >  1:7935515/1‑49 (MQ=255)
                  cAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTa                                  >  1:8262640/1‑49 (MQ=255)
                  cAGGCATAAGCAATAATATTTCGGCGGGAAAACCCTCCCCGCCGCACTa                                  >  1:4250671/1‑49 (MQ=255)
                   aGGCATAAGCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTaaa                                <  1:7664707/50‑1 (MQ=255)
                   aGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaa                                 >  1:6216065/1‑49 (MQ=255)
                   aGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaa                                 >  1:573491/1‑49 (MQ=255)
                   aGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaa                                 >  1:5109380/1‑49 (MQ=255)
                   aGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaa                                 >  1:6501999/1‑49 (MQ=255)
                      cATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaaaaaa                             <  1:5791607/50‑1 (MQ=255)
                      cATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaaaaaa                             >  1:1063937/1‑50 (MQ=255)
                      cATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaaaaaa                             >  1:4318001/1‑50 (MQ=255)
                      cATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaaaaaa                             >  1:3511590/1‑50 (MQ=255)
                       aTAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaaaaaa                             >  1:3722284/1‑49 (MQ=255)
                       aTAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTaaaaaa                             >  1:1066089/1‑49 (MQ=255)
                       aTAAGCAATAATATTTCGGCGGAAACCCCCCCCCCGCCGAACTAAAAAat                            <  1:5900013/50‑1 (MQ=255)
                        tAAGAAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAata                           <  1:8083491/50‑1 (MQ=255)
                           gCAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAatatat                        <  1:2158442/50‑1 (MQ=255)
                           gCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAatatat                        >  1:4383385/1‑50 (MQ=255)
                           gCAATAATATTTCGGCGGAAACACCCTCCCCGCCGAACTAAAAAatatat                        <  1:7408498/50‑1 (MQ=255)
                            cAATAATATTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATAtt                       <  1:5437596/50‑1 (MQ=255)
                            cAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATAtt                       >  1:4658259/1‑50 (MQ=255)
                             aataatATTTCGGCGGGAACACCCTCCCCGCCGGACTCAAAAATATAtt                       >  1:4828673/1‑49 (MQ=38)
                             aataatATTTCGGCGGGAACACCCTCCCCGCCGCACTAAAAAATATAtt                       >  1:2144153/1‑49 (MQ=255)
                             aataatATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATAtt                       >  1:7952067/1‑49 (MQ=255)
                              ataataTTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCa                     <  1:5187036/50‑1 (MQ=255)
                              aaaataTTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCa                     <  1:2287437/48‑1 (MQ=255)
                                 gtatTTCGGGGTGACCCCCCTCCCCGCCGAACTAAAAAATATATTCAATc                  <  1:5677898/49‑1 (MQ=255)
                                 atatTTCGGCGGGACCACCCCCCCCGCCGAACTAAAAAATATATTCAAt                   <  1:1041072/49‑1 (MQ=39)
                                 atatTTCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAAt                   <  1:4246173/49‑1 (MQ=255)
                                  tatTTCGGGGGGACCCCCCCCCCCGCCGAACTAAAAAATATATTCAATCg                 <  1:7861963/50‑1 (MQ=255)
                                  tatTTCGGCGGGAACCCCCCCCCCGCCGAACTAAAAAATATATTCAATCg                 <  1:369125/50‑1 (MQ=255)
                                   atTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCg                 >  1:1263319/1‑49 (MQ=255)
                                   atTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCg                 >  1:226927/1‑49 (MQ=255)
                                    tttGGGGGGAAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGt                <  1:2790364/49‑1 (MQ=37)
                                    tttCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTa               <  1:2426741/50‑1 (MQ=255)
                                    tttCGGCGGGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTa               <  1:7364252/50‑1 (MQ=255)
                                    tttCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCgga               >  1:5884601/1‑48 (MQ=255)
                                    tttCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTa               >  1:3740255/1‑50 (MQ=255)
                                    tttCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTa               >  1:1890545/1‑50 (MQ=255)
                                      ccggcggGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTAt              <  1:3525626/48‑1 (MQ=255)
                                       cggcggGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTAtt             <  1:593945/49‑1 (MQ=255)
                                       cggcggGAACCCCCCCCCCGCCGAACTAAAAAATATATTCAATCGTAtt             <  1:8170937/49‑1 (MQ=39)
                                         gcggGAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTtaa          <  1:5226427/50‑1 (MQ=255)
                                         gcggAAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTtaa          <  1:2503530/50‑1 (MQ=255)
                                         gcggAAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTtaa          <  1:4236424/50‑1 (MQ=255)
                                          cggGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTtaat         <  1:2337327/50‑1 (MQ=255)
                                           gggAACCCCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTtaata        <  1:2222793/50‑1 (MQ=255)
                                           gggAACCCCCCCCCCGCCGAACTAAAAAATATATTCAATCGTATTtaata        <  1:7211796/50‑1 (MQ=255)
                                           gggAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTtaat         >  1:1645277/1‑49 (MQ=255)
                                           gggAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTtaat         >  1:6712609/1‑49 (MQ=255)
                                           gggAACACCCCCCCCGCCGAACTAAAAAATATATTCAATCGTATTtaata        <  1:6644501/50‑1 (MQ=255)
                                           gggAACACCCCCCCCGCCGAACTAAAAAATATATTCAATCGTATTtaata        <  1:6402528/50‑1 (MQ=255)
                                            ggAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTtaata        >  1:2601927/1‑49 (MQ=255)
                                             gAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATaaa      <  1:4683912/50‑1 (MQ=255)
                                              aaCACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATaaa      >  1:6655270/1‑49 (MQ=255)
                                              aaCACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATaaa      >  1:5884874/1‑49 (MQ=255)
                                              aaCACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATaaa      >  1:3919981/1‑49 (MQ=255)
                                              aaCACCCCCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATaaaa     <  1:2108894/50‑1 (MQ=255)
                                                cacCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATaaaaa    >  1:396825/1‑49 (MQ=255)
                                                 acCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAata  <  1:6233736/50‑1 (MQ=255)
                                                 acCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAata  <  1:2162139/50‑1 (MQ=255)
                                                 acCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAata  >  1:3974455/1‑50 (MQ=255)
                                                 acCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAata  >  1:7173027/1‑50 (MQ=255)
                                                 acCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAata  >  1:8289517/1‑50 (MQ=255)
                                                 acCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAata  <  1:7108093/50‑1 (MQ=255)
                                                 |                                                 
GTATTAACTATTTTACACCAGGCATAAGCAATAATATTTCGGCGGGAACACCCTCCCCGCCGAACTAAAAAATATATTCAATCGTATTTAATAAAAATA  >  NC_000913/1775418‑1775516

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