Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,880,045	C→T	M25I (ATG→ATA)	gyrB ←	DNA gyrase subunit B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,880,045	0	C	T	96.4%	89.8 / NA	28	M25I (ATG→ATA)	gyrB	DNA gyrase subunit B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/15);  minor base A (1/0);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGC  >  NC_000913/3879799‑3880277                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         gTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATatatatac                                                                                                                                                                                                                                      <  1:321919/251‑1 (MQ=255)                   gCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCg                                                                                                                                                                                                                    <  1:308569/251‑1 (MQ=255)                   cccTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCTATATATATACCCGGGCGCTTACGCACCg                                                                                                                                                                                                                    <  2:415908/250‑1 (MQ=255)                    ccTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGc                                                                                                                                                                                                                   >  1:120216/1‑251 (MQ=255)                                         aTCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACttt                                                                                                                                                                                              <  2:112062/251‑1 (MQ=255)                                                   ccgccgATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTggag                                                                                                                                                                                    <  1:225878/251‑1 (MQ=255)                                                      ccgAGACGCCCTCTTCCGGATGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGGACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATAGACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGgaggag                                                                                                                                                                                 <  1:162213/251‑1 (MQ=255)                                                        gATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCCTGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCCTCCAGCCCTTTCAGGAATTTGATACTGGAGGAGTc                                                                                                                                                                               >  2:392087/1‑251 (MQ=255)                                                        gATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTc                                                                                                                                                                               <  1:161037/251‑1 (MQ=255)                                                           aCGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGTGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATa                                                                                                                                                                            <  1:390407/251‑1 (MQ=255)                                                                           gAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAAc                                                                                                                                                            <  1:286909/251‑1 (MQ=255)                                                                            aaTACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTc                                                                                                                                                                                                      <  2:75805/208‑1 (MQ=255)                                                                            aaTACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTc                                                                                                                                                                                                      >  1:75805/1‑208 (MQ=255)                                                                                 cGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGGATGGGGACGATAAATTCTTTACAGGGACCCGCGGGCGCTTTGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGGGCAGACCGGTGCCGTCATCCGGGTTGCCGATATATATACCCGGGCGCTTACGGACCGCATCCAGCCCTTTCAGGACATTGATACTGGGGGAAGCATAAGAATTAGAAATCAACGTTtct                                                                                                                                                      >  1:14463/1‑251 (MQ=255)                                                                                      ggAATGCCGCGCCCGTCATCCTGTTCAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGAGTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTc                                                                                                                                                 <  1:137499/251‑1 (MQ=255)                                                                                                    gTCATCCTGTACAGAGACAGAGTGATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGt                                                                                                                                   <  2:48085/251‑1 (MQ=255)                                                                                                         ccTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATc                                                                                                                              <  1:376451/251‑1 (MQ=255)                                                                                                                      agagTTATCGGCGTGAATTGTGACGATAATTTCTTTACCGTGACCCGCGAGCGCTTCGTCGAGAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGCGCCGTCCTCCGTGTCGCCGATATAAATAACCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATAATGGAGGAGTCATAAGAATTCGACATCAACGTTTCTTGCTCATTTATACTTGGGTTAATCCGGTATTTTAccc                                                                                                                 >  2:226372/1‑251 (MQ=255)                                                                                                                              cGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGtt                                                                                                                          <  1:283288/234‑1 (MQ=255)                                                                                                                              cGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATAATTGGGTTAATCCGtt                                                                                                                          >  2:283288/1‑234 (MQ=255)                                                                                                                                                               gACCCGCGAGCGCTTCGTCGGTAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTc                                                                        >  1:249917/1‑251 (MQ=255)                                                                                                                                                                   cgcgAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCtat                                                                     >  1:41804/1‑250 (MQ=255)                                                                                                                                                                    gcgAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATaa                                                                   <  2:41804/251‑1 (MQ=255)                                                                                                                                                                      gAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACg                                                                 >  2:83138/1‑251 (MQ=255)                                                                                                                                                                                                 caccTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACaa                                        >  1:213170/1‑249 (MQ=255)                                                                                                                                                                                                 caccTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAg                                      <  2:131532/251‑1 (MQ=255)                                                                                                                                                                                                    cTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGAcc                                   >  2:66096/1‑251 (MQ=255)                                                                                                                                                                                                       gAATACCATGGGGTGCAGAACGGTGCAGTCATCCGTGTAGCAGATATATATACACGGGCAATTACGCACAGAATCCAGCCCTTTCAGGAATTTGATACTGGAGGAGTTATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGAATTATACGTTATTTTACCGTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCAATGAGCGCGATGACAAAAACCtgg                                >  2:215220/1‑250 (MQ=255)                                                                                                                                                                                                                                      tCCGTGTCGCCGATATATATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAagc  >  2:296532/1‑250 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         GTTATCGTCAAATTTACCGCCTGCGTGCAGAACGGTCATGATCACTTCCGCCGCCGATACGCCCTCTTCCGGGTGAATACCGGTCGGAATGCCGCGCCCGTCATCCTGTACAGAGACAGAGTTATCGGCGTGAATGGTGACGATAATTTCTTTACAGTGACCCGCGAGCGCTTCGTCGATAGCGTTATCTACCACCTCGAATACCATGTGGTGCAGACCGGTGCCGTCATCCGTGTCGCCGATATACATACCCGGGCGCTTACGCACCGCATCCAGCCCTTTCAGGACTTTGATACTGGAGGAGTCATAAGAATTCGACATCAACGTTTCTCGCTCATTTATACTTGGGTTAATCCGTTATTTTACCCTTTTCCACGGTAAACATCTTCGAATTTTCGTCCGACATGTCTATAACGTGTTCAGCACTGATCGCGCTGACAAAGACCTGTGATTGCGTCGCTTTTAAGCGGCTGGCAAGC  >  NC_000913/3879799‑3880277

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑