Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913300,7110GC58.3% 5.1 / 27.9 24P75R (CCG→CGG) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base G (3/7);  new base C (14/0);  total (17/7)
Fisher's exact test for biased strand distribution p-value = 3.47e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-03

CGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAG  >  NC_000913/300566‑300843
                                                                                                                                                 |                                                                                                                                    
cgcCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCCGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTACTGGCAGCCCATACCCCCCgcgc                                                                                                                                  >  1:549052/1‑150 (MQ=255)
     gTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCtttttt                                                                                                                            >  1:96743/1‑150 (MQ=255)
     gTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGACAGCCCGTACCCCCCGGGCcttttt                                                                                                                            >  1:2055/1‑150 (MQ=255)
             gTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATTACAGCCCGTACCCCCCGCGCTTTTTTTgcggcgg                                                                                                                    >  1:116105/1‑150 (MQ=255)
                    ccACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATc                                                                                                             <  1:130134/151‑1 (MQ=255)
                                cAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAGCCAGGACGCCCGGCGCCTTTTTCGCGGCGGCCGGATCAAgggggggggg                                                                                                 >  2:172853/1‑147 (MQ=255)
                                    ggCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGTCCGTATCAAGGGCGGGGGGGCgg                                                                                             >  2:112531/1‑150 (MQ=255)
                                    ggCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGAAAGCCCGTACGCCCCGCGCTTTTTTCGCGGGGGCCGGATCAAgggggggggggggg                                                                                             >  2:188393/1‑143 (MQ=255)
                                        tgatgaTAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGGAATGAAAGCCAGTACGCCCGGGGCCTTTTTCGCGGCGTCCGGTTCAAAGGGGGGGGGGCGGcc                                                                                           >  1:126854/1‑149 (MQ=255)
                                                  gCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGGACGCCCCGCGCTTTTTTCGCGGCGGCCGTATCCAGGGGGGGGAGGGGTCCTTTGGCAATggg                                                                               >  2:537792/1‑150 (MQ=255)
                                                     cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGGCCGGATCAAGGGCGGGGAGGCGTCCTTTGGCAATGGCggg                                                                            >  2:253546/1‑150 (MQ=255)
                                                     cAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGGGAGGCGGCCTTTTGCAATGGCggg                                                                            >  2:551068/1‑150 (MQ=255)
                                                        tAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAAcc                                                                         <  1:240870/151‑1 (MQ=255)
                                                                  gccTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGGACGCCCCGCGCCTTTTTCGCGGCGTCCGTATCAAGGGCGGGGGGGCGTCCCTTGGGCATGGCGGGACCGGCGATATAg                                                               >  1:431124/1‑151 (MQ=255)
                                                                                  gtgtTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGGGAGGCGTCCTTTGGCAATGGCGGAAACGACGATATAGCCATAGGCGGCGTTgg                                               >  2:343494/1‑151 (MQ=255)
                                                                                              ccTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAACCAGTACGCCCCGCGCCTTTTGCGCGGCGCCCGTATCAAGGGCGGGGAGGCGGCCCTTGGCCATGGCGGGACCGACGATATAGCCATAGGGGGCGTTGGGGGcttcttt                                     >  1:493790/1‑148 (MQ=255)
                                                                                                tttGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGGGTCCGTCTCAAGGGCGGGGAGGCGCCCTTTTGCCATGGCGGAACCCGCGCTATAACCATAGGGGGCGGTGGGGGGTTCTTCATg                                  >  1:485430/1‑150 (MQ=255)
                                                                                                       aGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGGCCGTATCAAGGGGGGGGAGGCGGCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGGGTTGGGGGGCTCTTCCTGcccctcgg                          >  1:364256/1‑145 (MQ=255)
                                                                                                                gCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAg                  <  1:253542/150‑1 (MQ=255)
                                                                                                                gCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAg                  <  1:141708/150‑1 (MQ=255)
                                                                                                                     aCTGGCGGTAATGACAGCCAGTACGCCCCGCGCTTTTTTCGCGGCGTCCGTATCAAGGGCGGGGAGGCGTCCTTTTGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGGTTCTTTCTGCCATTCGTAGGCGTAGcgtgc             >  1:457572/1‑150 (MQ=255)
                                                                                                                      cTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGcgtgccg           <  2:2055/151‑1 (MQ=255)
                                                                                                                               aTGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAg  <  1:551058/151‑1 (MQ=255)
                                                                                                                               aTGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAg  <  2:364264/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                    
CGCCTGTTCGAAGGTCTCGGCCACTACCAGCGCAATGGCCTGATGATAGTGCTCAATAGTGGGGCCGCCTAACAGCCTGGCGGTGTTTTTGTCGCCTTTGCCGAGTGCCCCGGCGTTACTGGCGGTAATGACAGCCAGTACGCCCGGCGCTTTTTGCGCGGCGTCCGTATCAAGGGCGGTGAGGCGTCCTTTGGCAATGGCGGAACCGACGATATAGCCATAGGCGGCGTTGGGGGCTTCTTCATGCCATTCGTAGGCGTAGCGTGCCGTGCCGGTAG  >  NC_000913/300566‑300843

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