Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,381,715	G→C	8.4%	D156H (GAC→CAC)	dinB →	DNA polymerase IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,381,715	0	G	C	8.4%	133.7 / 2.1	45	D156H (GAC→CAC)	dinB	DNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base G (20/21);  new base C (2/2);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

CTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCGGCTGGGCATCGAAACCTGCCTG  >  NC_002947/1381574‑1381858                                                                                                                                              |                                                                                                                                                  cTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGgcg                                                                                                                                            >  1:163443/1‑149 (MQ=255)  tGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTCACTGgcga                                                                                                                                           >  8:224850/1‑148 (MQ=255)   ggATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGgcgc                                                                                                                                           >  6:233814/1‑148 (MQ=255)        tCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTCACTGgcgatcttgg                                                                                                                                     <  4:49833/149‑8 (MQ=255)        tCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTCACTGgcgatcttg                                                                                                                                      <  7:224850/148‑7 (MQ=255)          aGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCAAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCaa                                                                                                                                   >  7:290883/1‑149 (MQ=255)            cGAAAGCCATAGGTATTCGGGCAGCGCCACGCGAATTGCCGAGGATATTCGCCGGCGTGTAGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCTCCTAACAATTTTCTGGCAAAGATCGCCAGTCACTGGCGCAAACCCAATg                                                                                                                                 >  2:40266/1‑149 (MQ=255)              aCAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATggg                                                                                                                               >  7:25720/1‑149 (MQ=255)                    aGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGgcgc                                                                                                                                           <  4:83131/131‑1 (MQ=255)                    aGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGgcgc                                                                                                                                           >  3:83131/1‑131 (MQ=255)                              gggCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTAccc                                                                                                               >  2:249990/1‑149 (MQ=255)                               ggCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTAccc                                                                                                               >  8:240065/1‑148 (MQ=255)                                    cgcCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAAcc                                                                                                                                      >  4:27970/1‑120 (MQ=255)                                    cgcCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAAcc                                                                                                                                      <  3:27970/120‑1 (MQ=255)                                    cgcCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATg                                                                                                        <  2:106627/149‑1 (MQ=255)                                       cACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAg                                                                                                     <  7:18059/149‑1 (MQ=255)                                       cACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAg                                                                                                     <  5:233814/149‑1 (MQ=255)                                       cACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAg                                                                                                     <  1:257498/149‑1 (MQ=255)                                       cACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAg                                                                                                     <  2:171716/149‑1 (MQ=255)                                         cgcgCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTg                                                                                                   <  3:16894/149‑1 (MQ=255)                                            gcATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGa                                                                                                 >  5:139783/1‑148 (MQ=255)                                              aTTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAAc                                                                                              >  4:23588/1‑149 (MQ=255)                                                  ccGAGGATATTCGCCGGAGTGTTGCCCGAACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCATAACAAGTTTCTGGCCAAGATCGCCAGTGAATGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATCAAGTGGAAAccttc                                                                                          >  8:89065/1‑149 (MQ=255)                                                   cGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAAcgttcg                                                                                         <  6:139783/149‑1 (MQ=255)                                                          atTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGc                                                                                  >  4:153632/1‑149 (MQ=255)                                                                    gtgtTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCg                                                                                                                        >  4:246804/1‑102 (MQ=255)                                                                    gtgtTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCg                                                                                                                        <  3:246804/102‑1 (MQ=255)                                                                    gtgtTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTg                                                                       >  5:197570/1‑149 (MQ=255)                                                                       ttGCCCGCACCCTGCATATCACCATGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGc                                                                     >  3:57/1‑148 (MQ=255)                                                                                cccTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAACC‑TGGCCACCGGCAAGGCTGCaa                                                           <  7:240065/149‑2 (MQ=255)                                                                                 ccTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATg                                                          <  3:23588/149‑1 (MQ=255)                                                                                     cATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATgggg                                                       <  2:163443/148‑1 (MQ=255)                                                                                        atCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATGgggtggg                                                   >  5:290766/1‑149 (MQ=255)                                                                                         tgaCCGTGTCGGCTGGCGTGGCGCCTAAGATGTTTGTGGGCAAGATCGCCAGTGACTGGCGCAAACCCTATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATGgggtgggt                                                  <  7:89065/147‑1 (MQ=255)                                                                                             cGTGTCGGCTGGCGTGGCGCCTAACAAGGTTCTGGCCAAGATCGCCAGTGACTGTCGCAAACCCAATGGGCTGTTCGTGATTTCCC‑CAAATTAAGTGGAAACGTTAGTTGCATCCATTGCCGGTGGCCAGGTTGCATGGTGTGTGTAAg                                               <  1:40266/149‑1 (MQ=255)                                                                                               tgtCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTg                                            <  6:197570/149‑1 (MQ=255)                                                                                                     cTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCCTCCAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGc                                       <  1:249990/149‑1 (MQ=255)                                                                                                                    aCAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCgg                       <  6:290766/149‑1 (MQ=255)                                                                                                                       aGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCGGCTg                    <  8:25720/149‑1 (MQ=255)                                                                                                                               gCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCg                                     <  6:93270/124‑1 (MQ=255)                                                                                                                               gCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCg                                     >  5:93270/1‑124 (MQ=255)                                                                                                                                ccAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCGGCTGGGCATCGa            >  2:218699/1‑148 (MQ=255)                                                                                                                                    gATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACAC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCGGATGGGCATCGAAAcct       >  3:168404/1‑149 (MQ=255)                                                                                                                                         ccAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCGGCTGGGCATCGAAAcctgcctg  >  5:141984/1‑149 (MQ=255)                                                                                                                                          cAGTGACTGGCGCAAACCCAATGGGCGGATCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATGGGTTGGGTAAGGTGACGGCCGACAAACTGACGCGGCTGGGCATCGAAAcctgcct   <  3:208017/147‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                  CTGGATGTCAGCGACAGCCAGTGGTATTCGGGCAGCGCCACGCGCATTGCCGAGGATATTCGCCGGCGTGTTGCCCGCACCCTGCATATCACCGTGTCGGCTGGCGTGGCGCCTAACAAGTTTCTGGCCAAGATCGCCAGTGACTGGCGCAAACCCAATGGGCTGTTCGTGATTACCC‑CAAATGAAGTGGAAACGTTCGTTGCAGCC‑TTGCCGGTGGCCAGGTTGCATGGGGTGGGTAAGGTGACGGCCGACAAACTGACGCGGCTGGGCATCGAAACCTGCCTG  >  NC_002947/1381574‑1381858

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