Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913298,9600AC51.2% 42.2 / 70.4 41Y659D (TAT→GAT) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (20/0);  new base C (0/21);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 3.72e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-08

TATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACAC  >  NC_000913/298813‑299105
                                                                                                                                                    |                                                                                                                                                 
tataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACC‑GGCACCTCATAccc                                                                                                                                                >  2:620211/1‑151 (MQ=255)
         cgTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATc                                                                                                                                       >  2:661680/1‑151 (MQ=255)
               cGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATCGTGAtt                                                                                                                                 >  2:807177/1‑151 (MQ=255)
                 aTAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATCGTGATTa                                                                                                                                >  2:470914/1‑150 (MQ=255)
                     ccGCGCCCCCGCCGCACGGCCCCCGCTCCCCGACCCCTTTGCCCTTCATCGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTTTTTGGGGATACCCGCATGACCC‑GGCCCCTCATCCCCCGCCATATCGTGATTAACGa                                                                                                                            <  2:26938/150‑1 (MQ=255)
                      cgcgCTCCCGCCGCACAGCCCCAGCTCCCCGACCCTTTTGCCCTTCATGGGGGGGGATATGGGGTGGTTATCATCCAGGAAAACCCCCTCCTGTTTGGGAATACCCGAATGACCC‑GGCCCCTCATCCCCCGCCATATCGTGATTAACGAAg                                                                                                                          <  2:319990/151‑1 (MQ=255)
                        cgcTCACGCCGCACAGGCCCAGCCCACCGCCCCCTTTGGCCTTCATGGGGGGGGATATGGGGTCGGTATCATCCAGAAAAACCCCCTCCTGTTTTGGGATATCCGCATGAACC‑GGCCCCTCATCCCCCGCCATATCGTGATTAAc                                                                                                                              <  2:284513/145‑1 (MQ=255)
                          ccaccgcgcaacgggCCCACCCCCCCCACCCTTTTGGCTTTCATGGGGGAGATTAGGGGGTCGTTATCACCCAGAAAAACCCCCTCCTTTTTTGGGTTACCCGCATGACCC‑GGCACCTCATCCCCCGCCATATCGTGATTAACGAAGt                                                                                                                         <  2:599647/135‑1 (MQ=255)
                                ccgcACAGGCCCAGCTCCCCGACCCCTTTGGCCTTCAGCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAATCCCCTCCTGTTTTGGGATACCCGCATGACCC‑GGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                    <  1:213416/147‑1 (MQ=255)
                                ccgcACAGCCCCACCCCCCCGACCCCTTTGCCCTTCTTCGGGGAGGATAGGGGGTCGTTATCATCCAGAAAAACCCCCTCCTGTTTTGGAATACCCGCATGACCC‑GCCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc                                                                                                                <  2:302333/151‑1 (MQ=255)
                                  gcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACgg                                                                                                              >  1:153361/1‑151 (MQ=255)
                                   cacaGGCCCAGCTCACCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATGATCCAGAAAAATCCCCTCCTTTTTGGGGATATCCGCATGAACC‑GCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGt                                                                                                             <  1:363183/151‑1 (MQ=255)
                                     caGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc                                                                                                            >  2:177869/1‑150 (MQ=255)
                                               ccACCGACCCCTTTGGCCTTCATGGGGGGGGATATGGGGTCGGTATCATCCAGAAAAATCCCCTCCTTTTTTGGGATATCCGAATGAACC‑GCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa                                                                                                 <  1:520774/150‑1 (MQ=255)
                                                                       ggggAGGATGGGGGGTGGGTACCCCCCAGAAAAACCCCCCCTTTTTTGGGGATCCCCGCAGGACCC‑GCCACCCAACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccca                                                                         <  2:362799/151‑1 (MQ=255)
                                                                       ggggAGGATAGGGGGGCGTAATCACCCAGAAAAACCCCCCCTTGTTTGGGAATCCCCGCATGACCC‑GCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccc                                                                          <  2:514794/150‑1 (MQ=255)
                                                                              atatGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                                  >  1:808466/1‑151 (MQ=255)
                                                                                     gTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                              >  1:409763/1‑148 (MQ=255)
                                                                                             cacccAGAAAAACCCCCCCCTTTTTGGGAATACCCGCATGACCC‑GCCCCCCCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                              <  2:409721/137‑1 (MQ=255)
                                                                                              atcCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg                                                  >  1:817796/1‑151 (MQ=255)
                                                                                                       aaTCACCTCCTGTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg                                         >  1:203159/1‑151 (MQ=255)
                                                                                                        aTCACCTCCTGTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgc                                        >  2:55543/1‑151 (MQ=255)
                                                                                                          cccctcctTTTTGGGGATACCCGCATGACCC‑GGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCACCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg                                       <  1:161364/148‑1 (MQ=255)
                                                                                                            cctcctGTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATCGTGATTACCGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGACAATGACCTGGCTgcgcg                                     >  2:469363/1‑150 (MQ=255)
                                                                                                            ccccctGTGTGGGAATATCCGCAGGAACG‑GCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                                    <  1:116399/148‑1 (MQ=255)
                                                                                                               cttGTTTTGGGATACCCGCATGACCC‑GGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAg                                  <  1:355406/148‑1 (MQ=255)
                                                                                                                  gTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt                              >  2:24063/1‑151 (MQ=255)
                                                                                                                  gTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt                              >  1:459041/1‑151 (MQ=255)
                                                                                                                           ataccccaagAACCGGGCCCCACATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc                       <  1:381654/141‑1 (MQ=255)
                                                                                                                             cccccggcaaCC‑GGCACCTAACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAg                     <  1:284812/141‑1 (MQ=255)
                                                                                                                             ccccaaTGACCC‑GCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc                       <  2:136752/142‑1 (MQ=255)
                                                                                                                              ccGCATGACCC‑GCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGAt                  <  1:572636/151‑1 (MQ=255)
                                                                                                                                gCATGAACC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc                >  2:262308/1‑151 (MQ=255)
                                                                                                                                  aTGACCC‑GGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCg               <  1:305636/150‑1 (MQ=255)
                                                                                                                                  aTGAACC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGt              >  2:792397/1‑151 (MQ=255)
                                                                                                                                   tGAACC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTc             >  1:432533/1‑151 (MQ=255)
                                                                                                                                   aggacccgcccccTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg                 <  2:458998/137‑1 (MQ=255)
                                                                                                                                     aaCC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt           >  1:433006/1‑151 (MQ=255)
                                                                                                                                       cc‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc         >  1:85547/1‑151 (MQ=255)
                                                                                                                                        ccgcccccccATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc         <  2:114107/142‑1 (MQ=255)
                                                                                                                                               cTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacac  >  2:343717/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACC‑GGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACAC  >  NC_000913/298813‑299105

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: ‑