Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,004,747	A→G	intergenic (+90/‑21)	ycbF → / → pyrD	putative fimbrial chaperone YcbF/dihydroorotate dehydrogenase, type 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,004,747	0	A	G	100.0%	145.9 / NA	41	intergenic (+90/‑21)	ycbF/pyrD	putative fimbrial chaperone YcbF/dihydroorotate dehydrogenase, type 2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/26);  total (15/26)

TGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGT  >  NC_000913/1004612‑1004880                                                                                                                                        |                                                                                                                                      tGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCgg                                                                                                                                   <  1:560033/140‑1 (MQ=255) ggATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCgg                                                                                                                                   <  1:246871/139‑1 (MQ=255)   aTTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGa                                                                                                                                 <  1:1124920/140‑1 (MQ=255)   aTTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGa                                                                                                                                 <  1:442657/140‑1 (MQ=255)        ggcggGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCa                                                                                                                            <  1:448893/140‑1 (MQ=255)         gcggGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAg                                                                                                                           <  1:313334/140‑1 (MQ=255)         gcggGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAg                                                                                                                           <  1:420687/140‑1 (MQ=255)            ggTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGgag                                                                                                                        <  2:464986/140‑1 (MQ=255)            ggTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGgag                                                                                                                        <  2:271965/140‑1 (MQ=255)               tGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGt                                                                                                                     <  1:275001/140‑1 (MQ=255)                 aGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTc                                                                                                                   <  1:37540/140‑1 (MQ=255)                      aaaGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGta                                                                                                              <  2:957178/140‑1 (MQ=255)                       aagaagGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGtac                                                                                                             >  1:1148868/1‑140 (MQ=255)                       aagaagGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGtac                                                                                                             >  1:588591/1‑140 (MQ=255)                         gaagGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGtacta                                                                                                           <  1:359501/140‑1 (MQ=255)                             gAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTc                                                                                                                   >  1:992338/1‑128 (MQ=255)                             gAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTc                                                                                                                   <  2:992338/128‑1 (MQ=255)                             gAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTAcccc                                                                                                       >  1:80127/1‑140 (MQ=255)                                       tGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTa                                                                                             <  2:80127/140‑1 (MQ=255)                                       tGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTa                                                                                             <  2:1148868/140‑1 (MQ=255)                                       tGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTa                                                                                             <  2:26764/140‑1 (MQ=255)                                        gCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTaa                                                                                            <  2:843501/140‑1 (MQ=255)                                          ccTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAg                                                                                          >  1:679362/1‑140 (MQ=255)                                                  aGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTc                                                                                  <  2:679362/140‑1 (MQ=255)                                                    gTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCa                                                                                >  2:960637/1‑140 (MQ=255)                                                       cGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCt                                                                             <  2:588591/140‑1 (MQ=255)                                                        ggAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTc                                                                            >  1:961361/1‑140 (MQ=255)                                                               gTTTGCGTTTTTTTTGCCGCAGGTCAATTCCC‑TTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCaga                                                                    >  2:731777/1‑140 (MQ=255)                                                                tttGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCaga                                                                    >  1:826661/1‑140 (MQ=255)                                                                  tGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAgc                                                                  >  2:183224/1‑140 (MQ=255)                                                                                       cAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTc                                             <  1:183224/140‑1 (MQ=255)                                                                                       cAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTc                                             <  1:960637/140‑1 (MQ=255)                                                                                               ttttGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTa                                     >  1:203720/1‑140 (MQ=255)                                                                                                    gTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGc                                  <  2:521580/138‑1 (MQ=255)                                                                                                    gTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGc                                  >  1:521580/1‑138 (MQ=255)                                                                                                    gTCCGAACTCGCACATAATACGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCg                                <  2:826661/140‑1 (MQ=255)                                                                                                                        cGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTg            >  1:642167/1‑140 (MQ=255)                                                                                                                        cGCCCCCGGTTTGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTg            >  1:127867/1‑140 (MQ=255)                                                                                                                              cGGTTTGTACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTGTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCAc      >  2:1069823/1‑140 (MQ=255)                                                                                                                                  ttGTACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGt  <  1:1069823/140‑1 (MQ=255)                                                                                                                                  ttGCACGCCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGt  <  1:426997/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      TGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATCTGCCCTGAAACAGGTTCGGAAAACGTTTGCGTTTTTTTTGCCGCAGGTCAATTCCCTTTTGGTCCGAACTCGCACATAATACGCCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCCTTTTCCAGCTCGATCCAGAGCGCGCTCATGAGTTTACTTTTCAGCAATTACGCCGTATTACAGGAACGCCGTTTGAAGCACTGGT  >  NC_000913/1004612‑1004880

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