Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	3,781,120	0	T→C	28.6%	15.0 ‑7.5 0.0	36	intergenic (+103/‑95)	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 (0/10):  ref base (16/9):  total (16/20)
Fisher's exact test strand distribution p-value = 5.47e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

GGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGA  >  NC_000913/3780913‑3781322                                                                                                                                                                                                                |                                                                                                                                                                                                                ggCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTcc                                                                                                                                                                                                               <  1:491728/210‑1 (MQ=255)   cgcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCt                                                                                                                                                                                                             >  2:562564/1‑210 (MQ=255)            gATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTcc                                                                                                                                                                                                    >  2:575169/1‑210 (MQ=255)            caccagcgaatcttgcgtaatttCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTcc                                                                                                                                                                                                    >  2:489613/22‑210 (MQ=255)             aTCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTccccggag                                                                                                                                                                                                   >  1:273728/1‑206 (MQ=255)               cAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑ggg                                                                                                                                                                                                <  1:173679/210‑1 (MQ=255)                     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)                                                                                                                                                                                                                |                                                                                                                                                                                                                GGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCC‑GGGGAGAGGGCC‑GGGGTGA‑GGGG‑A‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCGCTAATACCGGCTTTCGTCTGCATATCATCGAACCGATGGGA  >  NC_000913/3780913‑3781322

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