Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,781,152	0	A→G	30.3%	12.5 ‑5.0 ‑0.0	33	intergenic (+135/‑63)	lldD/trmL	L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Rejected: Prediction has biased read strand distribution supporting polymorphism.
Reads supporting (aligned to +/- strand):  new base (1/9):  ref base (15/8):  total (16/17)
Fisher's exact test strand distribution p-value = 6.60e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGC  >  NC_000913/3780952‑3781359                                                                                                                                                                                                           |                                                                                                                                                                                                                   gTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAg                                                                                                                                                                                                        >  1:580413/1‑210 (MQ=255)   gccggggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGc                                                                                                                                                                                                      <  1:489638/207‑1 (MQ=255)    agggggCTGAGTAAGGAGTTGCTTGCGGCCCGGCGTCCCATGGCGAAAGGGAAGGCGCCATGGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGAGGGAACCGTTCGCCGCTGTATGTCCTCCCTCTCTCTAAAACGACCCCGCACTTTTTTTTCCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGcc                                                                                                                                                                                                     <  1:110309/208‑1 (MQ=255)       gggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCa                                                                                                                                                                                                  >  2:44707/1‑210 (MQ=255)        ggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCAc                                                                                                                                                                                                 >  2:514902/1‑210 (MQ=255)                 aggggttcccggcgcacgtGGCCCCCAGGGCGAAGGGGAATGGGGCATGGTCGTTGGCGCCCCTACACCCACCCCCCCCCCCTCGGGAGGAGGGGCCCTTTCGCCGTTTTTGGTACCCCCCCCCGCTGAACCGAACCCGCCTTTTTTTTTTCCCCCTCCCCCCCCG‑GGGGAGGGGGCG‑GGGGTGG‑GGGG‑A‑AAAGGCCGCACTGCTCACCa                                                                                                                                                                                        <  2:63583/192‑1 (MQ=255)                          aggccgccctgccccccATGCCGAACGGAAATGGGGAATGGTCTTTAGCCCCCCCCCACCATCCCCCTCTCCCTCGGGAAGAGGGAACTGTTCGGCACTGTTTATCCCCCCCCGTTCGAAAACGACACCGCATTCTTTTTTCCTCCCTCCCCCCGCC‑GGGGAGAGGGCG‑GGGGGGG‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGc                                                                                                                                                                               <  1:567082/196‑1 (MQ=255)                                                            ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                                                                                                                                             >  2:621036/1‑210 (MQ=255)                                                                 aGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACAt                                                                                                                                        >  2:210293/1‑210 (MQ=255)                                                                  gTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                                                                                                                                       >  1:412814/1‑210 (MQ=255)                                                                                      cTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAt                                                                                                                   >  2:520679/1‑210 (MQ=255)                                                                                          cccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑‑GGGAGGGGGCC‑GGGGGGG‑GGGG‑A‑AAAGGCCCCCCTGCCCACCAATTTTTTGCTATCCTTCCCCCCCACTAAGGGGGCAGTTTGCTAAACAACGTACTTTACGAACCAGAAATTcccc                                                                                                              >  1:472350/1‑208 (MQ=255)                                                                                           ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGc                                                                                                              >  2:508928/1‑210 (MQ=255)                                                                                                      gggcggggcgACCGTCGGGTGTTGTCTGTACCCCCTCACTGGGCACCGGACCCGCCAGCTGTTGTTCGCCCCCGCCCC‑CC‑GGGGGGGGGGCG‑GGGGGGAGGGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGc                                                                                                   <  2:483231/201‑1 (MQ=255)                                                                                                      gggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGT‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTACACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGc                                                                                                   >  1:211416/1‑210 (MQ=255)                                                                                                                   gTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTc                                                                                      >  2:311932/1‑210 (MQ=255)                                                                                                                       gccgcTTTTTGTCCCCCCCCCCTCAGAAACACCCCCGCCCTTTTTTTTCCCCCCCCCCCCCCCG‑GGGGAGGGGGCG‑GGGGGGG‑GGGA‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTg                                                                                  <  1:621073/208‑1 (MQ=255)                                                                                                                        cccccaccctctcccgtcctctcaGACACGACCCCCCCCTC‑TTTTTTCCCCCCCCCCCCCCC‑GGGGGGGGGGCCGGGGGGGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTgc                                                                                 <  2:299326/190‑1 (MQ=255)                                                                                                                           tctttattccccccccacctagAAACCACCCCCCGCTTTTTTTTTCCCCCCCCCCCCCCC‑GGGGGGGGGGCCGGGGGGGG‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTgcgc                                                                               <  1:225396/189‑1 (MQ=255)                                                                                                                             ttatgtcCCCCCCCACTCTAAAACCCCCCCCCACTTTTTTTTCCCCCCCCCCCCCCCG‑GGGGAGGGGGCG‑GGGGGGG‑GGGA‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTaa                                                                            <  1:520706/209‑1 (MQ=255)                                                                                                                                      ccccccccctataaccccacCCCACTCTTTTTTTTTCCCCCGCCCCCCCGGGGGGGGGGGCC‑GGGGGGG‑GGGGAA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGc                                                                     <  1:276326/193‑1 (MQ=255)                                                                                                                                        ccccctctAACACGCCCCCCCACTTTTTTTCCCCCCCCCCCCCCCCG‑GGGGAGGGGGCG‑GGGGGGG‑GGGA‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTc                                                                 <  1:508955/206‑1 (MQ=255)                                                                                                                                            ctctGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTgcagatcggaagagcacacgtc                                                             >  1:299338/1‑190 (MQ=255)                                                                                                                                                     caccccccctcttttttttCCCCCCCCCCCCCCC‑GGGGGGGGGGCG‑GGGGTGG‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATatcat                                                    <  1:44708/198‑1 (MQ=255)                                                                                                                                                       cgcccaatgttttttttCCACCCCGCCCCCCC‑GGGGGGGGGCCC‑GGGGGGA‑GGGG‑C‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCg                                                  <  2:224203/200‑1 (MQ=255)                                                                                                                                                          cccACTCTTTTTTCCCCCCCCCCCCCCCC‑GGGGGAGGGGCC‑GGGGGGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAAc                                               <  1:562592/208‑1 (MQ=255)                                                                                                                                                                   tttttGCCCCCCCCCCCCCC‑GGGGGGAAGGTC‑GCGGTTG‑GAGG‑AGAAATGCCGGACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGa                                       <  2:295925/210‑1 (MQ=255)                                                                                                                                                                          cccTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTg                               >  2:246815/1‑210 (MQ=255)                                                                                                                                                                                              agGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTgcgc            >  2:321271/1‑210 (MQ=255)                                                                                                                                                                                               gggCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGcgcc           >  2:534425/1‑210 (MQ=255)                                                                                                                                                                                                 gCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGcgccgc         >  2:68792/1‑210 (MQ=255)                                                                                                                                                                                                        gagA‑GGCG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCgggg   <  2:559900/208‑1 (MQ=255)                                                                                                                                                                                                         ggA‑GGGA‑A‑CAAGGCCGCACTGCTCACCTATTTTCTGCTATCCTGCACCCGGACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGc  <  1:670499/209‑1 (MQ=255)                                                                                                                                                                                                           |                                                                                                                                                                                                                   GTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGCCGCGCGGGGC  >  NC_000913/3780952‑3781359

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.