Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,932,640	+C	intergenic (‑107/+43)	rluA ← / ← minE	bifunctional tRNA pseudouridine(32) synthase TrmQ/ribosomal large subunit pseudouridine synthase RluA/cell division topological specificity factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,932,638	1	.	C	94.7%	120.1 / 0.3	38	intergenic (‑105/+45)	rluA/minE	bifunctional tRNA pseudouridine(32) synthase TrmQ/ribosomal large subunit pseudouridine synthase RluA/cell division topological specificity factor
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (20/16);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

GGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑G‑CC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTT  >  NC_002947/1932495‑1932781                                                                                                                                                  |                                                                                                                                                  ggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑G‑cccgc                                                                                                                                              <  8:329802/149‑4 (MQ=255)        tCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAggg                                                                                                                                        <  5:171954/148‑1 (MQ=255)        tCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAggg                                                                                                                                        >  6:171954/1‑148 (MQ=255)          gAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTg                                                                                                                                     <  5:145932/149‑1 (MQ=255)            agagGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGc                                                                                                                                    >  5:271268/1‑148 (MQ=255)              agGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCGACAAAGGATATGCGTCCCCTGTAGGGGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑Ag                                                                                                                                >  2:321995/1‑149 (MQ=255)                gATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑Agg                                                                                                                               >  8:73288/1‑148 (MQ=255)                  tCTGGACATTCGACAGCGGCATTGCAGGTATCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTACTACAAATGATATGCGTCCCATGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAg                                                                                                                            >  5:373033/1‑149 (MQ=38)                           tCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAAtt                                                                                                                   >  2:243834/1‑149 (MQ=255)                              aCAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTc                                                                                                                >  4:224598/1‑149 (MQ=255)                                          cAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGa                                                                                                     >  1:104248/1‑148 (MQ=255)                                                  aCGAAAATAGCAGTAAACGGTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCACAGGCAGCCCCAAATTGTCGACTACGCGGATCAGc                                                                                                <  1:243834/146‑1 (MQ=38)                                                    gAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCg                                                                                          >  8:191844/1‑149 (MQ=255)                                                            cAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCAACAAAGGATATGCCTCCCCAGTAGGAGCGGGCTT‑G‑CCACCGAACGGCTGCA‑AGGCAGCCCCAAATTGTCGACAACACGGATCATCGATCAGGCAGGGt                                                                                   >  8:100886/1‑148 (MQ=37)                                                              gTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGAt                                                                                <  3:224598/149‑1 (MQ=255)                                                               aaaccGCTGGGGCCGCTTGCCGGCCCTGCCGGGGCACGCCCGCCCCTACACCGCATCTGCGGCCCCTGTCGGCGGGGGTTTCC‑CC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATg                                                                               <  4:23168/145‑1 (MQ=255)                                                                        gggCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTTTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCt                                                                        <  5:108387/147‑1 (MQ=255)                                                                         ggCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCa                                                                     <  7:100886/149‑1 (MQ=255)                                                                             gCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCGGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGaa                                                                  >  6:328287/1‑148 (MQ=255)                                                                                        cTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGc                                                      >  5:383162/1‑149 (MQ=255)                                                                                           cgcgGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTCTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGccc                                                    <  8:262648/148‑1 (MQ=255)                                                                                            gcgGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTg                                                  <  7:73288/149‑1 (MQ=255)                                                                                               ggCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGtt                                               <  7:191844/149‑1 (MQ=255)                                                                                                      ccGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGc                                        >  2:139211/1‑149 (MQ=255)                                                                                                       cGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCt                                       <  1:139211/149‑1 (MQ=255)                                                                                                             tACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGtt                                               >  5:323987/1‑135 (MQ=255)                                                                                                             tACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGtt                                               <  6:323987/135‑1 (MQ=255)                                                                                                             tACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGatat                                 >  5:289921/1‑149 (MQ=255)                                                                                                              aCAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATg                                >  6:381629/1‑149 (MQ=255)                                                                                                               cAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATg                                >  3:361623/1‑148 (MQ=255)                                                                                                                aaaGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGca                              >  5:85509/1‑149 (MQ=255)                                                                                                                  aGGATATGCGTCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGcaca                            >  1:161049/1‑149 (MQ=255)                                                                                                                            tCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGc                  >  5:227176/1‑149 (MQ=255)                                                                                                                            tCCCCTGTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGACAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGc                  <  6:271268/149‑1 (MQ=255)                                                                                                                                  gTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGt                                                >  2:135985/1‑113 (MQ=255)                                                                                                                                  gTAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGt                                                <  1:135985/113‑1 (MQ=255)                                                                                                                                   tAGGAGCGGGCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTCGCCGATGtt           <  6:289921/149‑1 (MQ=255)                                                                                                                                   tAGGAGCGGGCTT‑GCCC‑GCGAAGGGCCGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGt            <  7:213104/148‑1 (MQ=255)                                                                                                                                            gCTT‑GCCC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATAttt  <  2:273649/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGCGGGCTT‑G‑CC‑GCGAAGGGCTGCA‑AGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCACATATTT  >  NC_002947/1932495‑1932781

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: ‑