Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,266,113	A→T	11.0%	intergenic (+30/+29)	PP_2863 → / ← PP_2864	membrane protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,266,113	0	A	T	11.0%	120.4 / 5.1	45	intergenic (+30/+29)	PP_2863/PP_2864	membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/17);  new base T (2/3);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATG  >  NC_002947/3265978‑3266256                                                                                                                                          |                                                                                                                                                ccGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTTCCGGCCCTTTCGCGGCTTTACCCGCGaagag                                                                                                                                    <  5:215029/149‑4 (MQ=38) ccGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTTCCGGCCCTTTCGCGGCTTTACCCGCGaagag                                                                                                                                    >  6:215029/1‑146 (MQ=38)   gCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTTCCGGCCCTTTCGCGGCTTTACCCGCGaagaggc                                                                                                                                  <  3:13874/149‑6 (MQ=38)      aGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCgg                                                                                                                               <  5:241236/149‑1 (MQ=255)        gTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGaa                                                                                                                             <  7:186843/149‑1 (MQ=255)         tttCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGaa                                                                                                                             >  5:115203/1‑148 (MQ=255)              ttCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTc                                                                                                                       >  4:32064/1‑149 (MQ=255)              ttCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTc                                                                                                                       >  6:30557/1‑149 (MQ=255)              ttCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTc                                                                                                                       >  5:133874/1‑149 (MQ=255)                       gtgtGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTCCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATg                                                                                                              <  1:85437/149‑1 (MQ=255)                                    ttGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGTCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCt                                                                                                 <  4:128180/149‑1 (MQ=255)                                       tgctgctgcACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCCGTTTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCAGGAACAGGTCGATCACATGACCAACCCCACCTg                                                                                                >  1:77264/1‑149 (MQ=255)                                          tgctgcACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCaa                                                                                             <  1:214610/147‑1 (MQ=255)                                          tgctgcACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCaa                                                                                             >  2:214610/1‑147 (MQ=255)                                            ctgcACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCg                                                                                          <  3:32064/148‑1 (MQ=255)                                            ctgcACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCg                                                                                          <  5:30557/148‑1 (MQ=255)                                                aCCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCtt                                                                                       >  3:172137/1‑147 (MQ=255)                                                aCCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCtt                                                                                       <  4:172137/147‑1 (MQ=255)                                                     gtgtTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAAc                                                                                >  6:275821/1‑149 (MQ=255)                                                     gtgtTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAAc                                                                                >  6:138532/1‑149 (MQ=255)                                                        ttCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGa                                                                             >  6:204627/1‑149 (MQ=255)                                                              gCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCAc                                                                        >  1:248024/1‑148 (MQ=255)                                                                  aCTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCa                                                                   <  6:133874/149‑1 (MQ=255)                                                                    tGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAg                                                                  >  4:204654/1‑148 (MQ=255)                                                                      gCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAg                                                                  <  3:204654/146‑1 (MQ=255)                                                                      gCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGc                                                                 <  2:77264/147‑1 (MQ=255)                                                                       cTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTgg                                                              >  1:117079/1‑149 (MQ=255)                                                                                      cGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAac                                               >  4:4925/1‑149 (MQ=255)                                                                                             tGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTca                                         >  4:164774/1‑148 (MQ=255)                                                                                              gCCCC‑‑TTGCTTGAAGTGCAGTTGTACCGGCCTCTACGCCCGTAAAGCCGCGAAAGCCCCGCAACACGTCGATCACAAGACCAACCCCACCACCCACCCCATCGGCAACCGACGCACTACCACCAGCACCGACCGCTCCAACACCTcaca                                       >  6:115186/1‑149 (MQ=255)                                                                                                  c‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGt                                   >  4:60229/1‑149 (MQ=255)                                                                                                     ttGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGc                                                                                         <  7:96352/94‑1 (MQ=255)                                                                                                     ttGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGc                                                                                         >  8:96352/1‑94 (MQ=255)                                                                                                         ttGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTtcatc                              <  5:204627/149‑1 (MQ=255)                                                                                                               tgCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGccc                        >  8:105668/1‑149 (MQ=255)                                                                                                                   gTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCg                                                      <  7:213075/115‑1 (MQ=255)                                                                                                                   gTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCg                                                      >  8:213075/1‑115 (MQ=255)                                                                                                                    ttGTGCCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgg                   >  8:161594/1‑149 (MQ=255)                                                                                                                    ctatctcGGCCCTTTCGCGGCTTTACCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATgg                   <  3:243763/143‑1 (MQ=255)                                                                                                                      gTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCAGACAGTTCATCACGCCCCATggg                  <  5:115186/148‑1 (MQ=255)                                                                                                                         ccGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTa               <  2:117079/148‑1 (MQ=255)                                                                                                                          cGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAgg             >  3:159482/1‑149 (MQ=255)                                                                                                                            gccctttCGCGGCTTTACCCGCGAAGAGGCCGGTACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGCCGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGggc           >  7:183744/6‑149 (MQ=38)                                                                                                                              ctctTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGggcgg         >  3:226853/1‑149 (MQ=255)                                                                                                                                     cgGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATg  <  5:138532/149‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                CCGCAAGGTTTCGTTCCAGGTGGTGTGCTGGGGGATTGTGCTGCTGCACCAGGTGTTCTGGGCAGACTGGCTGTTTCTTGACGGGCGTTATCTGCCCC‑‑TTGGTTGATGTGCAGTTGTACCGGCCTCTTCGCGGGTAAAGCCGCGAAAGGGCCGGAACAGGTCGATCACATGACCAACCCCACCTGCAACCGCTTGGGCAACCGACGCACTACCAGCTGGTGCGACCGCTGCAACACCTCACACAGTTCATCACGCCCCATGGGGTAGGGCGGCGCCATG  >  NC_002947/3265978‑3266256

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