Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,337,7230AC58.3% 3.6 / 11.3 13G566G (GGT→GGGgyrADNA gyrase (type II topoisomerase), subunit A
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base C (0/7);  total (4/9)
Fisher's exact test for biased strand distribution p-value = 1.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-02

CGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCG  >  NC_000913/2337628‑2337846
                                                                                               |                                                                                                                           
cgACCGCGGTTGGAAAAGCAAAGAATAGGGTGGGGTGGTTTCCCCCCCAGCGTTGGGTCGAAAAAGTTTTCTTTTTTATTAGGTGGGGAATATTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGt                                                                      <  2:68399/151‑1 (MQ=255)
           ggAGAAGCACAGAATATGGTCGTGAGTGTCCGCCCCCAGCAGTCGGTCGAAAAAGTCTTCTTTTTAAATACTGGCGGAGAATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTg                                                            <  1:340591/150‑1 (MQ=255)
                 gCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGACTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGagagag                                                      >  2:235764/1‑150 (MQ=255)
                      gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                 >  1:124581/1‑150 (MQ=255)
                      gAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCa                                                 >  1:136244/1‑150 (MQ=255)
                         tatGGTCGGGAGTGTTCGCCACCAGCAGTCGGTCGATAAGGTCTTCTTCTTAAATACGGGCGGAAGTTTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGAc                                             <  1:156478/151‑1 (MQ=255)
                                    gttttCGCCCCCAGCAGTCGGTGAAAAAAGTCTTTTTTTTTAATACGGGCGGAAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTcc                                   <  2:359242/147‑1 (MQ=255)
                                     tgtTCGCCCCCAGCGGGCGGGCGATAGAGTCTTCTTTTTTAATGGGTGCGGAAGATTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCt                                  <  1:48222/150‑1 (MQ=255)
                                     tgtCCGCCACCGGCGGTCGGTGAATAAAGTTTTTTTTTTTAATGCGTGCGGCAGTTTTCCTTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCt                                  <  1:38386/150‑1 (MQ=255)
                                       tccGCCCCCGGCGGTGGGTCGAAAAGGTTTTCTTCTTTAATGGGTGCGGGAGATTTTCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTggg                               <  2:172211/149‑1 (MQ=255)
                                                  ctgtcggtcgATAAAGTCTTCTTCTTTAAGCCGGGCGGAAGTTTTCCCTTCCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGa                       <  1:258342/146‑1 (MQ=255)
                                                              aaaGTTTTTTTTTTTAGTGAGTGCGGCATATTTCCCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATg        <  1:294281/151‑1 (MQ=255)
                                                                        tcttTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTgcg  >  2:333963/1‑147 (MQ=255)
                                                                                               |                                                                                                                           
CGACCACGGCTGGAGAAGCACAGAATATGGTCGTGAGTGTTCGCCACCAGCAGTCGGTCGATAAAGTCTTCTTCTTTAATACGTGCGGCAGATTTACCTTTCCCGCCACGACGCTGCGCTTCGTATTCAGAAAGCGGCTGATACTTAACGTAGCCCTGGTGAGAGAGCGTCACGACCACATCTTCCTGGGTGATCAGATCTTCCAGGTTGATGTCTGCG  >  NC_000913/2337628‑2337846

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