Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,042,172	C→T	10.5%	A25A (GCC→GCT)	PP_5537 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,042,172	0	C	T	10.5%	100.5 / 3.1	39	A25A (GCC→GCT)	PP_5537	membrane protein
Reads supporting (aligned to +/- strand):  ref base C (17/17);  new base T (2/2);  total (20/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 = 9.06e-01

CATGACCTCACGGGTCATTGAGTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGAGCCGGGATCGCGGTGC  >  NC_002947/3042024‑3042320                                                                                                                                                     |                                                                                                                                                     cATGACCTCACGGGTCATTGAGTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGcc                                                                                                                                                      <  2:164338/149‑1 (MQ=255)    gACCTCACGGGTCATTGAGTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAg                                                                                                                                                   >  8:137546/1‑149 (MQ=255)             ggTCATTGAGTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCgctg                                                                                                                                          <  2:237893/149‑1 (MQ=255)                     aGTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTACGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCTCGGCAATcagcat                                                                                                                                         <  1:305738/142‑6 (MQ=255)                     aGTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTACGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCTCGGCAATcagcat                                                                                                                                         >  2:305738/1‑137 (MQ=255)                      gTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTc                                                                                                                                 <  6:67898/149‑1 (MQ=255)                          aagcaaCCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTc                                                                                                                                 <  4:224865/145‑1 (MQ=255)                          aagcaaCCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTc                                                                                                                                 >  3:224865/1‑145 (MQ=255)                            gcaaCCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTg                                                                                                                           <  5:312991/149‑1 (MQ=255)                             caaCCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTgg                                                                                                                          <  8:165180/149‑1 (MQ=255)                                  cTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCAct                                                                                                                      >  3:96743/1‑148 (MQ=255)                                   tATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGATTGCGTATCTGCGCGCTGGCAATACTGGTGATGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCAGAGCAATCGCTGCTGGCATTCGGGCTGGCActgc                                                                                                                    >  1:149869/1‑149 (MQ=255)                                    aTCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCActgct                                                                                                                   >  2:322283/1‑149 (MQ=255)                                     tCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCActgc                                                                                                                    <  6:117711/147‑1 (MQ=255)                                       actactGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGgc                                                                                                                >  5:146430/1‑149 (MQ=255)                                           ctGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCa                                                                                                             <  5:43644/148‑1 (MQ=255)                                              tgCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGc                                                                                                          <  6:146430/148‑1 (MQ=255)                                                   gCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGacac                                                                                                    <  3:139904/149‑1 (MQ=255)                                                     aaaCCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACgg                                                                                                  <  1:322283/149‑1 (MQ=255)                                                         ccTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCTCGGCAATca                                                                                                                                             >  8:175098/1‑101 (MQ=255)                                                         ccTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCTCGGCAATca                                                                                                                                             <  7:175098/102‑2 (MQ=255)                                                           tATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAgg                                                                                            >  3:366477/1‑149 (MQ=255)                                                                 ggAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCa                                                                                      >  6:20244/1‑149 (MQ=255)                                                                    gagCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATc                                                                                    <  4:96743/148‑1 (MQ=255)                                                                      gCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGAc                                                                                 <  7:133490/149‑1 (MQ=255)                                                                       cTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGAcc                                                                                <  7:137546/149‑1 (MQ=255)                                                                             gggCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTAcc                                                                          >  5:212422/1‑149 (MQ=255)                                                                                       aTCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCgg                                                                >  3:42988/1‑149 (MQ=255)                                                                                           gcgcgcTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTg                                                            >  7:241796/1‑149 (MQ=255)                                                                                               gcTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACCCGGCCCAGGTGGTCATAGACCTGTACCTGATGGCCGGGCTGGCCg                                                        >  8:333192/1‑149 (MQ=255)                                                                                                     aaTACTTTTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGATGGCGCAGGCGGACACGGCCAAGCTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCt                                                  >  4:243801/1‑149 (MQ=255)                                                                                                          tGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACg                                                                                                   <  7:89292/95‑1 (MQ=255)                                                                                                          tGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACg                                                                                                   >  8:89292/1‑95 (MQ=255)                                                                                                            gTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGt                                           <  4:366477/149‑1 (MQ=255)                                                                                                                      cGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACaa                                  <  5:20244/148‑1 (MQ=255)                                                                                                                                           tGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGCCCCAGGTGGTCATCGACCTGTACCTGATCGCCGGGCTGGCCGGTTGCTGGATGGTCAGAGACAAACGCATGCGCGGCAGAGACgg             >  6:272443/1‑148 (MQ=255)                                                                                                                                               gATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGCACCGCGATcc        >  2:162272/1‑148 (MQ=255)                                                                                                                                               gATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGAGCCGGGATcg        >  1:169297/1‑149 (MQ=255)                                                                                                                                                 ttGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGAGCCGGGATcgcg      <  4:42988/149‑1 (MQ=255)                                                                                                                                                  tGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGAGCCGGGATCGCgg     >  3:365925/1‑149 (MQ=255)                                                                                                                                                     cGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGc                                                   >  4:262165/1‑100 (MQ=255)                                                                                                                                                     cGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGc                                                   <  3:262165/100‑1 (MQ=255)                                                                                                                                                     cGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGAGCCGGGATCGCGGtgc  >  4:6622/1‑149 (MQ=255)                                                                                                                                                     cGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGAGCCGGGATCGCGGtgc  >  3:338577/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CATGACCTCACGGGTCATTGAGTGCAAGCAACCCTATCACTACTGTGCCGGCAAACCCTATCAAGGAGAGCTATATGGGCTTGCGTATCTGCGCGCTGGCAATACTGGTGCTGTTTTCGCTGTATACCGGCGGGACCATGCTGATTGCCGAGCAATCGCTGCTGGCATTCGGGCTGGCACTGCTGGCGCAGGCGGACACGGCCCAGGTGGTCATCGACCTGTACCTGATGGCCGGGCTGGCCGGTTGCTGGATGGTGCGAGACAACCGCATGCGCGGCAGAGCCGGGATCGCGGTGC  >  NC_002947/3042024‑3042320

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