Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,362,324	G→C	10.4%	D322H (GAT→CAT)	PP_2957 →	acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,362,324	0	G	C	10.4%	128.8 / 3.8	46	D322H (GAT→CAT)	PP_2957	acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (24/17);  new base C (2/3);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

ATGTGCGGTCAGCGTTGCAAGACGGGTTCAGGCCGGTGGTGCTGATTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGTAAGGGCCGCTGTGCTGTTGCGCGGGCGATAGAACGGCTGTGCAGGTG  >  NC_002947/3362190‑3362471                                                                                                                                       |                                                                                                                                                    aTGTGCGGTCAGCGTTGCAAGACGGGTTCAGGCCGGTGGTGCTGATTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCACATCAGGATcatgca                                                                                                                                       <  8:313425/149‑5 (MQ=255) aTGTGCGGTCAGCGTTGCAAGACGGGTTCAGGCCGGTGGTGCTGATTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCACATCAGGATcatgc                                                                                                                                        <  4:359525/148‑4 (MQ=255)      cGGTCAGCGTTGCAAGACGGGTTCAGGCCGGTGGTGCTGATTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAgccgcc                                                                                                                                   >  8:151061/1‑148 (MQ=255)        gTCAGCGTTGCAAGACGGGTTCAGGCCGGTGGTGTTGATTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGAtgtg                                                                                                                                                >  7:157323/1‑133 (MQ=255)        gTCAGCGTTGCAAGACGGGTTCAGGCCGGTGGTGTTGATTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGAtgtg                                                                                                                                                <  8:157323/133‑1 (MQ=255)        gTCAGCGTTGCAAGACGGGTTCAGGCCGGTGGTGCTGATTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCaa                                                                                                                                 >  2:204679/1‑148 (MQ=255)             cGTTGCAAGACGGGTTCAGGCCGGTGGTGCTGATTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGc                                                                                                                           <  6:274896/149‑1 (MQ=255)                  cAAGACGGGTTCAGGCCGGTGGTGCTGGTTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCt                                                                                                                      >  1:15882/1‑149 (MQ=255)                           ttCAGGCCGGTGGTGCTGATTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGGTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCa                                                                                                             <  4:152567/149‑1 (MQ=255)                                    gtgcagCTGATTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCAATGCATCTTGGTCACGGGTTGCGATCACGAATTCGATTACCTGAATGAACCTCATGTGCATCAAAGACGCCAAAAACGCATCCTGGACTGCCAGCATCTGc                                                                                                     >  2:135908/6‑148 (MQ=255)                                        tgCTGATTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGt                                                                                                 >  5:211859/1‑148 (MQ=255)                                         gCTGATTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGt                                                                                                 >  8:306760/1‑147 (MQ=255)                                         gCTGATTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTaa                                                                                               >  6:64371/1‑149 (MQ=255)                                             aTTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCACATCAGGATcatgc                                                                                                                                        <  3:363304/104‑4 (MQ=255)                                             aTTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCACATCAGGATcatgc                                                                                                                                        >  4:363304/1‑101 (MQ=255)                                                 gcagcTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGc                                                                                        <  7:151061/148‑1 (MQ=255)                                                       aCCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAAttt                                                                                 >  5:76327/1‑149 (MQ=255)                                                           tttCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGTACTGCCAGCATCTGCCAGTAAGCCATGc                                                                                        >  7:267196/1‑138 (MQ=255)                                                           tttCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAgc                                                                             >  4:73068/1‑149 (MQ=255)                                                           tttCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGAGCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAgc                                                                             <  6:157741/149‑1 (MQ=255)                                                               aCTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCAt                                                                         >  6:190425/1‑149 (MQ=255)                                                                     cTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAg                                                                   >  7:189600/1‑149 (MQ=255)                                                                        aaaGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGt                                                                  <  5:64371/147‑1 (MQ=255)                                                                        aaaGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAACCATGCCGAATTTCAGCGCATGa                                                                       >  7:34537/1‑142 (MQ=255)                                                                        aaaGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAACCATGCCGAATTTCAGCGCATGa                                                                       <  8:34537/142‑1 (MQ=255)                                                                        aaaGACCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGAATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCATCGACTGCCAGCATCTGCCAGTAAGCCATGCTGAATTTCAGCGCATGAGCAGttt                                                                >  8:156102/1‑149 (MQ=255)                                                                            ccccGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCa                                                                                                   >  1:274239/1‑110 (MQ=255)                                                                            ccccGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCa                                                                                                   <  2:274239/110‑1 (MQ=255)                                                                                  acaggATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGgcaggcag                                                      >  4:95918/4‑149 (MQ=255)                                                                                   cTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGc                                                     <  5:190425/149‑1 (MQ=255)                                                                                    tGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCCGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCg                                                    >  7:336790/1‑149 (MQ=255)                                                                                      gATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCggg                                                  >  2:130310/1‑149 (MQ=255)                                                                                        tCTTGGTCACGGGTTGCGATGCCGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGt                                                 >  8:113741/1‑148 (MQ=255)                                                                                        tCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGt                                                 >  3:331691/1‑148 (MQ=255)                                                                                             gTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGTAAggg                                            <  4:156314/148‑1 (MQ=255)                                                                                                 cGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGTAAGGGCCgct                                       >  5:142300/1‑149 (MQ=255)                                                                                                    gTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGTAAGGGCCgctgtg                                    <  7:348205/149‑1 (MQ=255)                                                                                                        cGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGTAAGGGCCgctgtgctg                                 <  7:306760/148‑1 (MQ=255)                                                                                                         gATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCggcaggca                                                       <  5:131159/125‑1 (MQ=255)                                                                                                         gATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGTAAGGGCCGCTGTGCTGtt                               <  6:211859/149‑1 (MQ=255)                                                                                                              cGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGTAAGGGCCGCTGTGCTGTTGcgcg                          <  6:142300/149‑1 (MQ=255)                                                                                                                 aTTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGTAAGGGCCGCTGTGCTGTTGCGCGGGc                       <  6:76327/149‑1 (MQ=255)                                                                                                                     gtgtACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGTAAGGGCCGCTGTGCTGTTGCGCGGGCGAt                    <  4:331691/148‑1 (MQ=255)                                                                                                                                 gATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGACGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGTAAGGGCCGCTGTGCTGTTGCGCGGGCGATAGAACGGCtgtg        >  4:260951/1‑148 (MQ=255)                                                                                                                                      tGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGTAAGGGCCGCTGTGCTGTTGCGCGGGCGATAGAACGGCTGTGCAGGTg  >  2:58462/1‑149 (MQ=255)                                                                                                                                      tGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGTAAGGGCCGCTGTGCTGTTGCGCGGGCGATAGAACGGCTGTGCAGGTg  >  5:345954/1‑149 (MQ=255)                                                                                                                                       |                                                                                                                                                    ATGTGCGGTCAGCGTTGCAAGACGGGTTCAGGCCGGTGGTGCTGATTAGCAGCTACCGTTTCACTCGTCTGAAAGCCCCGCACTGGATCTTGGTCACGGGTTGCGATGACGAATTCGTGTACCTGCATGATCCTGATGTGGATCACAGCCGCCAAAAACGCATCCTGGACTGCCAGCATCTGCCAGTAAGCCATGCCGAATTTCAGCGCATGAGCAGTTTCGGCAGGCAGCGGGTAAGGGCCGCTGTGCTGTTGCGCGGGCGATAGAACGGCTGTGCAGGTG  >  NC_002947/3362190‑3362471

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