Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	374,784:1	+G	9.3%	coding (74/1233 nt)	PP_0312 →	electron transfer flavoprotein subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	374,784	1	.	G	9.3%	157.8 / 11.8	43	coding (74/1233 nt)	PP_0312	electron transfer flavoprotein subunit alpha
Reads supporting (aligned to +/- strand):  ref base . (22/17);  new base G (2/2);  total (24/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 = 2.52e-01

GCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGGCGCCCAGCAGGGCACCGGCGCGGGCGGGC  >  NC_002947/374639‑374930                                                                                                                                                   |                                                                                                                                                   gCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GAt                                                                                                                                                 >  1:21154/1‑149 (MQ=255)     aTCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATg                                                                                                                                                >  7:180866/1‑146 (MQ=255)         aggaggACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCA‑TCGCGGCATGCAGCGGGTGCAGACgg                                                                                                                                         <  4:228590/149‑2 (MQ=255)                 gCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCcgc                                                                                                                                                     <  1:251144/130‑1 (MQ=255)                   aTCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAg                                                                                                                               >  6:141254/1‑149 (MQ=255)                    tCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGc                                                                                                                              >  7:91333/1‑149 (MQ=255)                      cACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCt                                                                                                                             <  8:91333/148‑1 (MQ=255)                       aCTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTg                                                                                                                            >  3:88321/1‑148 (MQ=255)                       aCTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCA‑TCGCGGCATGCAGCGGGTGCAGACggttac                                                                                                                                     >  4:42194/1‑134 (MQ=37)                       aCTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCA‑TCGCGGCATGCAGCGGGTGCAGACggttac                                                                                                                                     <  3:42194/139‑6 (MQ=37)                         tGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGAt                                                                                                                         <  8:180866/149‑1 (MQ=255)                          gTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATg                                                                                                                        >  4:8843/1‑149 (MQ=255)                           tCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATgg                                                                                                                       <  2:21154/149‑1 (MQ=255)                              cAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGAGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGc                                                                                                                    >  4:278205/1‑149 (MQ=255)                               aaaGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATgggg                                                                                                                     <  3:100751/147‑1 (MQ=255)                               aaaGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATgggg                                                                                                                     >  4:100751/1‑147 (MQ=255)                                      cAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGc                                                                                                            <  6:69797/149‑1 (MQ=255)                                            gCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGAT‑CACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATcc                                                                                                     >  4:190459/1‑149 (MQ=255)                                              tAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGCCCAGCTGGATGGGGCCCAACGGCATCATCCg                                                                                                    >  6:93197/1‑149 (MQ=255)                                                aaCGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGc                                                                                                   <  4:88321/148‑1 (MQ=255)                                                          tGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCcaca                                                                                        >  6:359992/1‑149 (MQ=255)                                                            cGGAGGTGCACTGATGAGCGACATTATCCGCCGCGAT‑CACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGc                                                                                     <  3:190459/149‑1 (MQ=255)                                                              gAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCg                                                                                    <  5:141254/149‑1 (MQ=255)                                                                                aCATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCCCATGCGATTGCCGCAGGCTTTGTc                                                                  >  8:261968/1‑149 (MQ=255)                                                                                    tATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTcggcc                                                              >  4:124393/1‑149 (MQ=255)                                                                                      tCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCCTCATCCGCAAGACCCCACATGCGATTGCCGCAGGCTTTGTcggcccg                                                            >  4:38081/1‑149 (MQ=255)                                                                                        cgccgcGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTcggcccggc                                                          >  8:306955/1‑149 (MQ=255)                                                                                          ccgcGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCccggccg                                                        <  4:350297/149‑1 (MQ=255)                                                                                                  cACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAg                                                <  8:340206/149‑1 (MQ=255)                                                                                                   aCGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAgc                                               <  3:278205/149‑1 (MQ=255)                                                                                                   aCGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAgc                                               >  3:176553/1‑149 (MQ=255)                                                                                                   aCGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAgc                                               <  4:176553/149‑1 (MQ=255)                                                                                                             tGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCg                                     >  1:183931/1‑149 (MQ=255)                                                                                                               gATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCa                                   <  8:83535/149‑1 (MQ=255)                                                                                                                     ccGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAgcggcg                             >  4:18581/1‑149 (MQ=255)                                                                                                                            cGTCTGCA‑TCGCGGCATGCAGCGGGTGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGGCGCCcag                        >  7:193614/1‑147 (MQ=255)                                                                                                                                 gCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGgcgc                            <  4:32948/138‑1 (MQ=255)                                                                                                                                 gCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGgcgc                            >  3:32948/1‑138 (MQ=255)                                                                                                                                  cACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGc                                                                                     >  6:108847/1‑80 (MQ=255)                                                                                                                                  cACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGc                                                                                     <  5:108847/80‑1 (MQ=255)                                                                                                                                  cACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCAGCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGGCGCCCAGCAGGGCAc                <  3:38081/149‑1 (MQ=255)                                                                                                                                   aCCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGGCGCCCAGCAGGGCAcc               >  3:286049/1‑149 (MQ=255)                                                                                                                                     ccGCTGCATGCCGC‑GATGAAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCAGCAAGAACCCACATGCGATTGCCCCAGGCTTTGTCGGCCCGGCCGGCATCAAGCCCATCGAACGCAGCGGCCCCCAGCACGGCACCgg             >  7:140028/1‑149 (MQ=255)                                                                                                                                       gCTGCATGCCGC‑GATGCAGACGCAACAGACCAGATGGATGGGGCCCATCGGCATCATCCGCAAGAACCCTCATGCGTTTGCCGCAGGCTTTGTCGGCCCGGGCGGCATCAAGCGCATCGACCGCAGAGGCGCCCAGCAGGGCACCTgcg           <  8:140028/149‑1 (MQ=255)                                                                                                                                                cgc‑gATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGGCGCCCAGCAGGGCACCGGCgcgggcgggc  <  2:183931/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GCTGATCGAGGAGGACGCATCCACTGTCCCAAAGCCGCAAACGGCTAAACGTGAACCTGCGGAGGTGCACTGATGAGCGACATTATCCGCCGCGATCCACGCGCCGAGTGGATCGCCCGTAACCGTCTGCACCCGCTGCATGCCGC‑GATGCAGACGCAACAGACCAGCTGGATGGGGCCCAACGGCATCATCCGCAAGAACCCACATGCGATTGCCGCAGGCTTTGTCGGCCCGGCCGGCATCAAGCGCATCGACCGCAGCGGCGCCCAGCAGGGCACCGGCGCGGGCGGGC  >  NC_002947/374639‑374930

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