Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_000913	867,822	0	T→C	16.2%	10.5 ‑1.0 ‑0.9	37	V101V (GTT→GTC)	gsiA	glutathione ABC transporter ATPase
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands.
Reads supporting (aligned to +/- strand):  new base (1/5):  ref base (14/17):  total (15/22)
Fisher's exact test strand distribution p-value = 3.68e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.29e-01

CTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTC  >  NC_000913/867622‑868028                                                                                                                                                                                                         |                                                                                                                                                                                                               ctctctTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCg                                                                                                                                                                                                       <  2:291318/210‑1 (MQ=255)  tctctTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCgg                                                                                                                                                                                                      >  2:221632/1‑210 (MQ=255)          aGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATATCCGCACCGCGAACATGGCg                                                                                                                                                                                              <  2:436642/210‑1 (MQ=255)           gTCTGCAACGCGGTGAGAGGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATATCCGCACCGCGAACATGGCgc                                                                                                                                                                                             <  2:51984/210‑2 (MQ=255)           gTCTGAAACGCGGGGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATATCCGCACCGCGAACATGGCgc                                                                                                                                                                                             <  2:13992/210‑2 (MQ=255)             cTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCgatg                                                                                                                                                                                           <  2:424489/210‑1 (MQ=255)               gCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATATCCGCACCGCGAACATGGCgcattt                                                                                                                                                                                         <  2:169344/210‑6 (MQ=255)                 aaCGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGACATGGCGATGatat                                                                                                                                                                                       >  1:338981/1‑210 (MQ=255)                     cgGTGAGACGCTGGCAATTGTTGGCGAATCCGTCGCCGGTAAGTCAGTGGCTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGTGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCTGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGTCGATGATATTTCa                                                                                                                                                                                   <  1:177383/210‑1 (MQ=255)                             cGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGa                                                                                                                                                                           >  1:264663/1‑210 (MQ=255)                              gCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGTATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGAt                                                                                                                                                                          >  2:57954/1‑210 (MQ=255)                               cTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATg                                                                                                                                                                         <  1:277216/210‑1 (MQ=255)                                                   cGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGAGTATATTTTCGGAGCCGATCACATCGCTGAACTCGGTAtt                                                                                                                                                     >  2:367921/1‑210 (MQ=255)                                                                             tGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAAt                                                                                                                           <  1:196414/210‑1 (MQ=255)                                                                                     aTGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGt                                                                                                                   <  2:246118/210‑1 (MQ=255)                                                                                        cgcCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCTCTGAACCCGGTATTTACTGTGGGTGACCAGCTTGCCGACGCAAATCGTCTg                                                                                                                >  2:286288/1‑210 (MQ=255)                                                                                        cgcCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTg                                                                                                                <  1:9418/210‑1 (MQ=255)                                                                                        cgcCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTg                                                                                                                <  1:424121/210‑1 (MQ=255)                                                                                             gTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCa                                                                                                           <  1:242291/210‑1 (MQ=255)                                                                                                ggAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGGGAACAGATTGCCGAATCAATTCGTCTGCATCAGaa                                                                                                        >  1:82558/1‑210 (MQ=255)                                                                                                           gcggTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTg                                                                                             <  1:446893/210‑1 (MQ=255)                                                                                                            cggTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTga                                                                                            >  1:173697/1‑210 (MQ=255)                                                                                                              gTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGATCCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTgaag                                                                                          >  1:15511/1‑210 (MQ=255)                                                                                                                     aCAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACCGATTTCCGAATCAATTCGTCTGCATCAGTACGCCAGTCGTGAAGAAGCGAt                                                                                   >  2:18720/1‑210 (MQ=255)                                                                                                                        gTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGt                                                                                >  2:473347/1‑210 (MQ=255)                                                                                                                                        ttGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTCCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATg                                                                <  1:464136/210‑1 (MQ=255)                                                                                                                                         tGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTCCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAATCGAAGGTCGAGGCGAAGccgcaga                                                               >  2:464110/1‑204 (MQ=255)                                                                                                                                           cAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTg                                                             <  2:250780/210‑1 (MQ=255)                                                                                                                                            aGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTgg                                                            <  2:161103/210‑1 (MQ=255)                                                                                                                                                       cgcgAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGc                                                 >  2:120589/1‑210 (MQ=255)                                                                                                                                                                        tAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAAc                                <  1:221642/210‑1 (MQ=255)                                                                                                                                                                             aGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTc                           <  1:57956/210‑1 (MQ=255)                                                                                                                                                                                       tGCACACATGAGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTa                 <  2:312114/210‑1 (MQ=255)                                                                                                                                                                                          acaAATGCGCCATGTTAGCGGTGCGGATATGTCGATGATATTTCTGGAGCCGATGACATCGCTGAACACGGTATTTACTGTGGGTGAACAGATGGCCGAATCAATTCGTCTGCATCAGAACGCTAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATcc              <  1:126204/210‑1 (MQ=255)                                                                                                                                                                                             aaTGCGCCATGTTCGCGGAGCGGATATGACGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCa           <  1:61711/210‑1 (MQ=255)                                                                                                                                                                                                    cATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTCTTCGCATCAActc    >  2:469796/1‑210 (MQ=255)                                                                                                                                                                                                      tGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGACGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATCCTTTCACGTTACCCGCATCAActctc  >  2:104747/1‑210 (MQ=255)                                                                                                                                                                                                         |                                                                                                                                                                                                               CTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTC  >  NC_000913/867622‑868028

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