Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,370,019	T→C	*661R (TGA→CGA)  M1T (ATG→ACG) †	arnA → arnD →	fused UDP‑L‑Ara4N formyltransferase/UDP‑GlcA C‑4'‑decarboxylase undecaprenyl phosphate‑alpha‑L‑ara4FN deformylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,370,019	0	T	C	95.6%	144.6 / ‑2.9	45	*661R (TGA→CGA)  M1T (ATG→ACG) †	arnA arnD	fused UDP‑L‑Ara4N formyltransferase/UDP‑GlcA C‑4'‑decarboxylase undecaprenyl phosphate‑alpha‑L‑ara4FN deformylase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (18/25);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCATGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATG  >  NC_000913/2369885‑2370152                                                                                                                                       |                                                                                                                                      tAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGAcc                                                                                                                                   >  1:1519337/1‑139 (MQ=255)  agagCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCa                                                                                                                                  <  2:1241104/139‑1 (MQ=255)    agCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCaaa                                                                                                                                >  2:835590/1‑139 (MQ=255)           tAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCt                                                                                                                         >  2:316521/1‑139 (MQ=255)                  aGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCAt                                                                                                                  >  1:1550879/1‑139 (MQ=255)                   gCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCAtt                                                                                                                 <  1:225594/139‑1 (MQ=255)                           aaTGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGa                                                                                                         <  1:316521/139‑1 (MQ=255)                                ccACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCTCCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCt                                                                                                    <  1:1132469/139‑1 (MQ=255)                                    cGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCg                                                                                                >  2:1489009/1‑139 (MQ=255)                                       tGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTgg                                                                                             <  2:756603/139‑1 (MQ=255)                                          cTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCAc                                                                                          <  2:270430/139‑1 (MQ=255)                                                                     cAGGAAACCATCGACGAAACGCTGGATTTCTTGCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGa                                                               >  1:470006/1‑139 (MQ=255)                                                                       ggAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGaaa                                                             >  2:13180/1‑139 (MQ=255)                                                                       ggAAACCATCGACGAAACGCTGGATCTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGaaa                                                             >  2:867913/1‑139 (MQ=255)                                                                        gAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAAt                                                            >  2:114615/1‑139 (MQ=255)                                                                          aaCCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCa                                                                                                                   >  2:889650/1‑82 (MQ=255)                                                                          aaCCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCa                                                                                                                   <  1:889650/82‑1 (MQ=255)                                                                           aCCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAg                                                                                                                               >  2:462277/1‑69 (MQ=255)                                                                           aCCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAg                                                                                                                               <  1:462277/69‑1 (MQ=255)                                                                           aCCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCtt                                                         <  1:867913/139‑1 (MQ=255)                                                                              aTCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAg                                                      >  2:118998/1‑139 (MQ=255)                                                                               tCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGt                                                     <  1:114615/139‑1 (MQ=255)                                                                               tCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGt                                                     <  1:118998/139‑1 (MQ=255)                                                                               tCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGt                                                     <  2:497812/139‑1 (MQ=255)                                                                               tCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGt                                                     >  1:560149/1‑139 (MQ=255)                                                                               tCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGt                                                     <  1:1501414/139‑1 (MQ=255)                                                                                 gacgaAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGTGTCTGCTGGAAATCTTGAGTaa                                                   <  1:289904/139‑1 (MQ=255)                                                                                 gacgaAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTaa                                                   <  1:188012/139‑1 (MQ=255)                                                                                 gacgaAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTaa                                                   <  1:880239/139‑1 (MQ=255)                                                                                                    tgcTGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGc                                <  2:470006/137‑1 (MQ=255)                                                                                                       tGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGt                                                     <  2:1338121/115‑1 (MQ=255)                                                                                                       tGCGCACCGTTGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGt                                                     >  1:1338121/1‑115 (MQ=255)                                                                                                                ttGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGt                                                                       >  1:52555/1‑88 (MQ=255)                                                                                                                ttGATCTTACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGt                                                                       <  2:52555/88‑1 (MQ=255)                                                                                                                      ttACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGc              <  1:1489009/139‑1 (MQ=255)                                                                                                                      ttACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGc              <  1:13180/139‑1 (MQ=255)                                                                                                                      ttACGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGc              <  2:689698/139‑1 (MQ=255)                                                                                                                        aCGGATAAACCATCATGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCata                                              >  1:866706/1‑105 (MQ=255)                                                                                                                        aCGGATAAACCATCATGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCata                                              <  2:866706/105‑1 (MQ=255)                                                                                                                        aCGGATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGccc            >  2:498020/1‑139 (MQ=255)                                                                                                                           gATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGt                                                     >  2:65946/1‑95 (MQ=255)                                                                                                                           gATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGt                                                     <  1:65946/95‑1 (MQ=255)                                                                                                                           gATAAACCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGa         <  1:835590/139‑1 (MQ=255)                                                                                                                               aaCCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGc                                                                                            <  2:514264/52‑1 (MQ=255)                                                                                                                               aaCCATCACGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGc                                                                                            >  1:514264/1‑52 (MQ=255)                                                                                                                                  catcacGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATg  <  2:845806/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                      TAGAGCATCGTAAACCGAGCATCCGCAATGCCCACCGCTGCCTGGACTGGGAGCCGAAAATTGATATGCAGGAAACCATCGACGAAACGCTGGATTTCTTCCTGCGCACCGTTGATCTTACGGATAAACCATCATGACCAAAGTAGGCTTACGCATTGATGTCGATACCTTTCGTGGCACCCGTGAAGGCGTGCCGCGTCTGCTGGAAATCTTGAGTAAGCATAATATTCAGGCCAGCATTTTTTTCAGCGTCGGCCCGGACAATATG  >  NC_000913/2369885‑2370152

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