Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,9130GC53.3% 33.6 / 50.5 30G44G (GGC→GGG) ssuBaliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (11/3);  new base C (0/16);  total (11/19)
Fisher's exact test for biased strand distribution p-value = 6.66e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.56e-06

CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCA  >  NC_000913/993763‑994060
                                                                                                                                                       |                                                                                                                                                   
ccACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCgg                                                                                                                                                     <  2:297225/151‑1 (MQ=255)
         cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCcaccac                                                                                                                                             <  2:79979/150‑1 (MQ=255)
         cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCcaccac                                                                                                                                             <  2:655139/150‑1 (MQ=255)
                     cTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGa                                                                                                                                >  2:153771/1‑151 (MQ=255)
                       aaaaCCTCTTCCGCGTTTCTTCCTAAATTCCAGCCAGGGGTGTGGGCCCCGCAAACACACCCCCGGCGTTTGGCTTTCCCACCCCTGCCAGCAGGCCCAGCAGGTTACTTTTGCCCCCA‑CCGCTGCGGCCCCCCACCGCCACAAACTGAcc                                                                                                                              <  1:644189/150‑1 (MQ=255)
                                 tCGCTTGTCTTCCTGAATTTCAGCCGGGGGTGTGGTCCCCGCTAACACACCGCCGGCGGTTGGCGTTTCCAGCCCGGCCAGCAGGCGCAGCAGGTTACTTTTGCCACCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                    <  2:323046/151‑1 (MQ=255)
                                       gtCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTaaa                                                                                                              >  2:348172/1‑151 (MQ=255)
                                               gAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                      >  2:443257/1‑151 (MQ=255)
                                                              gtggGGTCCCCGCTACCACTCCGCCGGCGTTGGGCTTTCCCAGCCCTGCCGGCAGGCGCAGCAGGGTCTTTTTGCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                        <  2:448593/150‑1 (MQ=255)
                                                                   gtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcg                                                                                  >  2:163725/1‑151 (MQ=255)
                                                                        cGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                             >  1:237401/1‑151 (MQ=255)
                                                                                caTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                       >  1:358837/1‑149 (MQ=255)
                                                                                 accGCCTGCGGTGGGCGTTCCCAGCCCTGCCAGCAGGCGCAGCGGGTTTCTTTGCCCACCA‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                    <  2:315302/149‑1 (MQ=255)
                                                                                   cGCCTGGGTTGGGTGTTCCCAGCCCGCCAAGCAGGCGCAGAGGGTTTTTTTTGCCACCA‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                  <  1:660371/151‑1 (MQ=255)
                                                                                                  ttcccAGCCCGGCCAGCAGGCGCAGCAGGGTACTTTTCCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                   <  2:461494/148‑1 (MQ=255)
                                                                                                     ccAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                >  1:121168/1‑151 (MQ=255)
                                                                                                         cccggCCAGCGGGCGCAGCGGGGTTTTTTTCCCCCCC‑CCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGt                                             <  2:616413/146‑1 (MQ=255)
                                                                                                         cccggCCAGCGGGCGCAGCAGGGTATTTTTGCCCCCC‑CCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                              <  2:476679/145‑1 (MQ=255)
                                                                                                         cccggCCAGCAGGCGCAGAGGGGTACTTTGCCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                              <  2:38911/145‑1 (MQ=255)
                                                                                                                     gcgcAGCAGGGTACTTTTCCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                   <  2:624115/148‑1 (MQ=255)
                                                                                                                                 ttttttGCCCCCCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                     <  1:556707/149‑1 (MQ=255)
                                                                                                                                 ctttgccccacca‑CGGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                    <  1:401546/144‑1 (MQ=255)
                                                                                                                                     ttcccccca‑ccGCGCCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                <  2:441361/145‑1 (MQ=255)
                                                                                                                                     ttGCCACCC‑CCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCacac                <  2:441372/151‑1 (MQ=255)
                                                                                                                                      tccccccccccGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                <  1:624865/142‑1 (MQ=255)
                                                                                                                                      tccccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                <  1:190451/142‑1 (MQ=255)
                                                                                                                                      tGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                >  2:27424/1‑150 (MQ=255)
                                                                                                                                      tGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                >  2:44457/1‑150 (MQ=255)
                                                                                                                                          ccca‑ccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg           <  1:213838/150‑1 (MQ=255)
                                                                                                                                               ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                   >  1:21378/1‑107 (MQ=255)
                                                                                                                                                    gCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa  >  1:133077/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCA  >  NC_000913/993763‑994060

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