Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	991,554	G→A	8.9%	G59D (GGT→GAT)	der →	GTPase Der

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	991,554	0	G	A	8.9%	124.1 / 2.5	45	G59D (GGT→GAT)	der	GTPase Der
Reads supporting (aligned to +/- strand):  ref base G (18/23);  new base A (2/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

AACGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCG  >  NC_002947/991412‑991698                                                                                                                                                |                                                                                                                                                 aaCGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCAc                                                                                                                                            >  8:339315/1‑149 (MQ=255)   cGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCg                                                                                                                                          <  3:132516/149‑1 (MQ=255)       ggCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGa                                                                                                                                      <  1:320080/149‑1 (MQ=255)       ggCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGATACCACCGGTgt                                                                                                                                      <  5:105535/149‑2 (MQ=255)               caccaTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTgg                                                                                                                              <  7:131649/149‑1 (MQ=255)                accaTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTggg                                                                                                                             >  2:340596/1‑149 (MQ=255)                accaTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTggg                                                                                                                             <  6:61767/149‑1 (MQ=255)                   aTGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCAt                                                                                                                          <  2:266609/149‑1 (MQ=255)                         aaCCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCCCCGGTGACGAAGTGGGCATGGACg                                                                                                                     >  3:40008/1‑148 (MQ=255)                            cGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGaaa                                                                                                                  >  6:63485/1‑148 (MQ=255)                               cTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAAt                                                                                                              <  1:340596/149‑1 (MQ=255)                                  aCCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGc                                                                                                           <  7:326753/149‑1 (MQ=255)                                  aCCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGc                                                                                                           <  7:339315/149‑1 (MQ=255)                                    cAAGTCCCGCGATGCCATCGTGGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTTATGCGAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGaaaaaa                                                                                                               <  2:160243/143‑1 (MQ=255)                                    cAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGcc                                                                                                          >  6:372625/1‑148 (MQ=255)                                    cAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACAGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGTTTACGAAGTGGGCATGGACGaaaaaa                                                                                                               >  1:160243/1‑143 (MQ=255)                                       gACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAg                                                                                                       >  6:213309/1‑148 (MQ=255)                                       gACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAg                                                                                                       >  8:183754/1‑148 (MQ=255)                                       gACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGc                                                                                                      <  5:213309/149‑1 (MQ=255)                                        aCCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAg                                                                                                       >  6:49733/1‑147 (MQ=255)                                         cccGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGa                                                                                                                    <  7:314851/133‑1 (MQ=255)                                         cccGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGa                                                                                                                    >  8:314851/1‑133 (MQ=255)                                          ccGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCAGTGACCGCCAGTATGGTGATGCCAGCTAGAAGGGTCGTTCCTTCCTCCTGA‑TCGAAACAGGTGGTATCACCCGTGAAGAAGAGCCCATGGAAGAAAAAATGGACGAGCAGt                                                                                                   >  8:236534/1‑149 (MQ=255)                                            gcgATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCg                                                                                                 >  8:57869/1‑149 (MQ=255)                                                        tcgGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATGCTGA‑TCGACACCGGTTGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCAtt                                                                                     <  7:236534/149‑1 (MQ=255)                                                          ggTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTg                                                                                    <  5:63485/148‑1 (MQ=255)                                                           gTGACCTGTCCGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCACATGGCCATTgaa                                                                                  <  8:323079/149‑1 (MQ=255)                                                                    cGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTACAACCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGAc                                                                         <  4:3893/149‑1 (MQ=255)                                                                         cTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTAtgt                                                                    >  4:258809/1‑149 (MQ=255)                                                                                    ccGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTgg                                                         <  2:80030/149‑1 (MQ=255)                                                                                        cAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGa                                                     >  5:17559/1‑149 (MQ=255)                                                                                                   tGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCgg                                          <  5:372625/149‑1 (MQ=255)                                                                                                       aGCTGGCAGGGTCGTTCCTTCATCCTGAGTCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTa                                       >  1:204723/1‑149 (MQ=255)                                                                                                           ggCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACt                                  >  4:294858/1‑149 (MQ=255)                                                                                                                    gTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGAc                         <  7:57869/149‑1 (MQ=255)                                                                                                                    gTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGAc                         <  5:49733/149‑1 (MQ=255)                                                                                                                        cTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGAc                         >  8:162687/1‑145 (MQ=255)                                                                                                                        cTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAga                     <  7:183754/149‑1 (MQ=255)                                                                                                                            aTCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAgatgat                 >  6:16477/1‑149 (MQ=255)                                                                                                                             tCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATc                <  3:294858/149‑1 (MQ=255)                                                                                                                                  a‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGa           <  3:258809/149‑1 (MQ=255)                                                                                                                                   ttCGTCACCGGTGATACCACCGGTGTCGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTa                                       <  7:331633/120‑1 (MQ=255)                                                                                                                                   ttCGTCACCGGTGATACCACCGGTGTCGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTa                                       >  8:331633/2‑121 (MQ=255)                                                                                                                                       tcacCGGTGATACCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCAc       >  3:152831/2‑149 (MQ=255)                                                                                                                                         acCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCt     <  6:17559/149‑1 (MQ=255)                                                                                                                                            ggtggtATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCATGGTTGATGCCAGTGCGGGTATGACTGCGGCAGACAAGATGATCGCCGAGCACCTgcg  >  7:83992/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                 AACGTCGGCAAATCCACCATGTTCAACCGCCTGACCAAGACCCGCGATGCCATCGTCGGTGACCTGTCGGGCCTGACCCGTGACCGCCAGTATGGTGATGCCAGCTGGCAGGGTCGTTCCTTCATCCTGA‑TCGACACCGGTGGTATCACCGGTGACGAAGTGGGCATGGACGAAAAAATGGCCGAGCAGTCGCTCATGGCCATTGAAGAAGCCGACTATGTGCTGTTCCTGGTTGATGCCCGTGCGGGTATGACTGCGGCAGACCAGATGATCGCCGAGCACCTGCG  >  NC_002947/991412‑991698

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