Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,953,664:1	+T	9.5%	coding (834/1383 nt)	PP_2586 ←	putative amino acid transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,953,664	1	.	T	9.5%	153.2 / 10.4	42	coding (834/1383 nt)	PP_2586	putative amino acid transporter
Reads supporting (aligned to +/- strand):  ref base . (17/21);  new base T (2/2);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

CGCGGCCCATGGCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGTTTCCTCAGCCATGGTCGACACC  >  NC_002947/2953516‑2953811                                                                                                                                                      |                                                                                                                                                    cgcgGCCCCTGGCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGa                                                                                                                                                      <  4:57799/149‑1 (MQ=255)      cccATGGCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGc                                                                                                                                                <  7:220567/149‑1 (MQ=255)       ccATGGCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCg                                                                                                                                               <  8:8591/149‑1 (MQ=255)        cATGGCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCg                                                                                                                                               >  7:275297/1‑148 (MQ=255)            gCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCG‑ATCCGGATTCGGCgtcgg                                                                                                                                           <  6:126982/148‑5 (MQ=255)                           gACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGaa                                                                                                                           <  8:189657/149‑1 (MQ=255)                           gACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGaa                                                                                                                           <  3:66395/149‑1 (MQ=255)                              acGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAacac                                                                                                                        <  8:275297/149‑1 (MQ=255)                                gCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAAACGGGGAccacca                                                                                                                      >  6:166826/1‑149 (MQ=255)                                 cTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCAt                                                                                                                     <  5:134254/149‑1 (MQ=255)                                    gCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATc                                                                                                                    <  4:114461/147‑1 (MQ=255)                                         ggACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTgc                                                                                                               <  5:244115/147‑1 (MQ=255)                                             acCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCAGCCCCAATCGGGGAACACCATCTGTGCGCAAt                                                                                                          >  8:163148/1‑148 (MQ=255)                                               cATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAAt                                                                                                          >  5:225861/1‑146 (MQ=255)                                                         gCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGAcaccac                                                                                                >  6:197039/1‑146 (MQ=255)                                                           gAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGa                                                                                            >  6:155564/1‑148 (MQ=255)                                                                agcCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCACCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGccc                                                                                                                                      >  4:51190/1‑101 (MQ=255)                                                                agcCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCACCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGccc                                                                                                                                      <  3:51190/101‑1 (MQ=255)                                                                   ccGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAgg                                                                                   >  4:308772/1‑149 (MQ=255)                                                                               tagCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAAc                                                                       <  6:225861/149‑1 (MQ=255)                                                                                       aaGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGc                                                               >  5:267072/1‑149 (MQ=255)                                                                                       aaGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGc                                                               <  6:267072/149‑1 (MQ=255)                                                                                          gCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGc                                                               >  6:226313/1‑146 (MQ=255)                                                                                              ttACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAAc                                                                       >  8:184393/1‑134 (MQ=255)                                                                                              ttACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAAc                                                                       <  7:184393/134‑1 (MQ=255)                                                                                                    gcagcTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACg                                                  >  5:226533/1‑149 (MQ=255)                                                                                                            ccgccCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCg                                          <  5:197039/149‑1 (MQ=255)                                                                                                                  aCCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGt                                    >  4:259938/1‑149 (MQ=255)                                                                                                                  aCCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGt                                    >  1:270328/1‑149 (MQ=255)                                                                                                                     cGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCgg                                 <  6:226533/149‑1 (MQ=255)                                                                                                                        cgcATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGtt                              <  7:214039/149‑1 (MQ=255)                                                                                                                            tCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTg                          <  3:308772/149‑1 (MQ=255)                                                                                                                            tCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTg                          <  5:226313/149‑1 (MQ=255)                                                                                                                               cGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATGGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGtt                       >  5:255438/1‑149 (MQ=255)                                                                                                                                     ccGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGTTTCCTCa                 <  5:155564/149‑1 (MQ=255)                                                                                                                                       ttcGCCG‑ATCCGGATTCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTAc                            <  8:249111/134‑1 (MQ=255)                                                                                                                                       ttcGCCG‑ATCCGGATTCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTAc                            >  7:249111/3‑136 (MQ=255)                                                                                                                                       ttcGCCG‑ATCCGGATTCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGAACAACGACACCGCCAGCATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGTTTCATCAGc               >  1:363947/3‑149 (MQ=255)                                                                                                                                        aCGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCgg                                         <  8:132127/122‑1 (MQ=255)                                                                                                                                        aCGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCgg                                         >  7:132127/1‑122 (MQ=255)                                                                                                                                           ccGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATGGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGTTTCCTCAGCCATg           <  6:255438/149‑1 (MQ=255)                                                                                                                                              aaTCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGTTTCCTCAGCCATGGTc        <  7:163148/149‑1 (MQ=255)                                                                                                                                                    gA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGTTTCCTCAGCCATGGTCGACAcc  >  5:274591/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CGCGGCCCATGGCGAACAGTACGCGCGACACGCTGGCCTGGGACACCATGGCCGAGGCGAAGCAGCCCGCCACGTAGGTAGCGGTGAAGGCGGTTACCAGCAGCTCGCCGCCCACCCGGCGCATCACGTCAACCGACGCCGAATCCGGA‑TCGGCGAAGCTGCCCCAATCGGGGAACACCATCTGTGCGCAATAAGACACCACCAGAAATAGCAGGCCGCCGATCAACGACACCGCCAGAATCGCCACGGGAATGCGGTAGGTCGGGTTACTGGTTTCCTCAGCCATGGTCGACACC  >  NC_002947/2953516‑2953811

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