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

CGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACAC  >  NC_000913/298820‑299105
                                                                                                                                            |                                                                                                                                                 
cgcgTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCat                                                                                                                                          >  1:241975/1‑150 (MQ=255)
              ccGCGCTCACGCCGCACAGGCCCAGCTCACCGACCCCTTTGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGAAAAATCCCCTCCTGTTTGGGGATACCCGCAGGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGaa                                                                                                                           <  1:483143/151‑1 (MQ=255)
                       cgccgcGCAGGCCCAGCCCCCTCACCCTTTTGGCCTTCATGGGGGAGGATATGGGGTGGTTATCACCCAGAAAAACCCCCCCCTTTTTGGGAATACCCGCAGGACCCGACACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCaa                                                                                                                  <  2:523119/151‑1 (MQ=255)
                       cgccgcACAGGCCCAGCTCACCGACCCCTTTGCCCTTCATCGGGGAGGATATGGGGTGGGTATCATCCAGGAAAATCCCCTCCTGTTTTGGGATATCCGCATGAACCGGCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCaa                                                                                                                  <  1:491790/151‑1 (MQ=255)
                           gcACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACgg                                                                                                              >  1:506381/1‑151 (MQ=255)
                               agccccAGCTCCCCGCCCCCTTTGGCCTTCATGGGGGGGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTTTTTTGGGATATCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCa                                                                                                           <  2:505541/147‑1 (MQ=255)
                                ggCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAATCCCCTCCTGTTTTGGGATATCCGCATGACCCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCa                                                                                                           <  1:70627/149‑1 (MQ=255)
                                                          ttCATCGGGGAGGATAGGGGGTCGGTACCACCCAGAAAAACCCCCCCCTGTTTGGGAATACCCGCATAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTg                                                                                  <  1:478270/148‑1 (MQ=255)
                                                                ggggAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccc                                                                          >  1:373440/1‑150 (MQ=255)
                                                                ggggAGAAGAGGGGTTGGTTAACACCCAGAAAAACCCCCCCCTTTTGGGGGATCCCCGCATGACCCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccca                                                                         <  2:379957/151‑1 (MQ=255)
                                                                     ggaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGCCCATCCACCCCCACCCCCCCCATCAGCGCCCCccccatcccc                                                                    >  1:364577/1‑151 (MQ=255)
                                                                                ggttaccacccAGAAAAACCCCCCCCTGTTTGGGGATATCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTg                                                          <  1:82395/141‑1 (MQ=255)
                                                                                         ccAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc                                                >  1:704494/1‑151 (MQ=255)
                                                                                                 aTCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgc                                        >  1:224168/1‑151 (MQ=255)
                                                                                                      ctcctGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                                    >  2:344680/1‑150 (MQ=255)
                                                                                                           gTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt                              >  1:133983/1‑151 (MQ=255)
                                                                                                                gggTTACCCGCATGACCCGCCCCCCCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGa                          <  2:544533/150‑1 (MQ=255)
                                                                                                                gggATATCCGCATGAACCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                           <  1:191206/149‑1 (MQ=255)
                                                                                                                gggATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCgg                           <  1:387174/149‑1 (MQ=255)
                                                                                                                   atatCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTc                       >  2:35948/1‑150 (MQ=255)
                                                                                                                      tCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg                    >  1:555047/1‑150 (MQ=255)
                                                                                                                            tGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTc             >  1:37867/1‑151 (MQ=255)
                                                                                                                                cgggCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCg               <  2:484803/143‑1 (MQ=255)
                                                                                                                                ccGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt           <  1:263219/149‑1 (MQ=255)
                                                                                                                                  ggCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgcag       <  1:186948/151‑1 (MQ=255)
                                                                                                                                        tCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacac  >  2:606027/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
CGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACAC  >  NC_000913/298820‑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: ‑