Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,239,100	G→T	8.7%	intergenic (‑75/‑300)	argF ← / → PP_1080	ornithine carbamoyltransferase/molybdopterin‑binding oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,239,100	0	G	T	8.7%	120.6 / 2.1	45	intergenic (‑75/‑300)	argF/PP_1080	ornithine carbamoyltransferase/molybdopterin‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (21/20);  new base T (2/2);  total (23/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 = 3.89e-01

GTGGCGGATCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTAT  >  NC_002947/1238957‑1239241                                                                                                                                                 |                                                                                                                                              gTGGCGGATCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGc                                                                                                                                          >  5:186622/1‑149 (MQ=255)       gATCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCTTACGCCGCagg                                                                                                                                    <  8:278712/149‑3 (MQ=255)         tCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCTTACGCCGCaggtc                                                                                                                                  <  1:238200/149‑5 (MQ=255)           aCCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTAATGAGGCACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGc                                                                                                                                <  8:213636/149‑1 (MQ=255)            ccccGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTTCGGCGTAAGACGCACGGGGCg                                                                                                                               >  3:222181/1‑149 (MQ=255)               cGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCgg                                                                                                                              <  5:302222/147‑1 (MQ=255)               cGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCgg                                                                                                                              >  6:302222/1‑147 (MQ=255)                gAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGAc                                                                                                                            >  3:189557/1‑148 (MQ=255)                  gCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGa                                                                                                                          >  5:17040/1‑148 (MQ=255)                      ttCGTCGGTGGTTAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAg                                                                                                                     <  1:119880/149‑1 (MQ=255)                       tcgtcgGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGt                                                                                                                                              <  2:278767/123‑1 (MQ=255)                       tcgtcgGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGt                                                                                                                                              >  1:278767/1‑123 (MQ=255)                          tcgGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGa                                                                                                                  >  2:183007/1‑148 (MQ=255)                                gtgAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGgcg                                                                                                            <  3:81141/148‑1 (MQ=255)                                   aattCCAGCAGGGTGATAAAGTGCCTAGCGCTCATGAGTGAATACCTTATCTGCAACGGTTTGCAGGTCGACGGTA‑TGGTTTTTTTTTACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCaa                                                                                                         <  6:17040/145‑1 (MQ=255)                                        cagcagGGCGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTAATCTGCAACGGTATGCAGGTAGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACta                                                                                                    <  8:62424/148‑1 (MQ=255)                                        cagcagGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCTCAATAGta                                                                                                    >  6:212365/1‑148 (MQ=255)                                        cagcagGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACtat                                                                                                   <  7:296058/149‑1 (MQ=255)                                          gcagGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATaa                                                                                                 <  1:183007/149‑1 (MQ=255)                                                                  ctatctgTGACTACCTTATCTGCAAGGGTTTGCAGGTCGTCCGTT‑TGGTTTTTTTTGACTACGGGGGTTAACTGCGGCGTAAGCGGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAggg                                                                         <  4:23777/143‑1 (MQ=255)                                                                         tGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTCCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGgagatcg                                                                  >  3:23777/1‑146 (MQ=255)                                                                             tACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAgcgg                                                              >  2:226555/1‑149 (MQ=255)                                                                               ccTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAgcggc                                                             >  6:289867/1‑148 (MQ=255)                                                                               ccTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAgcggc                                                             <  5:289867/148‑1 (MQ=255)                                                                               ccTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAgcgg                                                              <  1:226555/147‑1 (MQ=255)                                                                               ccTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACATGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAATTAAATGTCGCCTGAAGCCAAGGGGGGAAAAgcggcg                                                            >  1:278772/1‑149 (MQ=255)                                                                                cttttgtGCAACGGTAGGCCGGTCCACCGTA‑TGTTTTTTTTTTACAACGGGGGTCACGTGCGGGGTACGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGt                                                           <  2:278772/143‑1 (MQ=255)                                                                                 ttATCTGCAACGGTATGCAGGTCGACCGTAGTGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGt                                                           >  1:192389/1‑149 (MQ=255)                                                                                    tCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAtt                                                       >  2:282036/1‑149 (MQ=255)                                                                                        catgggTATTCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttg                                                      <  4:189557/142‑1 (MQ=255)                                                                                        cAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGtctatgca                                                    >  8:296543/1‑141 (MQ=255)                                                                                         aaCGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAAAACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgca                                                    >  7:115747/1‑147 (MQ=255)                                                                                                      tCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgttgcaaag                                     <  6:186622/149‑1 (MQ=255)                                                                                                       ccaCCGTC‑TGGTTTTTTTTCACAACGGGGGGTACCTTCGGCGTAAGCGGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgttgcaaagg                                    <  8:115747/147‑1 (MQ=255)                                                                                                        gcccGTA‑TGTTTTTGTTTGACATCGGGGGTGCCCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGACAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGGGAttgcaaagggtgttgcaaaggg                                   <  4:222181/147‑1 (MQ=255)                                                                                                        gACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAATTAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgttgcaaaggg                                   >  5:259836/1‑149 (MQ=255)                                                                                                                           gACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGcc                 >  7:5847/1‑149 (MQ=255)                                                                                                                              aaCGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGt              >  3:295230/1‑149 (MQ=255)                                                                                                                                  ggggTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTg          <  2:192389/149‑1 (MQ=255)                                                                                                                                     gTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgtt                                           <  4:285560/113‑1 (MQ=255)                                                                                                                                     gTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGAttgcaaagggtgtt                                           >  3:285560/1‑113 (MQ=255)                                                                                                                                         cgtGCGGCTTACGCCGCACGGGGTGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAATGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTa   >  2:187764/3‑149 (MQ=255)                                                                                                                                         cgtGCGGCTTACGCCGCACGGGGTGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAATGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTa   >  2:192412/3‑149 (MQ=255)                                                                                                                                         ccggcggcgTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTATATGTCGCCTGAAGCCAAGGGGTGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTa   <  5:212365/146‑1 (MQ=255)                                                                                                                                          cTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTAt  <  1:282036/149‑1 (MQ=255)                                                                                                                                          cTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTAt  >  6:201179/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                              GTGGCGGATCACCCCGAGCAATTCGTCGGTGGTGAAGTCCAGCAGGGAGAGAAAGTGCCTAGCGCTCATGAGTGACTACCTTATCTGCAACGGTATGCAGGTCGACCGTA‑TGGTTTTTTTTGACAACGGGGGTGACCTGCGGCGTAAGCCGCACGGGGCGGACGAAATAGGGAAAGGCGCAATACTATAAGTAAATGTCGCCTGAAGCCAAGGGGGGAAAAGCGGCGTGATTGCAAAGGGTGTTGCAAAGGGTGTCTTTACCCGTTGGCCCGTGCTGTTCTTTAT  >  NC_002947/1238957‑1239241

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