Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,344,839	T→C	11.3%	D97G (GAC→GGC)	PP_2941 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,344,839	0	T	C	11.3%	111.4 / 5.3	46	D97G (GAC→GGC)	PP_2941	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/16);  new base C (3/2);  total (27/19)
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

TACCAGCCGTCGATCGGCGTTATGCCTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGCCAGGCCTCTGGCAGGCTCAGGGTTTCGACAAAGTCGTGTATGCCT  >  NC_002947/3344698‑3344982                                                                                                                                              |                                                                                                                                                 tACCAGCCGTCGATCGGCGTTATGCCTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGCccgcggg                                                                                                                                           >  4:167677/1‑146 (MQ=255) tACCAGCCGTCGATCGGCGTTATGCCTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGCccgcggg                                                                                                                                           >  4:164249/1‑146 (MQ=255)    cAGCCGTCGATCGGCGTTATGCCTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGTACGCCCAATGCGTTCGCCAGCCCCAAAACACCAGGTCCGCCCCCgg                                                                                                                                        >  2:249034/1‑149 (MQ=255)    cAGCCGTCGATCGGCGTTATGCCTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCgg                                                                                                                                        <  3:262200/149‑1 (MQ=255)    cAGCCGTCGATCGGCGTTATGCCTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCgg                                                                                                                                        >  2:109394/1‑149 (MQ=255)    cAGCCGTCGATCGGCGTTATGCCTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCgg                                                                                                                                        >  2:113337/1‑149 (MQ=255)         gTCGATCGGCGTTATGCCTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGCCCGCGGGCGGACC‑Tg                                                                                                                                  <  3:167677/149‑1 (MQ=255)         gTCGATCGGCGTTATGCCTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGCCCGCGGGCGGACC‑Tg                                                                                                                                  <  3:164249/149‑1 (MQ=255)            gATCGGCGTTATGCCTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCtt                                                                                                                               <  1:76129/149‑1 (MQ=255)               ggctggtATGCCTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCAt                                                                                                                            <  3:173930/143‑1 (MQ=255)                 gCGTTATGCCTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTCTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGACGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATcc                                                                                                                          >  2:96252/1‑149 (MQ=255)                         ccTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCGTGGCGTTGAACGGCGCGGCACGGTCCGTCGCCAAATGCGTACCCCACCCCCAGAACACCAGGGCCGCCCGCGGGCC‑TGGTACATCCAccgcc                                                                                                                    >  6:276713/1‑144 (MQ=255)                         ccTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGc                                                                                                                    >  6:251523/1‑147 (MQ=255)                             ggTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGt                                                                                                              <  1:113337/149‑1 (MQ=255)                             ggTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGt                                                                                                              <  1:109394/149‑1 (MQ=255)                              gTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTg                                                                                                             >  7:187363/1‑149 (MQ=255)                                     cAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAaga                                                                                                                >  6:154818/1‑139 (MQ=255)                                     cAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAaga                                                                                                                <  5:154818/139‑1 (MQ=255)                                      aGATCGCCCGGAAGAAGTGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCa                                                                                                          <  7:212165/144‑1 (MQ=255)                                      aGATCGCCCGGAAGAAGGGGCCGTGTGCCAACTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCTCCAGCCACATCTCACCAGGTACGCCCGCGGGCC‑TGCTTCATCCACCAGCAAGAATTCCAcaga                                                                                                      >  8:212165/1‑145 (MQ=255)                                             ccGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCg                                                                                              >  2:89903/1‑149 (MQ=255)                                                      ggggCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCAc                                                                                                                        >  4:119200/1‑114 (MQ=255)                                                       gggCCGTGTGCCACCATTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGtt                                                                                    <  1:249034/149‑1 (MQ=255)                                                                 ccACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGcc                                                                           >  1:59020/1‑148 (MQ=255)                                                                   aCCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAg                                                                        >  3:225387/1‑149 (MQ=255)                                                                    ccTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGt                                                                       >  7:213702/1‑149 (MQ=255)                                                                    ccTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGt                                                                       >  6:238455/1‑149 (MQ=255)                                                                        ttCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGTCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGc                                                                    <  1:96252/148‑1 (MQ=255)                                                                               cTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGcc                                                            >  8:12053/1‑149 (MQ=255)                                                                                    gTTGATCGGCGCGGCAGGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACgt                                                       <  5:276713/149‑1 (MQ=255)                                                                                    gTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACgt                                                       <  5:251523/149‑1 (MQ=255)                                                                                            gcgcGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGc                                               >  2:237622/1‑149 (MQ=255)                                                                                                                   tGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGCCAGGCCTCTGGCAGGCTCAGGGt                        >  6:95810/1‑149 (MQ=255)                                                                                                                     cGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGCCAGGCCTCTGGCAGGCTCAGGGttt                      >  1:21095/1‑149 (MQ=255)                                                                                                                       tccgcatgCCCCTCCTCACATGGGCCGCCCGCGGGCCGTCCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGCCAGGCCTCTGGCAGGCTCAGGGTTTc                     <  1:89903/142‑1 (MQ=255)                                                                                                                             aGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACgt                                                       <  4:213775/108‑1 (MQ=255)                                                                                                                             aGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACgt                                                       >  3:213775/1‑108 (MQ=255)                                                                                                                                 ccATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGCCAGGCCTCTGGCAGGCTCAGGGTTTCGACAAAGTCgt          <  4:225387/149‑1 (MQ=255)                                                                                                                                 ccATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGCCAGGCCTCTGGCAGGCTCAGGGTTTCGACAAAGTCgt          >  4:158192/1‑149 (MQ=255)                                                                                                                                    tctcACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGCCAGGCCTCTGGCAGGCTCAGGGtt                       >  2:226455/1‑133 (MQ=255)                                                                                                                                    tctcACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGCCAGGCCTCTGGCAGGCTCAGGGtt                       <  1:226455/133‑1 (MQ=255)                                                                                                                                      tcACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGCCAGGCCTCTGGCAGGCTCAGGGTTTCGACAAAGTCGTGTATg     >  2:180946/1‑149 (MQ=255)                                                                                                                                        agcAGGCCCGCGGGCGGACC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGCCAGGCCTCTGGCAGGCTCAGGGTTTCGACAAAGTCGTGTATGcc   >  6:68930/3‑149 (MQ=255)                                                                                                                                         ccAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGCCAGGCCTCTGGCAGGCTCAGGGTTTCGACAAAGTCGTGTATGCCt  <  1:180946/149‑1 (MQ=255)                                                                                                                                             gTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGCCAgg                                           <  5:241351/104‑1 (MQ=255)                                                                                                                                             gTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGCCAgg                                           >  6:241351/1‑104 (MQ=255)                                                                                                                                              |                                                                                                                                                 TACCAGCCGTCGATCGGCGTTATGCCTCGGTGCGGCCAGATCGCCCGGAAGAAGGGGCCGTGTGCCACCTTTTCGACCCTGGCGTTGATCGGCGCGGCACGGTCGGTCGCCCAATGCGGTCGCCAGCCCCATCTCACCAGGTCCGCCCGCGGGCC‑TGCTTCATCCACCCGCAAGAGTGCCACCGGCGTGGTCGGCGCGACGTTGTACCGGCCAAGTGGCTGGTCGCCGACGTTGTTCTGCCAGGCCTCTGGCAGGCTCAGGGTTTCGACAAAGTCGTGTATGCCT  >  NC_002947/3344698‑3344982

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