Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,008,313	A→G	11.3%	M261T (ATG→ACG)	PP_2628 ←	ABC transporter ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,008,313	0	A	G	11.3%	136.0 / 7.1	53	M261T (ATG→ACG)	PP_2628	ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (24/23);  new base G (3/3);  total (27/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 = 1.00e+00

CGAACTGGCTCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCA‑GTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGA  >  NC_002947/3008167‑3008460                                                                                                                                                    |                                                                                                                                                     cGAACTGGCTCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATc                                                                                                                                                    >  1:298562/1‑149 (MQ=255)   aaCTGGCTCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATcga                                                                                                                                                  <  4:318963/149‑1 (MQ=255)          tCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCgg                                                                                                                                           >  1:344639/1‑149 (MQ=255)            aCCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCggg                                                                                                                                          >  8:169267/1‑148 (MQ=255)            aCCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCggg                                                                                                                                          >  4:88757/1‑148 (MQ=255)            aCCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCgg                                                                                                                                           <  2:344639/147‑1 (MQ=255)             ccAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCa                                                                                                                                        >  4:309219/1‑149 (MQ=255)              cAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCGTCGATGACCGGGcgg                                                                                                                                       <  8:254434/149‑3 (MQ=255)              cAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCGTCGATGACCGGGcgg                                                                                                                                       <  7:119821/149‑3 (MQ=255)               aGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCAtt                                                                                                                                      <  4:175189/149‑1 (MQ=255)               aGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCAtt                                                                                                                                      <  6:314748/149‑1 (MQ=255)                ggCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTg                                                                                                                                     <  2:81427/149‑1 (MQ=255)                ggCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTg                                                                                                                                     <  3:115958/149‑1 (MQ=255)                  cgcATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGc                                                                                                                                   <  3:309219/149‑1 (MQ=255)                                      cTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTAGCCTGCTCCTGCCCGACCTTGc                                                                                                               <  7:237125/149‑1 (MQ=255)                                         cTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCa                                                                                                            >  2:111118/1‑149 (MQ=255)                                          ttGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAg                                                                                                           >  2:51645/1‑149 (MQ=255)                                          ttGAAATCTTCCATGGCTGCGCGCGCAAAGTGAGCTTCGCGGTTGGAGTGCGAAAACAGATTTACCCTGGTGATCTTGGTGTAAGCTTCAGAAACCCGCCCCG‑TCATCGACGACAGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCac                                                                                                           >  4:45985/1‑148 (MQ=255)                                            gAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGc                                                                                                                                   <  8:149096/123‑1 (MQ=255)                                            gAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGc                                                                                                                                   >  7:149096/1‑123 (MQ=255)                                              aaTCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCgg                                                                                                       <  4:359275/149‑1 (MQ=255)                                               aTCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCggg                                                                                                      <  4:258056/149‑1 (MQ=255)                                                  ttCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCTGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCgggg                                                                                                     >  8:17636/1‑147 (MQ=255)                                                  ttCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCTGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCgggg                                                                                                     <  7:17636/147‑1 (MQ=255)                                                    ccATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCt                                                                                                                           >  4:156385/1‑123 (MQ=255)                                                    ccATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCt                                                                                                                           <  3:156385/123‑1 (MQ=255)                                                      aTGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGaa                                                                                               <  1:51645/149‑1 (MQ=255)                                                             cgcgcgcAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGGCACATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACa                                                                                          <  1:111118/148‑1 (MQ=255)                                                                        tGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGACCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGcc                                                                             >  6:314708/1‑149 (MQ=255)                                                                            cTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGGGTAAGCGTCAGAAACCCGCCCAG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACCGCATGGCCAAGCCGAAc                                                                         <  8:201768/149‑1 (MQ=255)                                                                                            cGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCaa                                                         >  2:83381/1‑149 (MQ=255)                                                                                                 aCAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCTCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAg                                                    <  7:169267/149‑1 (MQ=255)                                                                                                 aCAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAg                                                    <  5:314708/149‑1 (MQ=255)                                                                                                     ttttACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAg                                                             >  4:96665/1‑136 (MQ=255)                                                                                                     ttttACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAg                                                             <  3:96665/136‑1 (MQ=255)                                                                                                       ttACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGgcag                                               <  2:298562/148‑1 (MQ=255)                                                                                                       ttACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGgcag                                               >  7:278809/1‑148 (MQ=255)                                                                                                         aCCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGgcagca                                             <  1:83381/148‑1 (MQ=255)                                                                                                                    gTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGACCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCA‑GTTTAAggg                                >  6:30299/1‑149 (MQ=255)                                                                                                                     ttGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCA‑GTTTCAggg                                >  8:158835/1‑148 (MQ=255)                                                                                                                     ttGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATAGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGTCAGCATCA‑GTTTCAggg                                >  5:172658/1‑148 (MQ=255)                                                                                                                           tAAGAGTCAGAAACCCGCCCGG‑TCATCGACGACAGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGCAATACCGCATGGCCAAGCCGAACAGCACCACCCAGGCAATCAAGGGCACCATCC‑GTTTCAGGGCAAAgcc                         >  8:153867/1‑149 (MQ=255)                                                                                                                                 tCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGACTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGCGCCGCATCA‑GTTTCAGGGCAAAGCCGCCGGc                   >  2:215212/1‑149 (MQ=255)                                                                                                                                    gAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCA‑GTTTCAGGGCAAAGCCGCCGGCCAg                >  8:83099/1‑149 (MQ=255)                                                                                                                                    gAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCA‑GTTTCAGGGCAAAGCCGCCGGCCAg                <  8:278809/149‑1 (MQ=255)                                                                                                                                    gAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCA‑GTTTCAGGGCAAAGCCGCCGGCCAg                >  6:305623/1‑149 (MQ=255)                                                                                                                                       caatgcccgG‑TCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCA‑GTTTCAGGGCAAAGCCGCCGGCCAggg              >  1:297689/5‑148 (MQ=255)                                                                                                                                       caatgcccgG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCA‑GTTTCAGGGCAAAGCCGCCGGCCAggg              >  1:81737/5‑148 (MQ=255)                                                                                                                                         ccgcccgG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCA‑GTTTCAGGGCAAAGCCGCCGGCCAGGGCaa           >  3:134361/1‑149 (MQ=255)                                                                                                                                          tgcccgG‑TCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCA‑GTTTCAGGGCAAAGCCGCCGGCCAGGGCAAt          >  8:295871/2‑149 (MQ=255)                                                                                                                                            cccgG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCA‑GTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGa        <  7:83099/149‑1 (MQ=255)                                                                                                                                             ccgG‑TCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCAGGTTTCAGGGCAAAGCCGCCGGCCAGGGCaa           <  1:33250/146‑1 (MQ=255)                                                                                                                                             ccgG‑TCGTCGATGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCA‑GTTTCAGGGCAAAGCCGCCGGCCAGGGCaa           >  2:33250/1‑145 (MQ=255)                                                                                                                                                    aTCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCA‑GTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGa  <  3:88757/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CGAACTGGCTCACCAGGCGCATCTGGCGAAAACCGGTCTGCTTGAAATCTTCCATGGCTGCGCGCGCAAAGTGCGCTTCGCGGTTGGAGTGCGAAAACAGTTTTACCGTGGTGATGTTGGTGTAAGCGTCAGAAACCCGCCCGG‑TCATCGACGACCGGGCATTGGCCTGCTCCTGCCCGACCTTGCCCAGCCGGGGGACGAAATACAGCATGGCCAAGCCGAACAGCACGACCCAGGCAATGAAGGGCAGCATCA‑GTTTCAGGGCAAAGCCGCCGGCCAGGGCAATGATCGCGA  >  NC_002947/3008167‑3008460

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