Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,536,197	0	A	T	41.9%	34.4 / 23.5	31	intergenic (+166/+417)	yjhR/nanS	pseudogene, helicase family; putative frameshift suppressor/9‑O‑acetyl N‑acetylneuraminic acid esterase
Reads supporting (aligned to +/- strand):  ref base A (1/17);  new base T (13/0);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 6.79e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGG  >  NC_000913/4536048‑4536343                                                                                                                                                      |                                                                                                                                                    aTGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAg                                                                                                                                                    <  1:349306/151‑1 (MQ=255)       cAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCGTTGCATAATGCGCCTGGGATTAATGGCGTGTGCTTTTATGAATg                                                                                                                                               <  1:70937/150‑1 (MQ=255)                           ttttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATTAAAAAAAAAtgaagcaa                                                                                                                              >  2:363298/1‑147 (MQ=255)                             ttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataata                                                                                                                         <  1:334226/150‑1 (MQ=255)                               aCCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatat‑a                                                                                                                      <  1:104469/150‑1 (MQ=255)                               aCCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAataatat‑a                                                                                                                      <  1:27980/150‑1 (MQ=255)                                 cAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAAtaatat‑aata                                                                                                                   <  2:330599/151‑1 (MQ=255)                                   aTGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAAAGAAAAAAAAAGAAAGAAATAATAt‑aataa                                                                                                                  >  1:146987/1‑150 (MQ=255)                                     gTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATgg                                                                                                               <  2:362152/151‑1 (MQ=255)                                        aaGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGcc                                                                                                             <  2:346434/150‑1 (MQ=255)                                               cTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGTTTTTTGGAATGAAAGAAAAAGCAAAAAATAAAAT‑AATAAAAGACTAATTTCt                                                                                                     >  2:81979/1‑151 (MQ=255)                                                  tCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCttt                                                                                                   <  1:7496/150‑1 (MQ=255)                                                        ttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAAcat                                                                                             <  2:120093/150‑1 (MQ=255)                                                            tGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCAtt                                                                                         <  1:83856/150‑1 (MQ=255)                                                            tGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGATTTTTTGAATTAAAAAAAAAGAAAGCAATAAAAT‑AAAAAATGCCTAATTTCTTTTAACATCAtt                                                                                         >  1:144106/1‑150 (MQ=255)                                                                aGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGCGTTTTTTTGGATTAAAAAAAAATCACGGAATAAAATAAAAAATGGCCTAATTTTTTTTAA‑ATCAtt                                                                                         >  1:17589/1‑146 (MQ=255)                                                                   cGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCatat                                                                                 <  1:64722/151‑1 (MQ=255)                                                                     ccTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTTATTGAAAGAAAAAGCAAAAAAAAATAT‑AAAAATGGTCTCATTTTTTTGAAAATCAATAAACa                                                                                    >  1:292733/1‑146 (MQ=255)                                                                      cTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTGGATGAAAAAAAAAGCAAGCAAAAATAT‑AAAAAAAGGCTTATTTTTTTTTACATAATTAAAGatat                                                                                 >  2:387835/1‑148 (MQ=255)                                                                          ccTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATatcatcat                                                                           <  2:20883/150‑1 (MQ=255)                                                                             gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGt                                                                        <  1:131403/150‑1 (MQ=255)                                                                                    aaTCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTGGAATGAAAAAAAAAACAAACAATAATAT‑AATAATGGACTAATTTTTTTTAAAAACATTAAAAATATTATTATCTTCTtttt                                                                  >  2:212564/1‑149 (MQ=255)                                                                                       cTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGGATTAAAAAAAAAGGAAACAAAAAAAA‑AAAAAAAGCCCAATTTTTTTTTAAATCATTAAAAATATAATTATCGTCTAtttt                                                                 >  2:11973/1‑146 (MQ=255)                                                                                          aTTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGTTTAAAGAAAATGCAAGCAATAATAT‑AATAAAAGCCTTATTTTTTTTAAAATAAATAAGAAAAAAATAATCGTCTATTTTCCCGaa                                                           >  1:110586/1‑150 (MQ=255)                                                                                                 aCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTc                                                     <  1:331200/149‑1 (MQ=255)                                                                                                 aCCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGAGTTTTTTTGGATGGAAAAAAAAGCCAGAAATAAAAT‑AATAATAGGCTCAATTCTTTGTAGATAATTAAACATATAATTAACTTCTAttttttaaaaatcgat                                                     >  1:61520/1‑136 (MQ=255)                                                                                                      aTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGGATGAAAGAAAAAGCAAACAATAATAT‑AATAATAGCCTAATTTCTTTTAAAATCATTAAGCATATCATCATCGTCTATTTTTCAGAATCCGTCGACTtat                                              >  2:279769/1‑151 (MQ=255)                                                                                                            ttCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGGTTTTGGGATGGAAAAAAAAGAAAGGAAAAAAAT‑AAAAAAAGGCTTATTTTTTTTTAAAACATTAAAAATATAATCATCGTCTATTTTTCCGAAACCGTCGACTtat                                              >  2:220535/1‑145 (MQ=255)                                                                                                                        aTGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCata                            <  2:132589/151‑1 (MQ=255)                                                                                                                                     tAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATGTGCCAGAATCCGTCGACTTATAATCGCCAGAACCGCATATCAGGTCGGtgc                <  1:2799/150‑1 (MQ=255)                                                                                                                                                  ttttAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTgg  >  1:188748/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                    ATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATAT‑AATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGG  >  NC_000913/4536048‑4536343

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: ‑