Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,033,298	G→T	10.9%	P312P (CCG→CCT)	PP_3555 →	transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,033,298	0	G	T	10.9%	116.7 / 4.4	43	P312P (CCG→CCT)	PP_3555	transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base G (13/25);  new base T (2/3);  total (15/28)
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.53e-01

ACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGAT  >  NC_002947/4033157‑4033440                                                                                                                                                |                                                                                                                                                aCAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTAcc                                                                                                                                         <  1:88067/149‑1 (MQ=255)    ggCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGa                                                                                                                                      <  5:44685/149‑1 (MQ=255)     gcgcTGCCGGTGATCACCCGCTACCTCGCCGCCCCGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGAc                                                                                                                                     <  8:230138/149‑1 (MQ=255)        cTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCg                                                                                                                                  <  4:89045/149‑1 (MQ=255)         tGCCGGTGATCAACCGCTACCTCGACGCGCCGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTAc                                                                                                                                          >  1:150824/1‑140 (MQ=255)         tGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTAc                                                                                                                                          <  2:150824/140‑1 (MQ=255)                  tCAACCGCTACCTCGACGCCCCGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGAc                                                                                                                        >  7:89941/1‑149 (MQ=255)                    aaCCGCTCCCTCGACGCCCCGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACgg                                                                                                                      <  6:259534/149‑1 (MQ=255)                     aCCGCTACCTCGACGCCCCGCTGCTACTGGCGCGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACgg                                                                                                                      >  1:165194/1‑149 (MQ=255)                     aCCGCTACCTCGACGCCCCGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGGCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACgg                                                                                                                      <  2:165194/148‑1 (MQ=255)                     aCCGCTACCTCGACGCCCCGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCgca                                                                                                                                <  8:67243/138‑1 (MQ=255)                     aCCGCTACCTCGACGCCCCGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCgca                                                                                                                                >  7:67243/1‑138 (MQ=255)                      ccGCTACCTCGACGCCCCGCTGCTACTGGC‑CGATGCCAGCACAGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGgca                                                                                                                    <  1:127700/149‑1 (MQ=255)                                gACGCCCCGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTg                                                                                                            <  7:202576/147‑1 (MQ=255)                                     cccGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATc                                                                                                     <  6:240634/149‑1 (MQ=255)                                     cccGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATc                                                                                                     >  4:101245/1‑149 (MQ=255)                                           gctACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATc                                                                                                     <  3:237394/143‑1 (MQ=255)                                           gctACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATc                                                                                                     >  4:237394/1‑143 (MQ=255)                                              aCTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTcgc                                                                                            >  8:219883/1‑149 (MQ=255)                                               cTGGC‑CGATGCCATCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAAAAACCTGCCCAAAGTCCTGCCCT‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTcgc                                                                                            <  1:66669/149‑1 (MQ=255)                                               cTGGC‑CGAAGCCAGCACCGGCCGGATTCGCATCCGCGGCATCTACAACAACCGGGAAGTGGGGCGCCTGGGCAACAACCTGCCCAAAGTCCTG‑CCGTGTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTcgc                                                                                            <  1:70131/149‑1 (MQ=255)                                                  gC‑CGGTGCTGGCATCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGTCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCT‑GTCTACCTGACCCGCAGCAAGGACGGCAGCAAAGTGCTCAATCGCATCTcgccgcc                                                                                        >  8:116220/1‑149 (MQ=255)                                                  gC‑CGGTGCTGGCATCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCT‑GTCTACCTGACACGCAGCAAGGACGGCAGCACCGTGCTCAATCGCGACTcgccgcc                                                                                        >  8:112742/1‑149 (MQ=255)                                                       aTGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATTTCGCCGCCGccc                                                                                    <  3:101245/149‑1 (MQ=255)                                                                 cGGACGCATTCGCATCACCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCt                                                                          <  6:189421/149‑1 (MQ=255)                                                                       cATTCGCATCAGCGGCATCTACAACACCCGGTAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGcc                                                                    <  2:65497/149‑1 (MQ=255)                                                                            gCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACACCCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCACTCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTccc                                                               >  6:205702/1‑149 (MQ=255)                                                                                   cGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTa                                                        <  7:219883/149‑1 (MQ=255)                                                                                       aTCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCT‑GTCTACCTGACCCGCATCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAg                                                    <  7:116220/149‑1 (MQ=255)                                                                                       aTCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCT‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAg                                                    <  7:112742/149‑1 (MQ=255)                                                                                               cacCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCTGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGc                                            <  1:49157/149‑1 (MQ=255)                                                                                                 cccGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTg                                          <  8:89941/149‑1 (MQ=255)                                                                                                    gggAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGcc                                        >  8:198434/1‑148 (MQ=255)                                                                                                           ggggCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGc                                 <  7:198434/148‑1 (MQ=255)                                                                                                              gcgcCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGaa                             >  1:26072/1‑149 (MQ=255)                                                                                                                                 cccAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCa                        <  1:105636/135‑1 (MQ=255)                                                                                                                                    aaaGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTaa       >  5:181361/1‑149 (MQ=255)                                                                                                                                     aaGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGGAAGGGCGGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAAt      <  6:15901/149‑1 (MQ=255)                                                                                                                                     aaGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAAt      <  1:49682/149‑1 (MQ=255)                                                                                                                                       gTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATgg    <  6:181361/149‑1 (MQ=255)                                                                                                                                          cTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGAt  >  3:228313/1‑148 (MQ=255)                                                                                                                                          cTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGAt  <  2:32785/148‑1 (MQ=255)                                                                                                                                          cTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGAt  >  1:32785/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                ACAGGCGCTGCCGGTGATCAACCGCTACCTCGACGCCCCGCTGCTACTGGC‑CGATGCCAGCACCGGCCGCATTCGCATCAGCGGCATCTACAACACCCGGGAAGTGGGGCGCCTGGTCAACAACCTGCCCAAAGTCCTG‑CCG‑GTCTACCTGACCCGCAGCAAGGACGGCAGCACCGTGCTCAATCGCATCTCGCCGCCGCCCGGCAAGGGCTGACGCCCTCCCCGCCTTACAAGGTCATCGCTGCCAACCAGCCGAACGCCATCAACGGCAGGTTGTAATGGAT  >  NC_002947/4033157‑4033440

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