Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,412,9060TA45.8% 26.2 / 10.3 24Y187F (TAT→TTT) acrSacrAB operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (9/4);  new base A (0/11);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 5.89e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAG  >  NC_000913/3412757‑3413050
                                                                                                                                                      |                                                                                                                                                
aaCTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTtat                                                                                                                                                  >  2:335587/1‑150 (MQ=255)
              gACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAAccc                                                                                                                                    <  2:109938/150‑1 (MQ=255)
               aCGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGc                                                                                                                                  <  2:455011/151‑1 (MQ=255)
                 gAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTATTTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCa                                                                                                                                <  2:72969/151‑1 (MQ=255)
                               ggAAAGAATGTAAATTTAATTCACACTTAATTCTTTTGTTTGATTAATTAATTTTTTTATTTTTTCATCTAGCATAAACATTCTTAATTCTTAATCACCCAGGGGGGGGGTTTTTTTAAAAAAATCATAACCCGCCATATTCATTAACCAg                                                                                                                  <  1:502327/151‑1 (MQ=255)
                                                 cATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTg                                                                                                >  2:452447/1‑151 (MQ=255)
                                                  atcacatTTATTTCATTCTTTTAAAGAATTATTTTCTTTTTTTTTTCATCGGAAATAAACATTTTTAATCCGTTATCGACCAGAGGGGGGTTTTTTTTAAAAAAATCATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCt                                                                                                 <  1:694137/146‑1 (MQ=255)
                                                           aaTTCTTTCTTTTGATAAATTAATTTCTTTTTTTTTTCATCTGGCATAAACATTCTTAATACTTTATCGACCAGGGCGGGGGTTTTTTTAAAAAAATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCa                                                                                       <  2:659382/150‑1 (MQ=255)
                                                                      ttaattaaataattTCCTTTATGTTTTAATCTGGCATAAACTTTTTTAATCCTTTTTCGCCCGGGGGGGGGGTTTTTTTAAAAATACCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCaataat                                                                               <  2:199783/139‑1 (MQ=255)
                                                                      ttGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAatc                                                                           >  2:553460/1‑151 (MQ=255)
                                                                         atgaATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAatcatc                                                                        >  1:75481/1‑151 (MQ=255)
                                                                         atgaATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAatcatc                                                                        >  2:395848/1‑151 (MQ=255)
                                                                            atttatttTCTTTTTTTTTTCTTTTGGAAAAAACATTTTTAATTCGTTAACGCCCAGGGGGGGGTTTTGTTTAAAAAGATCAAAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAAtcatca                                                                       <  1:425104/147‑1 (MQ=255)
                                                                               taattTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCacaac                                                                   >  1:409560/1‑150 (MQ=255)
                                                                                         tttGTTTTCATTTGCCATAAACATTATAAATTCTCTAGCACCCAGAGCGGGGTTTTTTTAAAAAAGTTAACAGCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATc                                                          <  1:265783/148‑1 (MQ=255)
                                                                                           tGTTTTCATCTGTAAAAAACATTCTAAATCCTTTATCGCCCAGGGCGGGGTTTTGTTAAAAAAGATAAAAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCg                                                         <  1:544519/149‑1 (MQ=255)
                                                                                           tGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCg                                                         >  2:689805/1‑149 (MQ=255)
                                                                                             ttttCATCGGGAAAAAACATTTTTAATACTTTATCGACCAGAGGGGGGTTTTTTTAAAAAAGATAAAACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAg                                                       <  1:439367/149‑1 (MQ=255)
                                                                                                       acaaaaaaaTTTTTAATACTTTATCGACAAGGGGGGGGTTTTTTTTATAAAAATAATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTac                                           <  2:277553/143‑1 (MQ=255)
                                                                                                             aCATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAAcc                                     >  1:552098/1‑151 (MQ=255)
                                                                                                                   ttAATTCTTTATCAACCAGGGCGGGTTTTTTTTAAAAAAAATAAAAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCt                                    <  1:276572/146‑1 (MQ=255)
                                                                                                                                   cAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAAt               >  2:405738/1‑151 (MQ=255)
                                                                                                                                       gggggggTTTTTTAAAAAAAATAATACCCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACtt           <  1:79678/149‑1 (MQ=255)
                                                                                                                                            ggCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCgg      >  2:458298/1‑151 (MQ=255)
                                                                                                                                                  tttAAAAATATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGagag  <  2:440045/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
AACTGGTTAACTGTGACGAACTGAATTTTCAGGACAGAATGTGAATTTACATGACACTTAATTCATTCGTTTGATGAATTAATTTCGTTATGTTTTCATCTGGCATGAACATTCTTAATACGTTATCGACCAGAGCGGGGGCTTGTTTATAAAGATCATAACCCGCCATATTCATTAACCAGTTTTGAACAATTCCGCTGAAGGCACCATCAATAATAATCATCACAACATCTAAATCGAGGTTATTTGCTACACAACCTTGTTGCTGACACGCCTGCAATACTTCGCGGAGAG  >  NC_000913/3412757‑3413050

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