Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,9170AC63.0% 14.3 / 28.5 27V43G (GTG→GGG) ssuBaliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (9/1);  new base C (0/17);  total (9/18)
Fisher's exact test for biased strand distribution p-value = 2.13e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-05

GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993768‑994064
                                                                                                                                                     |                                                                                                                                                   
gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCcca                                                                                                                                                     >  1:476802/1‑150 (MQ=255)
    cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCcaccac                                                                                                                                                 <  2:80263/150‑1 (MQ=255)
                                ttgtCTTCCTGAATTCCAGCCGGCGGTGTGGTCCCCGCTAACCCACCGCCTGCGGTTGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCAGGGTTCTTTTCCCACCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                    <  2:378490/150‑1 (MQ=255)
                                    ctcccTGAATTCCGCCCGGGGGTGGGGTCCCCCCTAACACACCCCCTGGGGTGGGCTTTCCCACCCCTGCCAGCGGGCGCAGCAGGGTTCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                                <  1:10672/148‑1 (MQ=255)
                                    cTTCCTAAATTTCACCCGGGGGTGTGGTGCCCGCTAACACATCGCCTGCGTTGGGCGTTCCCAGCCCTGCCAGCAGGCGCAGCAGGGTACTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                                <  2:37290/151‑1 (MQ=255)
                                     ttCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGAGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGACACCACCGCTGCGGCCCACCACCGACACAAACTGACCTTCAGGAATCTGTAAATc                                                                                                                >  1:17228/1‑150 (MQ=255)
                                                      gcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat                                                                                              >  2:16819/1‑151 (MQ=255)
                                                             ggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTc                                                                                        >  1:216664/1‑150 (MQ=255)
                                                             ggcccccGCAAAACCACCGCCGGCGGTTGGCTTTCCCAGCCCTGCCGGCGGGCGAAGCGGGGTCTTTTCCCCCCCCCCGCTGGGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTc                                                                                        <  1:588873/146‑1 (MQ=255)
                                                                      aaaCACCCCGCGTGCGGTTGGCGTTCCCAGCCCTCCCAGAGGCCGCAGCAGGGTCCTTTCCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                              <  1:399391/150‑1 (MQ=255)
                                                                                 cggCGTTTGGCTTTCCCAGCCCTGCCAGCGGGCGCAGCGGGTTACTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                   <  2:728586/149‑1 (MQ=255)
                                                                                    cGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAtt                                                                >  2:638637/1‑151 (MQ=255)
                                                                                             tttCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                       >  1:597328/1‑151 (MQ=255)
                                                                                             tttCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                       >  2:392339/1‑151 (MQ=255)
                                                                                             ttcccaccccTGCCAGCAGGCGCAGCAGGGTATTTTGCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                       <  2:388405/143‑1 (MQ=255)
                                                                                             ttcccAGCCCTGCCAGCAGGCGCAGCAGGGTCCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                       <  2:718361/148‑1 (MQ=255)
                                                                                             ttcccAGCCCTCCCAGCAGGCGCAGGGGGGTTCTTTTCCCCCCCCCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                       <  2:140795/148‑1 (MQ=255)
                                                                                                    cccTGCCAGCAGGCGCAGCAGGTTATTTTTGCCCCCACCGCTGCGCCCCCCCACCGCCGCAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                  <  2:288953/148‑1 (MQ=255)
                                                                                                     ccTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                                >  2:422583/1‑150 (MQ=255)
                                                                                                        gCCAGCGGGCGCAGCGGGTTACTTTGCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGccc                                             <  2:309485/150‑1 (MQ=255)
                                                                                                           agcagGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                         >  1:691566/1‑151 (MQ=255)
                                                                                                                   cagcgggttctttttGCCCCCCCCGCTGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                  <  1:298299/139‑1 (MQ=255)
                                                                                                                   cagcagGTTACTTTCCCCACCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                  <  1:583382/150‑1 (MQ=255)
                                                                                                                          gTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                          >  1:588555/1‑151 (MQ=255)
                                                                                                                                tgcccccccccGCTGCGGCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                    <  1:118755/142‑1 (MQ=255)
                                                                                                                                tccccccccccGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                    <  1:148029/142‑1 (MQ=255)
                                                                                                                                tccccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                    <  1:8597/142‑1 (MQ=255)
                                                                                                                                  cccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg                  <  2:566640/144‑1 (MQ=255)
                                                                                                                                  cccaccaCCGGTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg                  <  1:584723/150‑1 (MQ=255)
                                                                                                                                                  ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATg  <  1:70123/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993768‑994064

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