Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,006,468	T→G	14.5%	F1017V (TTC→GTC)	PP_2627 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,006,468	0	T	G	14.5%	101.3 / 9.5	41	F1017V (TTC→GTC)	PP_2627	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (17/18);  new base G (4/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GCTGACGTGTACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTT  >  NC_002947/3006321‑3006610                                                                                                                                                    |                                                                                                                                                   gCTGACGTGTACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGttt                                                                                                                                                    >  1:231447/1‑148 (MQ=255)      cGTGTACAGTTTTTTCCTTCTACGCCTGGCAGGCCGCTAGACCTTAGTGGGGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGGTGTCCttc                                                                                                                                              <  6:248209/149‑3 (MQ=255)                       aGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGcc                                                                                                                             >  6:335360/1‑149 (MQ=255)                        gCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCt                                                                                                                            <  7:81329/149‑1 (MQ=255)                            cGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTg                                                                                                                          >  2:71360/1‑147 (MQ=255)                            cGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCa                                                                                                                        <  2:231447/149‑1 (MQ=255)                            cGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCa                                                                                                                        <  5:210242/149‑1 (MQ=255)                                  gCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCg                                                                                                                   <  3:87054/148‑1 (MQ=255)                                         gCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCAc                                                                                                           >  3:270696/1‑149 (MQ=255)                                             ccccTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACg‑g                                                                                                        <  5:6435/147‑1 (MQ=255)                                                cTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCt                                                                                                                            <  6:185644/125‑1 (MQ=255)                                                cTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCt                                                                                                                            >  5:185644/1‑125 (MQ=255)                                                    gtggtgAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATcc‑‑c                                                                                              <  4:270696/148‑1 (MQ=255)                                                                cccGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGa                                                                                 >  8:8392/1‑149 (MQ=255)                                                                 ccGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGAc                                                                                <  5:335360/149‑1 (MQ=255)                                                                 ccGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGAc                                                                                >  4:148559/1‑149 (MQ=255)                                                                   gATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGt                                                                              >  7:56411/1‑149 (MQ=255)                                                                     tGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGa                                                                            <  8:56411/149‑1 (MQ=255)                                                                               gATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCt                                                                  >  5:63794/1‑149 (MQ=255)                                                                                 tGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCg                                                                <  4:167063/149‑1 (MQ=255)                                                                                        gAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACGCGGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTgg                                                          <  1:71360/149‑1 (MQ=255)                                                                                             gTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAg‑c                                                   <  3:148559/149‑1 (MQ=255)                                                                                                        cGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTTGAGCTGGCGGAG‑CGCATCCGCGa                                         >  6:25318/1‑148 (MQ=255)                                                                                                         gAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGTGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGaaa                                       >  7:285463/1‑149 (MQ=255)                                                                                                               gcgcCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGTGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACcttct                                 <  8:285463/149‑1 (MQ=255)                                                                                                                     gcagcaGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGAt                            >  7:104175/1‑148 (MQ=255)                                                                                                                           gTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTg                     <  7:8392/149‑1 (MQ=255)                                                                                                                            tttAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAg‑c                                                   >  8:270513/1‑118 (MQ=255)                                                                                                                            tttAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAg‑c                                                   <  7:270513/118‑1 (MQ=255)                                                                                                                                 tGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGc               <  8:50280/149‑1 (MQ=255)                                                                                                                                       gTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCCGACTGATTCGTGCCGAGGTGTGCGGCCGGCAGGAGAT‑GCGGCGCCGCTTCCGCGAACCCTTCTTCGATCACGGTGGAATCCTAAg           >  1:276040/1‑149 (MQ=255)                                                                                                                                          tctcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGtgtg       <  8:104175/148‑1 (MQ=255)                                                                                                                                           ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGa                                                               <  8:141181/92‑1 (MQ=38)                                                                                                                                           ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGa                                                               >  7:141181/1‑92 (MQ=38)                                                                                                                                           ctcAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGtgtg       >  6:83343/1‑147 (MQ=255)                                                                                                                                           ctcAAGGTTGTCCTTCAGGAAAACCTTGAGACCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGAGCCGACGAGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCCCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGtgtg       >  7:256601/1‑147 (MQ=255)                                                                                                                                            tcAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAg    >  1:133696/1‑149 (MQ=255)                                                                                                                                             cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCg                                                                <  7:56179/89‑1 (MQ=255)                                                                                                                                             cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCg                                                                >  8:56179/1‑89 (MQ=255)                                                                                                                                             cAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt   >  4:99122/1‑149 (MQ=255)                                                                                                                                             cAAGGTTGTCCTTCAGGAAAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGTCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGt   >  6:246671/1‑149 (MQ=255)                                                                                                                                                ggTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGtt  >  5:7966/1‑147 (MQ=255)                                                                                                                                                  ttttCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGa                                                                        >  8:30641/1‑76 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCTGACGTGTACAGTTTTTTCCAGCAACGCCTGGCAGGCCGCTACCCCTTTGTGGTGAAGGCACCCGATGCTTCTTTGGATGACTTCGAGGCGTTCTTCGGGCCGAAAGGGCGCCTGCAGCAGTTTAATGATCGGTATCTCAAGGTTTTCCTGAAGGACAACCTTGAGGCCTTGCAGTCCGCTCAGCACG‑GGCAATCC‑‑CTGATTCGTGCCGACGTGATCGAGCAGCTCGAGCTGGCGGAG‑CGCATCCGCGAAACCTTCTTCGATCAGCGTGGCAAGCTAAGTGTGCAGTT  >  NC_002947/3006321‑3006610

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