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	93.3%	90.8 / 0.3	30	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 (14/14);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTG‑CC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCAC  >  NC_002947/1932491‑1932775                                                                                                                                                      |                                                                                                                                           tGGCGGCGTCTTCGAAGAGGATCTGCACATTCGACAGCTGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTG‑c                                                                                                                                           <  8:184971/149‑1 (MQ=255)   gcggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTG‑ccc                                                                                                                                         >  4:213659/1‑148 (MQ=255)   gcggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCCCTACTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTG‑ccc                                                                                                                                         >  3:291481/1‑148 (MQ=255)   gcggcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCGTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTG‑ccc                                                                                                                                         <  1:149449/149‑2 (MQ=255)      gcgTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTG‑cccgcg                                                                                                                                      >  8:265914/1‑145 (MQ=255)        gTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTG‑cccgcgta                                                                                                                                    <  7:265914/149‑7 (MQ=255)        gTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCGCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGaa                                                                                                                                    <  8:98487/149‑1 (MQ=255)             cGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCt                                                                                                                               >  8:1473/1‑149 (MQ=255)                 gagGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGc                                                                                                                             >  3:74576/1‑147 (MQ=255)                     aTCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAgg                                                                                                                         >  3:59217/1‑147 (MQ=255)                     aTCTGCACATTCCACAGCGGCATTGCAGGGCTCACCAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGCATATGCGCCCCCTGTCGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCa                                                                                                                       <  3:213659/149‑1 (MQ=37)                                 gACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCCCTACTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGACGGGCTGCCATGCAGCCCCAAATTGt                                                                                                           >  1:130424/1‑149 (MQ=38)                                    cgCGGCATGCCAGGCCTCACGAAAATCGCAGTAACCGCTGGGGCCGCTTTGCGGCCCTGCGCGGGCAACCCCGCTCCTACAAAGCATATGCGTCCCCTGTCGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGc                                                                                                                     <  3:111770/135‑1 (MQ=37)                                    aGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCGCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCGTGCCC‑GCGAAGGGCTGCAAGGCCGc                                                                                                                     >  4:111770/1‑136 (MQ=38)                                      cGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGAc                                                                                                       >  8:57396/1‑148 (MQ=255)                                           ttGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTAAAAAGGATAAGGGTGCCCAGTAGGACC‑GGGCGTGCCC‑GCGAACGGCAGCAAGGCAGCCCCAAATTTTAGACAAagcg                                                                                                 >  2:3251/1‑149 (MQ=25)                                              cAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGAt                                                                                              >  5:317628/1‑149 (MQ=255)                                                    tcACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCg                                                                                         >  1:31852/1‑148 (MQ=255)                                                      aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATc                                                                                      <  6:317628/149‑1 (MQ=255)                                                      aCGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCCGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGGTCAGCGATc                                                                                      >  6:262331/1‑149 (MQ=255)                                                          aaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTg                                                                                                                                             >  8:148256/1‑91 (MQ=255)                                                          aaaTAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTg                                                                                                                                             <  7:148256/91‑1 (MQ=255)                                                           aaTAGCAGTAACCGGGGGGGCCGTTTTGGGGCCGTTCGCGGGAAAGCCCGCTCATACAAAGTATATGCGGCCCCTTTCGGCGG‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCa                                                                                 <  4:59217/149‑1 (MQ=255)                                                                             ggcggcTTTGGGGCCCTTCGCGGTCAAGCCCGCTCATAAAAAGTATATGCGTCCCCTGTAGTAGG‑GGGTTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCa                                                               <  2:289957/149‑1 (MQ=255)                                                                                ggttttGGGGCCTTTCGCGGGCAAGCCGGTTCCTACAAAGTATATGCTTCCCTTGTTGGAGG‑GGGTTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGaa                                                            <  4:291481/146‑1 (MQ=255)                                                                                         gCCCTTCGCGGGCAAGCCCCCGCCTAAAAAGGATATGCGTACCCGGTAGGAGCGGGGGTTG‑CC‑CCGAACGGCTGCAAGGCAGCCCCAAAATGTCGACAACGCCGCTAAGCCATACGGCAGGGTGATTTTGAGCGACAGAATCCaacagc                                                   >  2:192977/1‑149 (MQ=255)                                                                                                    gCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGttt                                        >  6:244832/1‑149 (MQ=255)                                                                                                    gCAAGCCCGCTCCTACAAAGGATATGCGTCACCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCAAAATTTTCGACTACGCGGATCAGCGATACGTCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGccaagggta                                        >  8:268717/1‑142 (MQ=38)                                                                                                      acccccGCTCCTACAAAGCTTTTGCGTCCCCTGTAGGCGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCTGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                      <  7:51867/146‑1 (MQ=38)                                                                                                      aaGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTc                                      <  7:57396/149‑1 (MQ=255)                                                                                                           cGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCt                                 <  7:1473/149‑1 (MQ=255)                                                                                                            gCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAACTc                                >  8:145019/1‑149 (MQ=255)                                                                                                                cTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGa                              <  2:31852/147‑1 (MQ=255)                                                                                                                    aaaGGATATGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGAt                                                                          <  2:106062/99‑1 (MQ=255)                                                                                                                            tGCGTCCCCTGTAGGAGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACatcatc                 <  4:74576/148‑1 (MQ=255)                                                                                                                                       tagtagC‑GGGCTTG‑CCAGCGAAGGTCTGCAAGGAAGACCCAAATTGTCGACTACGCGGATCAGAGATCCGTCAGGTTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATAATCGTTGCCGATGt      <  7:268717/148‑1 (MQ=255)                                                                                                                                          gAGC‑GGGCTTGCCC‑GCGAAGGGCCGCAAAGCGGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGTTATGCACATCATCGTTGCCGATGTTcac  <  1:19957/149‑1 (MQ=39)                                                                                                                                           aGC‑GGGCTTGCCC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGtt     >  5:275136/1‑145 (MQ=255)                                                                                                                                                      |                                                                                                                                           TGGCGGCGTCTTCGAAGAGGATCTGGACATTCGACAGCGGCATTGCAGGTCTCACGAAAATAGCAGTAAACGCTGGGGCCGCTTTGCGGCCCTTCGCGGGCAAGCCCGCTCCTACAAAGGATATGCGTCCCCTGTAGGAGC‑GGGCTTG‑CC‑GCGAAGGGCTGCAAGGCAGCCCCAAATTGTCGACTACGCGGATCAGCGATCCGGCAGGGTGATGTTCAGCTCCAGAATCGAGCAGCTGCCCTGGTTTTCCAGCTCGATATGCACATCATCGTTGCCGATGTTCAC  >  NC_002947/1932491‑1932775

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