Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,340,076	G→C	13.3%	intergenic (+61/‑270)	ribAB‑II → / → PP_3814	3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,340,076	0	G	C	13.3%	134.7 / 8.9	51	intergenic (+61/‑270)	ribAB‑II/PP_3814	3,4‑dihydroxy‑2‑butanone‑4‑phosphate synthase/polyamine ABC transportersubstrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (18/25);  new base C (4/3);  total (22/29)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGT  >  NC_002947/4339929‑4340220                                                                                                                                                     |                                                                                                                                                 gggCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAgg                                                                                                                                                 <  7:38312/149‑1 (MQ=255)   gccgccGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGt                                                                                                                                               <  5:101894/149‑1 (MQ=255)          cTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGAcacca                                                                                                                                        <  2:349266/149‑1 (MQ=255)             aaaTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGCGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACaccacc                                                                                                                                      <  3:316537/148‑1 (MQ=255)               aTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTTCCCGCGCATGGGAGGGTGCCCCCCCCggc                                                                                                                                   >  8:166345/1‑147 (MQ=255)               aTACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGAc                                                                                                                                   <  3:356526/149‑1 (MQ=255)                       ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGTGTGACACCACCGACGCATTcgc                                                                                                                           <  7:233900/149‑1 (MQ=255)                       ggCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc                                                                                                                           <  6:174514/149‑1 (MQ=255)                       ggCTGGCGGCCTATGAGCTGGAGGTGCTGCAGAGCAACCCGTTCGACCCGGGTATGGATGGCTCACAACACACAACCAATTTATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGATGGTGACACCACCGACGCATTcgc                                                                                                                           <  7:232037/149‑1 (MQ=38)                        gCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGAATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTcgc                                                                                                                           >  2:219851/1‑148 (MQ=255)                         cTGGCGGGCTATGAGCTGGAGGTGGTGCATAGCATCCCGTTCGTCCCGGGGTTGGATGGCTGAGAACACACAAGGAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                         <  8:283595/149‑1 (MQ=255)                         cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGAGCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCgg                                                                                                                         <  7:330574/149‑1 (MQ=25)                         cTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCgg                                                                                                                         <  7:98182/149‑1 (MQ=255)                                gCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATCGCTCAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAAcc                                                                                                                  <  5:24870/149‑1 (MQ=255)                                  tATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGCCGCATTCGCCGGTAAACCCg                                                                                                                >  1:45492/1‑149 (MQ=255)                                    tGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCt                                                                                                              <  3:91952/149‑1 (MQ=255)                                          ggAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGt                                                                                                                                               >  3:10923/1‑110 (MQ=255)                                          ggAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGt                                                                                                                                               <  4:10923/110‑1 (MQ=255)                                          ggAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCCCCGGCGCCTTCGCGGGTAAACCCGCTCCTACa                                                                                                        >  5:177471/1‑149 (MQ=255)                                           gAGGTGGGCCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTCAGAACAAACAAGCAATTTATTTCGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCTAATGCGAGGGTTACACCACCGACGCATTCGCGGTTAAACCCGCTCCTACa                                                                                                        <  2:55793/148‑1 (MQ=37)                                                             ccGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTAACCGCGAATGCGAGGGTGAAACCCCCGACGCATTAGCGGGTAAAACCCCAGCTACAGGGCCCTTGGTGCTTTTGa                                                                                     >  5:332120/1‑149 (MQ=38)                                                               gTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTAACCGCGACTGCGCGGGTGACACCACCGACGCATTAGCGGGGAAACCCGCGCCGACAGGGCCCGTGGTGGGTTTGa                                                                                     >  6:73418/1‑147 (MQ=37)                                                                      ccGGGGATGGATGGCTTAGAACACACAAGAAATTGTTTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCTAATGCGAGGGTGACACCACCGACGCATTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGa                                                                            <  2:45492/149‑1 (MQ=38)                                                                             tggatggCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTAACCGCGAATGCGAGGGTGGCACCACCGACGCAATAGCGGGTAAACCCGCTCCTACAGGGCCCGTGGTGGTTTTGACACGGCCGAGTGAttg                                                                     >  5:107971/1‑147 (MQ=37)                                                                               gatggCTGAGAACACACACGCAATTTATTTCGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGcc                                                                   <  5:73418/149‑1 (MQ=255)                                                                                    cTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGGTTTTACCCCCAAAT‑CGAGGGGGGCCCCACCCCCGCCTTAGCGGGGAAACCCCCGCCGACAAGGGCCGTGGGGGTTTTTAGAGGGGGGGGTGGTGGGCggg                                                                >  8:182590/1‑146 (MQ=255)                                                                                    cTGACAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGCGCGGG‑TTTACCCGCTAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGt                                                                <  7:182590/147‑1 (MQ=255)                                                                                      gataaCACACATGCATTTTATTTCGGGGGGCCGTGTTGGCGGGGG‑TTTACCGGCTATTGCTTGGGTTACACCACCGACGCTTTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGc                                                             <  6:332120/145‑1 (MQ=255)                                                                                            acacaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGTGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAAcc                                                                                                                  >  8:319508/1‑89 (MQ=255)                                                                                            acacaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAAcc                                                                                                                  <  7:319508/89‑1 (MQ=255)                                                                                              acaAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTAACCGCGCATGGGAGGGTGGCACCCCCGCCGCATTAGCGGGGAAACCCCCGGCGACAGGGCCCGTGGTGGTTTTGACGGGGCGGAGTGATGGCCGGTGGCCCCCTaac                                                     >  4:111120/1‑146 (MQ=255)                                                                                                aaGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTAACCGCGAATGCGAGGGTGACACCCCCGACGCATTAGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACctct                                                  >  5:173908/1‑149 (MQ=255)                                                                                                            accaagggCCCTGTAGGAGCGGG‑TTTAACCGCGAATGGGTCGGTGGTGTGACCCACGCCTTAGCGGGTGAACCCGCTCCGACAGGGCCCGTGGTGGTTTTGACAGGGCGGGGTGATGGGCGGTCGCCCCCAACCTCTTGCCAAAATAtt                                      >  2:262928/6‑149 (MQ=9)                                                                                                             aGCGGGGCCCTGTAGGCGGGTG‑TTTACCCGCTAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCtt                                      <  7:283107/148‑1 (MQ=255)                                                                                                              ggggggCCCTGTTGGTGCGGG‑TTTACCCGGTAATGGGAGGGTGACACCACCGACGAATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                     <  6:107971/146‑1 (MQ=38)                                                                                                               cGGGGCCGTTTAGGAGGGGT‑TTTACCCGCTAATGGGTGGGTTACCCCACCGACCAATTCGGGTGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGa                                   <  1:120607/149‑1 (MQ=255)                                                                                                                    ccgtGTAGGCGCGGG‑TTTACCCGCTAAGGCGTGGGTGACACCACCGACGCATTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATAcc                              <  3:123612/146‑1 (MQ=38)                                                                                                                          aGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTAGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGc                                                                    >  8:61259/1‑105 (MQ=255)                                                                                                                           ggAGCGGG‑TTAACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGTTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACGGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCataata                        <  3:111120/148‑1 (MQ=255)                                                                                                                                         ccGTGAATGCGTGGGTGACACCACCGACGCTTTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa          <  6:177471/149‑1 (MQ=255)                                                                                                                                         ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                     >  6:374933/1‑122 (MQ=21)                                                                                                                                         ccGCGAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTg                                     <  5:374933/122‑1 (MQ=21)                                                                                                                                         ccGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGa          >  5:147229/1‑149 (MQ=255)                                                                                                                                          cgcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCCCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAAAACCATAATATGATAGACGGTAGAc         >  1:128626/1‑149 (MQ=21)                                                                                                                                          cgcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGAAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAc         <  1:262928/149‑1 (MQ=25)                                                                                                                                           gcgAATGCGTCGGTGGTGTCACCCTCGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc        >  5:228545/1‑149 (MQ=25)                                                                                                                                           gcgAATGCGAGGGTGACACCACCGATGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGAcc        >  1:104921/1‑149 (MQ=255)                                                                                                                                           gcgAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCGTTTGCCAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATAAGATATCCCGTAGccc        >  6:11838/1‑149 (MQ=255)                                                                                                                                            cgAATGCGAGGGTGACACCACCGATGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCt       <  2:104921/149‑1 (MQ=255)                                                                                                                                                 gCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATAc                               >  2:95595/1‑120 (MQ=255)                                                                                                                                                 gCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATAc                               <  1:95595/120‑1 (MQ=255)                                                                                                                                                 gCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGt  >  5:5961/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GGGCCGCCGCTGAAATACGCGGGGCTGGCGGGCTATGAGCTGGAGGTGGTGCAGAGCATCCCGTTCGACCCGGGGATGGATGGCTGAGAACACACAAGCAATTGATTTAGCGGGGCCCTGTAGGAGCGGG‑TTTACCCGCGAATGCGAGGGTGACACCACCGACGCATTCGCGGGTAAACCCGCTCCTACAGGGCCCTTGGTGCTTTTGACAGGGCGGAGTGATGGCCGGTCGCCCCATACCTCTTGCCAAAATCTTGGAATACCATAATATGATATCCCGTAGACCTTGAGT  >  NC_002947/4339929‑4340220

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑