Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,115,668	G→A	R305R (CGC→CGT)	lpxL ←	lauroyl acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,115,668	0	G	A	96.3%	73.9 / NA	27	R305R (CGC→CGT)	lpxL	lauroyl acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (17/9);  total (17/10)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCATGATGCATTTTTCGACCACTTTGTTCATCCACGCGGCGGTAGTTTCGGCATCATCCAGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCTTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATC  >  NC_000913/1115430‑1115905                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              gCCATGCGATGATCAGGTTTCGTTTCCAGGTTATAGGGGTGTCAGTTTCACAGGGTGACATAGCAATCCGCGGTTGGTGCGCCAGGCGAGGTGAACATAAGAAGAAAAGATAAGCACTCTAATTAGGCTCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTGTTATCGCTAAATCTGGGGGATGTGCTGTGGGAATGTCTGGTGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGAAGGAACg                                                                                                                                                                                                                                   <  1:370648/251‑1 (MQ=255) gCCAGGCGACGATCATGTTTCGTTTCCCGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTCATAACGTGAAGGAACg                                                                                                                                                                                                                                   >  1:134814/1‑251 (MQ=255)                ggTTTCGTTGCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACGCTAATTATGCTCCCGACTGCCAGGGGGCGCAATCCAGATAGCTTTTAGCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAACAGCCGGCATGTGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGtt                                                                                                                                                                                                                    <  1:491998/251‑1 (MQ=255)                                                                 aaTCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAAc                                                                                                                                                                                                  >  1:17257/1‑220 (MQ=255)                                                                                      gcgcGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCATGATGCATTTTTCGACCACTTTGTTc                                                                                                                                               >  1:13934/1‑250 (MQ=255)                                                                                       cgcgGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCACTCCAGAGAGCTTTTATAGCTAAATCAGGGGGATTTGCTGTGGTACTGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTACTAACGTGAAGGAACGCCTTCCGGGCGTGtttt                                                                                                                                                                                                                  >  1:334538/1‑182 (MQ=255)                                                                                          ngTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTACAGGGGGAGCAATCAAGAGAGCTTTTATAGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGAAGGAACGCCTTCCGGGCGTGttt                                                                                                                                                                                                                   >  1:55241/2‑178 (MQ=255)                                                                                           gTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAg                                                                                                                                                                                                              <  1:215/182‑1 (MQ=255)                                                                                                cATAAGAAGCAAAGATAAGCACACTCATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGAAGGAACg                                                                                                                                                                                                                                   >  1:391495/1‑156 (MQ=255)                                                                                                  tAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGAAGGAACGCCTTCCGGGCGTGGTTTAAAGCGACGGTGTAACCACATAGACTGCTCTGGTGCCATCATGATGCAGTTTTCGAc                                                                                                                                                         <  1:423509/228‑1 (MQ=255)                                                                                                                 aGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTg                                                                                                                                                                                                                                           <  1:466831/131‑1 (MQ=255)                                                                                                                   cacaCTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCNGAGAGCGGTTNTCGNNAAATNANGNGGATTTGCTGTGGTAATGCCGGATGCCATNCTGAAGCATCCGGCATGGGAGATNNAANNNCCNNAAGGa                                                                                                                                                                                                                                      <  2:378404/134‑1 (MQ=255)                                                                                                                   cacaCTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGCGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGCTGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGAAGGa                                                                                                                                                                                                                                      >  1:378404/1‑134 (MQ=255)                                                                                                                     cacTAATTATGCGCCCGACTTCCAGGGGGCGCAATACAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATCACGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCatgat                                                                                                                                                                     >  1:220041/1‑197 (MQ=255)                                                                                                                                        ttCCAGGGGGCGCAATCCAGATAGCTTTTATCGCTAAATCAGGGGGATGTGCGGTGGTAATTCCGGATTCCATTCTGAAGCGGCCGGCATGGGGGATTTAATAACGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACg                                                                                                                                                                                                         <  1:289129/142‑1 (MQ=255)                                                                                                                                             gggggCGAAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCATGATGCATTTTTCGACCACTTTGTTCATCCACGCGGCGGTAGTTTCGGCATCATCCAGTGGCGGAGAACACTCTggcggc                                                                                         >  1:136991/1‑249 (MQ=255)                                                                                                                                                 gcgcAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCATGATGCATTTTTCGACCACTTTGTTCATCCACGCGGCGGTAGTTTCGGCATCATCCAGTGGCGGAGAACACTCTGGCGGCAGCAt                                                                                    >  1:120215/1‑250 (MQ=255)                                                                                                                                                                    tATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGACATTCTGACGCCTCCGGCATGGGCGATTTAATAACGAGAAGGAAGGCTTTCAGGGCGTGTTTTAAAGCGACGGTGTAACCCCATATACTGc                                                                                                                                                                                      >  1:222215/1‑133 (MQ=255)                                                                                                                                                                    tATCGCTAAATCAGGGGGATTTGCTGCGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGCGACATAAACACGCGTTGGAACGCCTTCCGGGGGTGTTTTAAAg                                                                                                                                                                                                              >  1:229855/1‑109 (MQ=255)                                                                                                                                                                        gCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGCAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGGAACCACATATACTGCTCTGGTGCCATCATGATGCTTTTTTCGACCACTTTGTTCATACACGCGGCGGTAGTTTc                                                                                                                           >  1:457866/1‑188 (MQ=255)                                                                                                                                                                        gCTAAATCAGGGGGATTTGCTGTGGGAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGGGAAGGAACGCCGTCCGGGCGTGTTTTAAAGCTACTGTGTAACCACATATACTGCTCTGGTGCCATCTTGATGCATGTTTCGACCTCTGTCGTCAGCCACGCGGCGGTAGTTTCGGCATCATCCAGTGGCGGAGAACACTCTgg                                                                                             <  1:341548/218‑1 (MQ=255)                                                                                                                                                                            aaTCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCATGATGCATTTTTCGACCACTTTGTTCATCCACGCGGCGGTAGTTTCGGCATCATCCAGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGc                                                        <  1:468268/251‑1 (MQ=255)                                                                                                                                                                                     naTTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGCGATTTAATAACGTGATGGAACGCCTTACGGGCGTGTTTTAAAGCGACGGTGTCACCACATATCCTGCTCTGGTGCCATCATGCTGCATTTTTCGACCACTTTGTTCATCCACGCGGCGGtcg                                                                                                                               >  1:48497/2‑169 (MQ=255)                                                                                                                                                                                          gNNGTGGTANNGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAAc                                                                                                                                                                                                  <  2:17257/99‑1 (MQ=255)                                                                                                                                                                                             ntgGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCATGATGCATTTTTAGACCACTTTGTTCATCCACGCGGCGGTAGTTTCGGCATCATCCAGTGGCGGAGAACACTCTggc                                                                                            >  1:148946/2‑198 (MQ=255)                                                                                                                                                                                                            gCCATTCTGAAGCATCCGGCATGGGAGATTTAATAACGTGCAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCAAATATACTGATCTGGTGCCATCATGATGCATTTTTCGCCCACTTTGTTCCTCCACGCGGCGGTAGTTTCGGCATCATCCAGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCTTACGGCGTGGACCGACGGGCAcc                                >  1:344341/1‑243 (MQ=255)                                                                                                                                                                                                                              gcngggnANNTNNAATNACGTGAAGGAACGCCTGCCGGNNGGGTTTTAANGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCATGATGCATTTTTCGANCACTTTGTTNNTCNACGCGNCNNNANNNTCNGCATCATCCAGTGGCGGAGNACNNTCTGGCNGCAGCATAATCAATTGATACCCTTTg                                                                 <  2:120215/188‑1 (MQ=255)                                                                                                                                                                                                                                  gggAGATTTAATAACGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCATGATGCATTTTTCGACCACTTTGTTCATCCACGCGGCGGTAGTTTCGGCATCATCCAGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCTTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATc  >  1:434284/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              GCCAGGCGACGATCAGGTTTCGTTTCCAGTTTATAGGGGTGTCATTTTCACAGGGTGACATAGCAATCCGCTGTTGGTGCGCCAGGCGCGGTGAACATAAGAAGAAAAGATAAGCACACTAATTATGCGCCCGACTTCCAGGGGGCGCAATCCAGAGAGCTTTTATCGCTAAATCAGGGGGATTTGCTGTGGTAATGCCGGATGCCATTCTGAAGCATCCGGCATGGGAGATTTAATAGCGTGAAGGAACGCCTTCCGGGCGTGTTTTAAAGCGACGGTGTAACCACATATACTGCTCTGGTGCCATCATGATGCATTTTTCGACCACTTTGTTCATCCACGCGGCGGTAGTTTCGGCATCATCCAGTGGCGGAGAACACTCTGGCGGCAGCATAATCAATTGATACCCTTTGCCATCTGGCTTACGGCGTGGAACGAAGGGCACCAGACATGCGCCGGACATCCGTGCCAGCATC  >  NC_000913/1115430‑1115905

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