Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,336,901	0	T	A	43.3%	29.6 / 13.5	30	Y44F (TAT→TTT)	adiC	arginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (13/4);  new base A (0/13);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 2.17e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

TAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGC  >  NC_000913/4336754‑4337051                                                                                                                                                    |                                                                                                                                                       tAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATaa                                                                                                                                                      >  1:113468/1‑150 (MQ=255)    agGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATgg                                                                                                                                                  <  1:61022/151‑1 (MQ=255)      gACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGGCGAGGAAAGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGc                                                                                                                                                 >  1:427625/1‑150 (MQ=255)        cGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCaa                                                                                                                                               >  2:225057/1‑150 (MQ=255)        cGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAAt                                                                                                                                              >  2:679310/1‑151 (MQ=255)                       cGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAgg                                                                                                                                >  1:23965/1‑150 (MQ=255)                           aaaCGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAgg                                                                                                                           >  2:708906/1‑151 (MQ=255)                                    aaaGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGt                                                                                                                  >  1:501049/1‑151 (MQ=255)                                     aaGCGGGGGCGGGCGTAAGCGAAAAACCCCCCAGGATTTGGGTCGAGAAACGACATTTTGGCTTACCCCAGGGAGGGCCCCGGTGCACCGAAAATCGTCACCACCCATCCAAAAGGGACATTCCCGCCAGTAGAGGCCAGGTTTGCAGGt                                                                                                                  <  1:470641/150‑1 (MQ=255)                                                     aaacataaaaccccccaGGACTTGGGTCGAGGAACGACATTTTGGCGTATCCCATCGAGAACCCCTGTCCCCCAATAATCGTCACCAACCATCAAAAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACAcc                                                                                                    <  1:577653/134‑1 (MQ=255)                                                         gtaagAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAcccc                                                                                             >  2:700970/1‑151 (MQ=255)                                                           aaaaccccccaGGATTTGGGTAGAGAAACGACTTTTTGGTGTAACCCACCGAGCGCCCCAGTGCACAGAAAACCATCACCAACCATCAAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGAcccc                                                                                             <  1:400588/141‑1 (MQ=255)                                                           aaaaaCCCCCAGGACTGGGGTAGAGGAACGAAATTTGGGCGTACACAATCGAGCGCCCCGGGGCACGAATAACCGACCCAAACCACCAAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGa                                                                                                 <  1:541945/142‑1 (MQ=255)                                                              accccccaGGACTTGGGTCGAGGAACGACATTTTGGCGTATCCCATCGAGAGCCCCAGTGCACCGATAATCGCCACCAACCATCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCataat                                                                                        <  2:85849/146‑1 (MQ=255)                                                                         cTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGac                                                                              >  1:261983/1‑150 (MQ=255)                                                                          ttGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCCCCGATAATCGTCACCAACCATCCAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGaca                                                                             <  2:391074/150‑1 (MQ=255)                                                                                   ggAACGACATTTTGGCGTATCCCATCGAGAGCCCCAGGGCACCGATAATCGTCCCCACCCATCAAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAgg                                                                     <  2:262038/149‑1 (MQ=255)                                                                                        gACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgg                                                              >  1:558615/1‑151 (MQ=255)                                                                                         actttttGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCACCCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACggg                                                             <  2:238052/148‑1 (MQ=255)                                                                                          cATTTTGGCTTATACCACGGAGAGCCCCAGCGCCCCAATAATCGCCACCAACCACCAAAAAAGGCCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACgggg                                                            <  2:254004/151‑1 (MQ=255)                                                                                             tttGGCGTATACCATCGAGACCCCCAGTGCCCCGATAATCGTCCCCAACCACCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGAt                                                          <  2:289693/150‑1 (MQ=255)                                                                                                     atcccATGGAGCCCCCCAGTGACCCAATATCCGCCACCAACCACCAATAAAGGCCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCAc                                                 <  1:425324/148‑1 (MQ=255)                                                                                                       aCCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACtt                                               >  2:624571/1‑151 (MQ=255)                                                                                                                  ccccGGGGCCCCGAAAACCGTCCCCACCCACCAAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatca                                    <  1:708820/151‑1 (MQ=255)                                                                                                                   cccAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGagcatc                                     >  2:368273/1‑149 (MQ=255)                                                                                                                       gagccccGAAAATCGTAACCAACCATCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCACCGACACCATCAGGGTAACAGGAATTAAGCCCACTTTGTGagcatcagcatc                               <  2:718551/146‑1 (MQ=255)                                                                                                                            ccGAAAATCGTCCCAACCCACCCAAAAAGGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAg                          <  1:211818/151‑1 (MQ=255)                                                                                                                                 aatctgccccaaccaCCCAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATa                     <  1:632658/143‑1 (MQ=255)                                                                                                                                   tCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAt                   >  1:321232/1‑151 (MQ=255)                                                                                                                                    ctccaccaccacaccatAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAg                          <  2:537218/130‑1 (MQ=255)                                                                                                                                         ccaacaaACAATAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTa                 <  1:531228/147‑1 (MQ=255)                                                                                                                                                    tAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGc  >  2:459590/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                       TAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGC  >  NC_000913/4336754‑4337051

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