Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,303,924	T→A	9.1%	V187D (GTC→GAC)	cysS →	cysteine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,303,924	0	T	A	9.1%	120.4 / 3.5	44	V187D (GTC→GAC)	cysS	cysteine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (21/19);  new base A (2/2);  total (23/21)
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.48e-01

AATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGA  >  NC_002947/3303779‑3304071                                                                                                                                                  |                                                                                                                                                    aaTGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTc                                                                                                                                                    >  4:19107/1‑147 (MQ=255)      cGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCgg                                                                                                                                             <  4:350165/149‑1 (MQ=255)       gACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCGGCTAGAATTCGTCCTGTGGAAGGGCGTCAAGCCggg                                                                                                                                            >  5:356360/1‑149 (MQ=255)          gtgtACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGa                                                                                                                                         <  4:134337/149‑1 (MQ=255)              actacGGGGTTGGCTAGTTCGTCGGCTACGGCATGGTGTCGCGGAAGAGGATCGAAGTCCTGCGTTTGGGTGCACGGTTGGAGGTCGTCGAAGGCATGCAGGACGCGGTCGACTTGGTGGTGTGGTAGGGCGTCAAGCCGGGCGAGcc                                                                                                                                      <  6:355050/148‑1 (MQ=255)                tacCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAg                                                                                                                                   >  6:273013/1‑149 (MQ=255)                    gggTTGGGACGTTCGTGGGCGTCGGCATGCTGTCGCCCAATAGGATCGCACACCTGCGTTTCTGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGTGCGTCAAGCCGGGCGAGCCAAGCTgg                                                                                                                               <  6:304781/149‑1 (MQ=255)                    gggTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTgg                                                                                                                               <  2:97694/149‑1 (MQ=255)                                  gtcgGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATgggg                                                                                                                 <  3:359077/149‑1 (MQ=255)                                            gCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGc                                                                                                       <  1:195442/149‑1 (MQ=255)                                              aaGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGccc                                                                                                     <  8:180471/149‑1 (MQ=255)                                                    tCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTg                                                                                               <  5:273013/149‑1 (MQ=255)                                                        gcAAGAGGATCGAAGAACTGAGTATCGGTTAAAGCATCGAGGTCGACGAAGAAAAGCAGGACCCGCACGACTTCGTCCTGTCGAAGGGCGTCAAGCCGGGCGAGCAAAGATTGGAAACGCCATGGGGCCAGAGTAGCCCGGGATGGca                                                                                            >  4:242660/1‑148 (MQ=255)                                                              ggATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAg                                                                                     >  4:217441/1‑149 (MQ=255)                                                              ggATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAg                                                                                     <  3:19107/149‑1 (MQ=255)                                                                    aaGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCg                                                                               >  4:52388/1‑149 (MQ=255)                                                                     aGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGGTCgg                                                                              >  8:323292/1‑149 (MQ=255)                                                                           gagTATCGGAGCACCCATCGAGCTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTCGAAGGGCGTCAAGCCCGGCGAGCCAAGCTGCGAAACCCCATGGGCCCCGAGTCGCCCCGGCTGGCACAACGACTGCTCGGACATGt                                                                        >  4:78587/3‑149 (MQ=255)                                                                                   gTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACctgc                                                                >  8:151467/1‑149 (MQ=255)                                                                                   gTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACctgc                                                                >  7:264064/1‑149 (MQ=255)                                                                                       aCGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGACAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGcc                                                            >  3:348085/1‑149 (MQ=255)                                                                                       aCGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGACAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGcc                                                            >  3:342682/1‑149 (MQ=255)                                                                                          cATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTgg                                                         <  3:217441/149‑1 (MQ=255)                                                                                                ggTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCgaga                                                   <  3:52388/149‑1 (MQ=255)                                                                                                      cgaAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTc                                              >  4:298485/1‑148 (MQ=255)                                                                                                            cAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGGCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTc                                       <  3:78587/149‑1 (MQ=255)                                                                                                             aaGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGCGCTTCGACATTCa                                      >  7:297304/1‑149 (MQ=255)                                                                                                                cAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACgg                                   >  8:327896/1‑149 (MQ=255)                                                                                                                     cccGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTggcg                              >  3:267346/1‑149 (MQ=255)                                                                                                                         cTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGt                                  <  7:199508/141‑1 (MQ=255)                                                                                                                         cTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGt                                  <  7:204844/141‑1 (MQ=255)                                                                                                                         cTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGt                                  >  8:199508/1‑141 (MQ=255)                                                                                                                         cTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGt                                  >  8:204844/1‑141 (MQ=255)                                                                                                                               ttCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCt                    >  8:157399/1‑149 (MQ=255)                                                                                                                                     cTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGtt              >  1:227300/1‑149 (MQ=255)                                                                                                                                        tgGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTccc           <  7:151467/149‑1 (MQ=255)                                                                                                                                        tgGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTccc           <  4:267346/149‑1 (MQ=255)                                                                                                                                        tgGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTccc           >  3:40927/1‑149 (MQ=255)                                                                                                                                          gAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGc         >  7:207374/1‑149 (MQ=255)                                                                                                                                            aGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGcac       >  6:234414/1‑149 (MQ=255)                                                                                                                                             gggCGACAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGcacc      <  4:348085/149‑1 (MQ=255)                                                                                                                                             gggCGACAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGcacc      <  4:342682/149‑1 (MQ=255)                                                                                                                                                 gTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCacga  <  4:40927/149‑1 (MQ=255)                                                                                                                                                 gTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCacga  <  7:327896/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    AATGGCGACGTGTACTACCGGGTTGGCAAGTTCGTCGGCTACGGCAAGCTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGA  >  NC_002947/3303779‑3304071

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