Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,004,746	C→A	11.8%	intergenic (+89/‑22)	ycbF → / → pyrD	putative periplasmic pilin chaperone/dihydro‑orotate oxidase, FMN‑linked

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,004,746	0	C	A	11.8%	160.5 / 10.6	51	intergenic (+89/‑22)	ycbF/pyrD	putative periplasmic pilin chaperone/dihydro‑orotate oxidase, FMN‑linked
Reads supporting (aligned to +/- strand):  ref base C (26/19);  new base A (3/3);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

GGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGTGC  >  NC_000913/1004603‑1004882                                                                                                                                                |                                                                                                                                         ggcgATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCgg                                                                                                                                     <  2:210550/149‑1 (MQ=255)    gATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCCGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGAACACCGGGaa                                                                                                                                  >  3:182615/1‑149 (MQ=255)         aTGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAg                                                                                                                             >  1:437815/1‑149 (MQ=255)         aTGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgggggc                                                                                                                                 <  3:212555/145‑4 (MQ=255)         aTGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCgggggc                                                                                                                                 >  4:212555/1‑142 (MQ=255)          tGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAgg                                                                                                                            >  3:31861/1‑149 (MQ=255)                      gTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCggggg                                                                                                                                  >  3:412477/1‑129 (MQ=255)                      gTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTGTGCAAACCggggg                                                                                                                                  <  4:412477/131‑3 (MQ=255)                       ttGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGtac                                                                                                               <  3:196728/149‑1 (MQ=255)                       ttGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGtac                                                                                                               <  3:198101/149‑1 (MQ=255)                              cAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCt                                                                                                        >  2:110962/1‑149 (MQ=255)                                         cAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTaaa                                                                                             >  4:432544/1‑149 (MQ=255)                                              tCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCt                                                                                        <  4:182615/149‑1 (MQ=255)                                                     tGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGgaga                                                                                                                         <  2:184917/109‑1 (MQ=255)                                                     tGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGgaga                                                                                                                         >  1:184917/1‑109 (MQ=255)                                                      gAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTc                                                                                    >  2:261993/1‑145 (MQ=255)                                                      gAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTc                                                                                    <  1:261993/145‑1 (MQ=255)                                                            ggTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATc                                                                          <  1:110962/149‑1 (MQ=255)                                                              ttCGGAAAACGTTTGCG‑TTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTc                                                                              >  1:55092/1‑142 (MQ=255)                                                              ttCGGAAAACGTTTGCG‑TTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTc                                                                              <  2:55092/142‑1 (MQ=255)                                                                 ggAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGaa                                                                                                                                  <  1:419353/88‑1 (MQ=255)                                                                 ggAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGaa                                                                                                                                  >  2:419353/1‑88 (MQ=255)                                                                  gAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATGCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgc                                                                    >  1:1905/1‑149 (MQ=255)                                                                  gAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgc                                                                    <  1:10091/149‑1 (MQ=255)                                                                              gTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGc                                                                                >  1:313930/1‑125 (MQ=255)                                                                              gTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGc                                                                                <  2:313930/125‑1 (MQ=255)                                                                              gTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGt                                                        >  4:226795/1‑149 (MQ=255)                                                                              gTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGt                                                        <  4:31861/149‑1 (MQ=255)                                                                               ttttttttGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGtt                                                       >  3:124255/1‑149 (MQ=255)                                                                                       gCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGTATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTc                                               <  4:59009/149‑1 (MQ=255)                                                                                        ccGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCa                                              <  4:358407/149‑1 (MQ=255)                                                                                                 aaTTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACg                                     >  4:42493/1‑149 (MQ=255)                                                                                                  aTTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGc                                    >  3:94757/1‑149 (MQ=255)                                                                                                  aTTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGc                                    >  3:180433/1‑149 (MQ=255)                                                                                                   ttCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGcc                                   >  4:63888/1‑149 (MQ=255)                                                                                                   ttCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGcc                                   >  2:468233/1‑149 (MQ=255)                                                                                                   ttCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGcc                                   >  2:382501/1‑149 (MQ=255)                                                                                                      ccTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTa                                <  2:437815/149‑1 (MQ=255)                                                                                                       cTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTAt                               >  3:73163/1‑149 (MQ=255)                                                                                                                             aataCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTAc                                      <  2:347004/120‑1 (MQ=255)                                                                                                                             aataCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTAc                                      >  1:347004/1‑120 (MQ=255)                                                                                                                             aataCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCa         <  4:364664/149‑1 (MQ=255)                                                                                                                               taCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGaa            >  2:302042/1‑144 (MQ=255)                                                                                                                               taCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGaa            <  1:302042/144‑1 (MQ=255)                                                                                                                               taCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACt       >  2:121672/1‑149 (MQ=255)                                                                                                                               taCGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACt       >  1:250939/1‑149 (MQ=255)                                                                                                                                    ccccGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGTGc  >  1:214650/1‑149 (MQ=255)                                                                                                                                     cccGGTTTGCACACCGGGAACCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCAt                                                            >  2:283280/1‑90 (MQ=255)                                                                                                                                     cccGGTTTGCACACCGGGAACCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCAt                                                            <  1:283280/90‑1 (MQ=255)                                                                                                                                      ccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgcgcgc                                                                <  3:36946/85‑1 (MQ=255)                                                                                                                                      ccGGTGTGCAAACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgcgcgc                                                                >  4:36946/1‑85 (MQ=255)                                                                                                                                                |                                                                                                                                         GGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGTGC  >  NC_000913/1004603‑1004882

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