Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,840,132	G→A	8.7%	R352H (CGC→CAC)	PP_3389 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,840,132	0	G	A	8.7%	119.8 / 2.6	47	R352H (CGC→CAC)	PP_3389	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (24/18);  new base A (2/2);  total (27/20)
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.12e-01

GCTCAGCCTGCTGCTGTGCTGGCGACCGCGCTGGCCGAAGGTCATCCGCGAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCT  >  NC_002947/3839989‑3840275                                                                                                                                                |                                                                                                                                                gCTCAGCCTGCTGCTGTGCTGGCGACCGCGCTGGCCGAAGGTCATCCGCGAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGCCACGCCAGCTCGCCAACGCCGACCCCTgcc                                                                                                                                            >  8:211412/1‑147 (MQ=12)           ctgctgTGCTGGCGACCGCGCTGGCCGAAGGTCATCCGCGAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAg                                                                                                                                    <  7:198844/147‑1 (MQ=32)              ctgtgctgGCGACCGCGCTGGCCGAAGGTCATCCGCGAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCgg                                                                                                                               >  7:43808/1‑149 (MQ=32)                   ctgGCGACCGCGCTGGCCGAAGGTCATCCGCGAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGc                                                                                                                          >  4:171595/1‑149 (MQ=32)                     ggCGACCGCGCTGGCCGAAGGTCATCCGCGAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGcgg                                                                                                                        >  7:222206/1‑149 (MQ=32)                      gCGACCGCGCTGGCCGAAGGTCATCCGCGAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGc                                                                                                                                                <  4:188716/124‑1 (MQ=32)                      gCGACCGCGCTGGCCGAAGGTCATCCGCGAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGc                                                                                                                                                >  3:188716/1‑124 (MQ=32)                      gCGACCGCGCTGGCCGAAGGTCATCCGCGAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGcggc                                                                                                                       >  6:57207/1‑149 (MQ=32)                      gCGACCGCGCTGGCCGAAGGTCATCCGCGAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGATGGCAGCTCAAGCCCGGCAAGcggc                                                                                                                       >  6:73797/1‑149 (MQ=17)                        gACCGCGCTGGCCGAAGGTCATCCGCGAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGcggccg                                                                                                                     >  3:220393/1‑149 (MQ=32)                            gcgcTGGCCGAAGGTCATCCGCGAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGa                                                                                                                 <  8:43808/149‑1 (MQ=32)                                               cgcgAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAAcgcg                                                                                              <  8:222206/149‑1 (MQ=39)                                                  gAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAg                                                                                                                           >  2:100688/1‑117 (MQ=32)                                                    ggccggccTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAAcgcg                                                                                              >  3:247592/1‑144 (MQ=39)                                                                   cAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGa                                                                          <  3:171595/149‑1 (MQ=255)                                                                       ccAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAAcc                                                                      >  2:85967/1‑149 (MQ=255)                                                                           tCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCt                                                                  >  4:168921/1‑149 (MQ=255)                                                                                     gACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCACTGGCAGCTCAAGCCCGGCAAGCGTCCGGCGAAAAAGCCACGCAAAACgcgc                                                                                             <  6:178993/112‑1 (MQ=255)                                                                                     gACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCACTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACgcgc                                                                                             <  6:175068/112‑1 (MQ=255)                                                                                     gACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCACTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACgcgc                                                                                             >  5:178993/1‑112 (MQ=255)                                                                                     gACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCACTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACgcgc                                                                                             >  5:175068/1‑112 (MQ=255)                                                                                       cAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACg                                                        >  2:89327/1‑147 (MQ=255)                                                                                       cAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGaa                                                      <  7:211412/149‑1 (MQ=255)                                                                                          cccTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCaaa                                                                                    <  3:168428/116‑1 (MQ=255)                                                                                          cccTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCaaa                                                                                    >  4:168428/1‑116 (MQ=255)                                                                                           ccTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGAtt                                                  >  4:218457/1‑149 (MQ=255)                                                                                            cTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATtg                                                 >  1:265181/1‑149 (MQ=255)                                                                                             tGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATtgt                                                >  2:155084/1‑149 (MQ=255)                                                                                               aGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTgtgt                                              <  3:168921/149‑1 (MQ=255)                                                                                                 cGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCCGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCt                                            <  1:85967/149‑1 (MQ=255)                                                                                                        ccTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATtg                                                 >  5:258142/1‑137 (MQ=255)                                                                                                        ccTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATtg                                                 <  6:258142/137‑1 (MQ=255)                                                                                                        ccTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTgtgt                                              <  8:87001/140‑1 (MQ=255)                                                                                                        ccTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTgtgt                                              >  7:87001/1‑140 (MQ=255)                                                                                                          tCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAg                                    >  8:120219/1‑148 (MQ=255)                                                                                                          tCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACCCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAg                                    >  4:149073/1‑148 (MQ=255)                                                                                                                    cAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGACAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGACCCAAAGCGCAATAGCAACCGGCTGCGCCCAACGAAGATTCTGTCTGAATTCAGAATTCGAACgt                         >  8:39632/1‑149 (MQ=25)                                                                                                                     aGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTc                                             >  2:189610/1‑128 (MQ=255)                                                                                                                     aGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTc                                             <  1:189610/128‑1 (MQ=255)                                                                                                                       gCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTa                      <  3:218457/149‑1 (MQ=255)                                                                                                                                 tCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGAt             >  3:51372/1‑148 (MQ=255)                                                                                                                                       aCGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCgg      <  5:73797/149‑1 (MQ=255)                                                                                                                                       aCGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCgg      <  5:57207/149‑1 (MQ=255)                                                                                                                                          ccGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTtc   <  3:149073/149‑1 (MQ=255)                                                                                                                                          ccGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTtc   <  1:89327/149‑1 (MQ=255)                                                                                                                                          ccGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTtc   <  4:247592/149‑1 (MQ=255)                                                                                                                                            gACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTtct  >  1:67605/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                GCTCAGCCTGCTGCTGTGCTGGCGACCGCGCTGGCCGAAGGTCATCCGCGAGGCCGGCCTGCCCGGCAAGCCAGTCCCACCCGTGACAACCCTGAGCGACTACCCTCCCGCAATCCAGGCACGCCTGCTCGCCAACGCCGACCGCTGGCAGCTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCT  >  NC_002947/3839989‑3840275

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