Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,781,140	0	A→G	38.7%	12.9 ‑5.7 0.0	31	intergenic (+123/‑75)	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/11):  ref base (12/7):  total (13/18)
Fisher's exact test strand distribution p-value = 3.38e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGC  >  NC_000913/3780933‑3781349                                                                                                                                                                                                                 |                                                                                                                                                                                                                      aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCGGA‑GGGGCGAgg                                                                                                                                                                                                                     <  2:273717/210‑1 (MQ=255)  aTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAggg                                                                                                                                                                                                                    <  1:432572/210‑1 (MQ=255)  aTTACGCAAGAGTCGCTGGTGCAGGGGCTGGGCAAAGAGTTGCCTGCGGCACTGGGTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAggg                                                                                                                                                                                                                    <  2:180829/210‑1 (MQ=255)                    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)                                                                                                                                                                                                                 |                                                                                                                                                                                                                      AATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGATTTGCCTGGGACGATAAGCGCCTGCGC  >  NC_000913/3780933‑3781349

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 35 ≤ 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.