Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,031,372	T→A	13.3%	intergenic (+45/+30)	PP_1806 → / ← kdsA‑II	arabinose 5‑phosphate isomerase/2‑dehydro‑3‑deoxyphosphooctonate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,031,372	0	T	A	13.3%	113.3 / 9.1	45	intergenic (+45/+30)	PP_1806/kdsA‑II	arabinose 5‑phosphate isomerase/2‑dehydro‑3‑deoxyphosphooctonate aldolase
Reads supporting (aligned to +/- strand):  ref base T (19/20);  new base A (4/2);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-01

AAGCTTACATGCTGGAGAACAAAATCCGTGCCTTGGCGGTGGTCGACGATTGCGGGGCTATTGTTGGCGTCGTCGAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAATGGCAGTGCGCTCGGGCCATCACACTTGGCCTGA  >  NC_002947/2031224‑2031519                                                                                                                                                     |                                                                                                                                                    aaGCTTACATGCTGGAGAACAAAATCCGTGCCTTGGCGGTGGTCGACGATTGCGGGGCTATTGTTGGCGTCGTCGAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGt                                                                                                                                                     >  5:44933/1‑149 (MQ=255)           gCTGGAGAACAAAATCCGTGCCTTGGCGGTGGTCGACGATTGCGGGGCTATTGTTGGCGTCGTCGAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCgg                                                                                                                                           >  5:325881/1‑149 (MQ=255)                    cAAAATCCGTGCCTTGGCGGTGGTCGACGATTGCGGGGCTATTGTTGGCGTCGTCGAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCa                                                                                                                                  >  5:248496/1‑149 (MQ=255)                            gTGCCTTGGCGGTGGTCGACGATTGCGGGGCTATTGTTGGCGTCGTCGAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGc                                                                                                                                          >  1:158956/1‑133 (MQ=255)                            gTGCCTTGGCGGTGGTCGACGATTGCGGGGCTATTGTTGGCGTCGTCGAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGc                                                                                                                                          <  2:158956/133‑1 (MQ=255)                                   ggCGGTGGTCGACGATTGCGGGGCTATTGTTGGCGTCGTCGAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTa                                                                                                                     <  8:127916/147‑1 (MQ=255)                                    gCGGTGGTCGACGATTGCGGGGCTATTGTTGGCGTCGTCGAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAgg                                                                                                                   <  5:90900/148‑1 (MQ=255)                                       gtggtCGACGATTGCGGGGCTATTGTTGGCGTCGTCGAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCaa                                                                                                                >  1:142067/1‑148 (MQ=255)                                              acgaTTGCGGGGCTATTGTTGGCGTCGTCGAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAAt                                                                                                        >  3:154304/1‑149 (MQ=255)                                                gaTTGCGGGGCTATTGTTGGCGTCGTCGAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTg                                                                                                      >  6:199693/1‑149 (MQ=255)                                                    gCGGGGCTATTGTTGGCGTCGTCGAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCgg                                                                                                   <  2:12771/148‑1 (MQ=255)                                                          cTATTGTTGGCGTCGTCGAGATTTTTGACGGCAATCCCGGAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGAt                                                                                            <  6:248496/149‑1 (MQ=255)                                                                      tcgtcgAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGatca                                                                                >  3:291297/1‑149 (MQ=255)                                                                        gtcgAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGatcatc                                                                              >  8:283803/1‑149 (MQ=255)                                                                        gtcgAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGatcatc                                                                              >  2:62529/1‑149 (MQ=255)                                                                         tcgAGATATTTGACGGCAATCCAGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGAAATAGTCAATTGCGTGAATGATTTTACAAGatcatc                                                                              >  8:129572/1‑148 (MQ=255)                                                                           gagaTTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCa                                                                            <  4:291297/148‑1 (MQ=255)                                                                               tttttGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCt                                                                       <  1:322692/149‑1 (MQ=255)                                                                                 tttGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCtt                                                                      <  4:154304/148‑1 (MQ=255)                                                                                   tGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTg                                                                    <  2:142067/148‑1 (MQ=255)                                                                                      cGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGAttt                                                                <  6:44933/149‑1 (MQ=255)                                                                                       ggCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTtg                                                               >  5:228278/1‑149 (MQ=255)                                                                                          aaTCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGc                                                            <  4:292199/149‑1 (MQ=255)                                                                                                 tAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGAtttt                                                                                         <  6:227626/113‑1 (MQ=255)                                                                                                 tAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGAtttt                                                                                         >  5:227626/1‑113 (MQ=255)                                                                                                   gCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAgg                                                    <  7:283803/148‑1 (MQ=255)                                                                                                    cGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGt                                                  <  6:325881/149‑1 (MQ=255)                                                                                                                  tAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCaa                                    <  6:228278/149‑1 (MQ=255)                                                                                                                   aaTTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAAt                                   >  8:102682/1‑149 (MQ=255)                                                                                                                           ttGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAATGGCAGTgc                           <  1:325236/149‑1 (MQ=255)                                                                                                                           ttGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAATGGCAGTgc                           >  2:311191/1‑149 (MQ=255)                                                                                                                                      ggCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAATGGCAGTGCGCTCGGGCCAt                <  7:102682/149‑1 (MQ=255)                                                                                                                                        aaGCCGCTTCTGGACAAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGc                                            >  7:241035/2‑119 (MQ=255)                                                                                                                                        aaGCCGCTTCTGGACAAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGc                                            <  8:241035/118‑1 (MQ=255)                                                                                                                                        aaGCCGCTTCTGGACAAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAATGGCAGTGCGCTCGGGCCATca              >  1:10138/2‑149 (MQ=255)                                                                                                                                         aGCCGCTTCTGGACAAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAATGGCAGTGCGCTCGGGCCATca              >  8:232606/1‑148 (MQ=255)                                                                                                                                          gCCGCTTCTGGACAAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAATGGCAGTGCGCTCGGGCCATcaca            <  7:232606/149‑1 (MQ=255)                                                                                                                                           ccGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTtgt                                                              >  4:145866/1‑98 (MQ=255)                                                                                                                                           ccGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTtgt                                                              <  3:145866/98‑1 (MQ=255)                                                                                                                                           ccGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAATGGCAGTGCGCTCGGGCCATCacac           >  6:70249/1‑149 (MQ=255)                                                                                                                                           ccGCTTCTGGACAAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAATGGCAGTGCGCTCGGGCCATCacac           >  7:148509/1‑149 (MQ=255)                                                                                                                                            cGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGAAATAGTCAATTGCGGGAATGATTTTACCAGATCATAGACCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAATGGCAGTGCGCTCGGGCCATCacac           <  7:129572/148‑1 (MQ=255)                                                                                                                                                  ttGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAATGGCAGTGCGCTCGGGCCATCACACTTGGCCt    >  7:298013/1‑149 (MQ=255)                                                                                                                                                    gTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAATGGCAGTGCGCTCGGGCCATCACACTTGGCCt    <  5:70249/147‑1 (MQ=255)                                                                                                                                                    gTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAATGGCAGTGCGCTCGGGCCATCACACTTGGCCTGa  >  3:259137/1‑149 (MQ=255)                                                                                                                                                    gTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTACAATGTCAGTGATCTCGGGCAATCACACTTGGCCTGa  >  2:155201/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AAGCTTACATGCTGGAGAACAAAATCCGTGCCTTGGCGGTGGTCGACGATTGCGGGGCTATTGTTGGCGTCGTCGAGATTTTTGACGGCAATCCCGTAGCGTGAAACAGCCGGTAATTATCGTTGCCTGGGTAGGCAGCCGCTTCTTGTCCAGAAGCGGCTTTACTCACCGTCGTAGATTAGGCAATAGTCAATTGCGGGAATGATTTTACCAGATCATCCAGCGCTTTGATTTGTGCCAGGAAGGGTTCCAGCTTGTCCAATGGCAGTGCGCTCGGGCCATCACACTTGGCCTGA  >  NC_002947/2031224‑2031519

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