Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,763,095	A→G	14.6%	D113G (GAT→GGT)	PP_3326 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,763,095	0	A	G	14.6%	111.3 / 11.7	48	D113G (GAT→GGT)	PP_3326	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/19);  new base G (3/4);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

ATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCGGGCCAGGCAAGACGCTG  >  NC_002947/3762947‑3763234                                                                                                                                                     |                                                                                                                                            aTCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGa                                                                                                                                             <  7:54043/149‑1 (MQ=255)  tCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGAt                                                                                                                                            >  5:84649/1‑149 (MQ=255)      gggcACGGTCGAGCGCCTGGCGGAGCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGGCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCATCCTTGACATGGGTGCGa                                                                                                                                        <  5:109939/146‑1 (MQ=255)        tCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCATCCTTGACATGGGTGCGATc                                                                                                                                      <  4:250323/149‑1 (MQ=255)                       tGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGc                                                                                                                       >  3:22237/1‑149 (MQ=255)                          cGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTg                                                                                                                    <  5:274425/149‑1 (MQ=255)                             aTCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGaaa                                                                                                                 <  4:24260/149‑1 (MQ=255)                                  gCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAg                                                                                                            >  1:26892/1‑149 (MQ=255)                                    cGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCATCCTTGACATGGGTGCGa                                                                                                                                        >  4:1547/1‑119 (MQ=255)                                    cGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCATCCTTGACATGGGTGCGa                                                                                                                                        <  3:1547/119‑1 (MQ=255)                                         tgcAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGc                                                                                                      >  5:250708/1‑148 (MQ=255)                                         tgcAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGc                                                                                                      >  5:273834/1‑148 (MQ=255)                                                aCTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTAcc                                                                                              <  4:22237/149‑1 (MQ=255)                                                       ccGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCgg                                                                                       >  2:94008/1‑149 (MQ=255)                                                                gCACTTGCCGTGCTGCTGGACAATGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCt                                                                              >  2:181784/1‑149 (MQ=255)                                                                  aCTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGcc                                                                                                          <  2:211737/119‑1 (MQ=255)                                                                  aCTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGcc                                                                                                          >  1:211737/1‑119 (MQ=255)                                                                      gCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGa                                                                         >  1:198982/1‑148 (MQ=255)                                                                      gCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGa                                                                         >  3:98166/1‑148 (MQ=255)                                                                        cGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAAc                                                                      <  8:84220/149‑1 (MQ=255)                                                                        cGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAAc                                                                      <  2:26892/149‑1 (MQ=255)                                                                          tgctgcTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTg                                                                    >  2:177977/1‑149 (MQ=255)                                                                              gctgGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAAc                                                                      <  8:212991/143‑1 (MQ=255)                                                                              gctgGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAAc                                                                      >  7:212991/1‑143 (MQ=255)                                                                               ctgGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAAAGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGAc                                                                 <  5:272017/147‑1 (MQ=255)                                                                                  gACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACAtgct                                                            >  3:13880/1‑149 (MQ=255)                                                                                     aaCGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCg                                                         <  1:94008/149‑1 (MQ=255)                                                                                        gggCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGttt                                                      >  3:88907/1‑149 (MQ=255)                                                                                            cTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGc                                                                                                           >  4:1003/1‑92 (MQ=255)                                                                                            cTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGc                                                                                                           <  3:1003/92‑1 (MQ=255)                                                                                               ggTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGAt                                                <  1:177977/148‑1 (MQ=255)                                                                                                 tCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACGTGCTGCGTTTGGCGa                                                 >  5:293004/1‑145 (MQ=255)                                                                                                 tCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACGTGCTGCGTTTGGCGa                                                 <  6:293004/145‑1 (MQ=255)                                                                                                     aGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACAtgctgc                                                          >  8:297999/1‑132 (MQ=255)                                                                                                     aGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACAtgctgc                                                          <  7:297999/132‑1 (MQ=255)                                                                                                     aGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTccc                                         <  6:250708/149‑1 (MQ=255)                                                                                                              aCAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGAc                                <  6:84649/149‑1 (MQ=255)                                                                                                                aaGCAATGGGTCAATCTGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGAccc                              >  6:158405/1‑149 (MQ=255)                                                                                                                            aaTCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCgg                  >  6:156995/1‑149 (MQ=255)                                                                                                                              tCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACAtgct                                                            >  6:7463/1‑105 (MQ=255)                                                                                                                              tCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCg                                                         >  5:149259/1‑108 (MQ=255)                                                                                                                              tCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCg                                                         <  6:149259/108‑1 (MQ=255)                                                                                                                                  gATCGCATCCTTGACATGGGTGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGa                                                 >  7:230006/1‑112 (MQ=255)                                                                                                                                  gATCGCATCCTTGACATGGGTGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGa                                                 <  8:230006/112‑1 (MQ=255)                                                                                                                                  gATCGCATCCTTGACATGGGTGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCGGGCCAgg            >  2:97240/1‑149 (MQ=255)                                                                                                                                  gATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCGGGCCAgg            <  4:13880/149‑1 (MQ=255)                                                                                                                                    tCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCGGGCCAGGc           >  2:13304/1‑148 (MQ=255)                                                                                                                                        aCCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCGGGCCAGGCAAGAc      <  5:158405/149‑1 (MQ=255)                                                                                                                                            aTGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCGGGCCAGGCAAGACGCTg  <  5:74067/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                            ATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAGGATGCGATCGCCTCCACCCACAGCTTGAAATACAGCCTGGCGATCTACCAGAGCGGTGCACGCCTGCCGAAACTGGACATGCTGCGTTTGGCGATCGTTCCCCTGAGCGACCCTCTGCAAGTCGGGCCAGGCAAGACGCTG  >  NC_002947/3762947‑3763234

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