Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,855,446	T→A	8.3%	intergenic (+18/+114)	tcyP → / ← folA	sodium/cystine symporter/dihydrofolate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,855,446	0	T	A	8.3%	138.6 / 5.4	48	intergenic (+18/+114)	tcyP/folA	sodium/cystine symporter/dihydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base T (20/24);  new base A (2/2);  total (22/26)
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.32e-01

CATGGGCCGTACGGCGTTGAACGTGAACGGTTCGATGACGGCGGGCGTTGTCACCAGCCAGCTGTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAGGGCTGCAAAGCAGACCCCAACAATCTCAACCCTTATCCCAAACCTCGAAGTGATACGCCG  >  NC_002947/5855302‑5855594                                                                                                                                                 |                                                                                                                                                     cATGGGCCGTACGGCGTTGAACGTGAACGGTTCGATGACGGCGGGCGTTGTCACCAGCCAGCTGTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTgcgc                                                                                                                                                 <  2:1219400/150‑1 (MQ=255)    gggCCGTACGGCGTTTAACGTGAACGGTTCGATGACGGCGGGCGTTGTCACCAGCCAGCTGTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGcccc                                                                                                                                              <  1:1199414/150‑1 (MQ=255)                         gaacgGTTCGATGACGGCGGGCGTTGTCACCAGCCAGCTGTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTcgc                                                                                                                                        >  1:592672/1‑135 (MQ=255)                         gaacgGTTCGATGACGGCGGGCGTTGTCACCAGCCAGCTGTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTcgc                                                                                                                                        <  2:592672/135‑1 (MQ=255)                         gaacgGTTCGATGACGGCGGGCGTTGTCACCAGCCAGCTGTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCgt                                                                                                                         >  2:453696/1‑150 (MQ=255)                              gTTCGATGACGGCGGGCGTTGTCACCAGCCAGCTGTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCtt                                                                                                                    >  2:939385/1‑150 (MQ=255)                                 cGATGACGGCGGGCGTTGTCACCAGCCAGCTGTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCt                                                                                                                 <  1:429400/150‑1 (MQ=255)                                  gATGACGGCGGGCGTTGTCACCAGCCAGCTGTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTc                                                                                                                <  2:1177351/150‑1 (MQ=255)                                      aCGGCGGGCGTTGTCACCAGCCAGCTGTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTAAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCaaaa                                                                                                            >  1:84575/1‑150 (MQ=255)                                                    cACCATCCAGCTGTTGAAAGAGACTGCCCCGCCAGTCCTGGTGGTTGCTGATGTTGCTGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTgc                                                                                              <  2:852117/150‑1 (MQ=255)                                                     aCCAGCCAGCTGTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTgcg                                                                                             >  1:762677/1‑150 (MQ=255)                                                      ccagccagCTGTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTc                                                                                                                                          <  2:681438/104‑1 (MQ=255)                                                      ccagccagCTGTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTc                                                                                                                                          >  1:681438/1‑104 (MQ=255)                                                                gTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAAc                                                                                                           >  2:1128994/1‑125 (MQ=255)                                                                gTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAAc                                                                                                           <  1:1128994/125‑1 (MQ=255)                                                                 ttGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCgt                                                                                                                         <  2:1185049/110‑1 (MQ=255)                                                                 ttGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCgt                                                                                                                         >  1:1185049/1‑110 (MQ=255)                                                                 ttGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAgg                                                                                 >  1:474485/1‑150 (MQ=255)                                                                    aaaGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGc                                                                              <  2:919537/150‑1 (MQ=255)                                                                    aaaGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCAATGCGccccg                                                                                                                                             >  2:1075138/1‑86 (MQ=255)                                                                    aaaGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCAATGCGccccg                                                                                                                                             <  1:1075138/87‑2 (MQ=255)                                                                     aaGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGATCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCg                                                                             >  2:862480/1‑150 (MQ=255)                                                                            tGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTa                                                                                   >  1:154160/1‑137 (MQ=255)                                                                            tGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTa                                                                                   <  2:154160/137‑1 (MQ=255)                                                                              aCAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCAccc                                                                    <  2:474485/150‑1 (MQ=255)                                                                                  ggAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTg                                                                                                                   >  2:228129/1‑99 (MQ=255)                                                                                  ggAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTg                                                                                                                   <  1:228129/99‑1 (MQ=255)                                                                                  ggAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGa                                                                >  1:1110611/1‑150 (MQ=255)                                                                                   gAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGaa                                                               <  2:262376/150‑1 (MQ=255)                                                                                   gAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCAACAATTTACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGaa                                                               <  1:336970/150‑1 (MQ=255)                                                                                    aaGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAg                                                              <  2:98694/150‑1 (MQ=255)                                                                                    aaGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAg                                                              <  2:762677/150‑1 (MQ=255)                                                                                              cGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAGGGCTGCAAAg                                                    >  2:970595/1‑150 (MQ=255)                                                                                                gTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAGGGCTGCAAAGCa                                                  <  1:294565/150‑1 (MQ=255)                                                                                                   atgaACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGaa                                                               >  2:763309/1‑134 (MQ=255)                                                                                                   atgaACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGaa                                                               <  1:763309/134‑1 (MQ=255)                                                                                                             cGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAGGGCTGCAAAGCAGACCCCAACAAtc                                     >  2:315816/1‑150 (MQ=255)                                                                                                                      ccACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAGGGCTGCAAAGCAGACCCCAACAATCTCAAc                                <  2:1075960/146‑1 (MQ=255)                                                                                                                      ccACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAGGGCTGCAAAGCAGACCCCAACAATCTCAAc                                >  1:1075960/1‑146 (MQ=255)                                                                                                                      ccACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAGGGCTGCAAAGCAGACCCCAACAATCTCAACCCtt                            >  2:107607/1‑150 (MQ=255)                                                                                                                        acacCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAGGGCTGCAAAGCAGACCCCAACAATCTCAACCCTTAt                          >  2:885821/1‑150 (MQ=255)                                                                                                                                  gATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAGGGCTGCAAAGCAGACCCCAACAATCTCAACCCTTATCCCAAACCTc                <  1:453696/150‑1 (MQ=255)                                                                                                                                     gaaaggggCGCAATGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAgg                                                             >  2:828101/5‑102 (MQ=255)                                                                                                                                     gaaaggggCGCAATGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAgg                                                             <  1:828101/98‑1 (MQ=255)                                                                                                                                      gCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAGGGCTGCAAAGCAGACCCCAACAATCTCAACCCTTATCCCAAACCTCGAAg            <  1:315816/150‑1 (MQ=255)                                                                                                                                        cGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAGGGCTGCAAAGCAGACCCCAACAATCTCAACCCTTATCCCAAACCTCGAAGTg          >  1:429656/1‑150 (MQ=255)                                                                                                                                               cATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAGGGCTGCAAAGCAGACCCCAACAATCTCAACCCTTATCCCAAACCTCGAAGTGATACGcc   <  1:220432/150‑1 (MQ=255)                                                                                                                                                aTTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAGGGCTGCAAAGCAGACCCCAACAATCTCAACCCTTATCCCAAACCTCGAAGTGATACGCCg  <  1:862480/150‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CATGGGCCGTACGGCGTTGAACGTGAACGGTTCGATGACGGCGGGCGTTGTCACCAGCCAGCTGTTGAAAGAGACTGACAAGGAAGTGCTGGCCGGTGATGAACATGCCGAGCTGAGCCACACCTGACAGATCGCCGGGGCGCATTGCGCCCCTTTCGCGACAAGGCCGCTCGTGTCTTGCTCAAAACCTGAACACCCTGCGCCCCCCCGTAGGAGCGGGTTCACCCGCGAAGGGCTGCAAAGCAGACCCCAACAATCTCAACCCTTATCCCAAACCTCGAAGTGATACGCCG  >  NC_002947/5855302‑5855594

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑