Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,314,800	A→C	11.2%	Q438P (CAG→CCG) ‡	rcsD →	RcsD phosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,314,800	0	A	C	11.2%	77.5 / 2.5	43	Q438P (CAG→CCG) ‡	rcsD	RcsD phosphotransferase
Reads supporting (aligned to +/- strand):  ref base A (19/17);  new base C (3/2);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

CATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGA  >  NC_000913/2314658‑2314934                                                                                                                                               |                                                                                                                                       cATTTCAACCTGCAAAACATCCCCACCATGGCGGCACCGCATCAGGTGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCa                                                                                                                                        <  1:302970/143‑1 (MQ=255)  aTTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAg                                                                                                                                       >  1:11716/1‑143 (MQ=255)      gactCTGGAAATCATCACCTCCATGGCGGAACAGCATGAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATc                                                                                                                                   <  2:28906/140‑1 (MQ=255)        aTCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATcgc                                                                                                                                 >  1:329874/1‑143 (MQ=255)         tCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGCCGCTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCGGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGCTAGGGATcgcg                                                                                                                                >  1:480028/1‑143 (MQ=255)           tGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGATATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaa                                                                                                                              <  2:290295/143‑1 (MQ=255)           tGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGCTCAGGAGCGCGaa                                                                                                                              >  2:679081/1‑143 (MQ=255)           tGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGaa                                                                                                                              >  1:620861/1‑143 (MQ=255)           ggCAAAACATCACCACCAGGGGGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGaa                                                                                                                              <  1:731163/142‑1 (MQ=255)            gCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAAGTTTCACTATTCGCGATCAGCACCGCGAAg                                                                                                                             >  2:812850/1‑143 (MQ=255)               aaaCATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCCGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCTGATCGCGaataat                                                                                                                          <  1:420456/143‑5 (MQ=255)                 aCATCACCACCATGGCGGAACAGCATCAGGGGAATATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGGACACAAATTTTCATTATTCGCGATCAGGATCGGGGAGTGCTg                                                                                                                        >  1:235213/1‑143 (MQ=255)                   aTCCCCACCATGCCGGAACCGCATCAGGGGGCTATTCAGGCGACGATCAATAACGAGCTGTCTGAAATCCGCATGTTCCGCAGCCAGGCCGCGCCGAGCACACAAACTTGCATTATTCGCGATCAGGATCGCGCAGTGCTggg                                                                                                                      >  2:647214/1‑142 (MQ=255)                     caccacCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTaa                                                                                                                    >  1:356006/1‑143 (MQ=255)                       ccaccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAAc                                                                                                                  >  2:442286/1‑143 (MQ=255)                       ccaccaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCAGCAGCCAGGTCGCGCCGCGCACACAAATTGTCATTATTCGCGCACGGGCTCGCGAAgtggtgacaaag                                                                                                                  >  2:237928/1‑137 (MQ=255)                         accaTGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTCTTCGCGATCAGGATCGCGAAGTGCTGGTAAACaa                                                                                                                >  2:245443/1‑143 (MQ=255)                             tGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGaaa                                                                                                            >  1:529205/1‑143 (MQ=255)                             tGGAGGAACAGCATCAGGGGATTATTCATGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGaaa                                                                                                            <  2:63668/143‑1 (MQ=255)                                 ggAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCa                                                                                                        <  1:395313/143‑1 (MQ=255)                                         aTCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGACACTCAAGCAGgcg                                                                                                >  1:598250/1‑143 (MQ=255)                                          tCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCCGGATCGCGAAGTGCTGGTAAACAAGAAACTCTAGCAGgcgc                                                                                               >  1:381761/1‑143 (MQ=255)                                                       cAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTTTga                                                                                  <  1:696977/143‑1 (MQ=255)                                                       cAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATga                                                                                  <  1:442286/143‑1 (MQ=255)                                                             acgaTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGACGCGCCGCGCACACAAATTTTCATTCTTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGGATGACaaaaa                                                                            >  2:570636/1‑143 (MQ=255)                                                               gaTCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAAcc                                                                          <  1:748751/143‑1 (MQ=255)                                                                             tGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGgatg                                                            <  2:784168/143‑1 (MQ=255)                                                                                aTGATATCCGCATGTTCCGCAGCCAGGCCGCGCCGCGCACACAAATTTGCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAATCTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATc                                                         <  2:284086/143‑1 (MQ=255)                                                                                aTGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATc                                                         <  2:268132/143‑1 (MQ=255)                                                                                aTGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATc                                                         <  1:621478/143‑1 (MQ=255)                                                                                 tGAGATCCGCATGTTCCGCAGCCAGGTAGCGCGGCGCACACAAATTTTCATTATTCGCGCTCAGGATCGCGAAGTGCAGGTAAACAAGAAACTCAAGCAGGCGCAGCGCCTGCATGAGAAAAACAAGCCGCGGCGGATGATcc                                                        >  2:149121/1‑142 (MQ=255)                                                                                    gctcggCATGTTCCGCAGCCAGGTCGCGGCGCGCACACAAATTTGCATTATTCGCGATCAGGAGCGCGAAGTGCTGGTAAACAAGAAGCTCAGGCAGGGGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTa                                                     <  2:185075/138‑1 (MQ=255)                                                                                    gctCCGCATTGTGCGCAGCCGGGTCGCGCCGCGCACACAAACTATCCTGCTTGGCGTGCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTa                                                     <  2:757745/141‑1 (MQ=255)                                                                                                         ggTCGCGCCGCGCACACAAATTGTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCCGATGATCTTTAGGAAAAACATTGGCGATgcgc                                <  2:430021/143‑1 (MQ=255)                                                                                                          gTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTGAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGgcgctgcgct                               >  1:573993/1‑143 (MQ=255)                                                                                                                     cacaaaAATGTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCg                    <  1:625448/143‑1 (MQ=255)                                                                                                                          aaaTTTTCATTATTCGCGATCAGGATGGCGAAGTGCTGGTAAACAAGAGACGCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGCTGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAg               <  2:527137/143‑1 (MQ=255)                                                                                                                                tCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAGAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTg         >  1:302951/1‑143 (MQ=255)                                                                                                                                  attattCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGc       >  2:124494/1‑143 (MQ=255)                                                                                                                                   cacttCGCGATCCTGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCg      >  1:479431/4‑143 (MQ=255)                                                                                                                                    tattCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCACAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTAGCGATGCGCTGAAAGAAACCGCACAGTCCCTGCCgg     >  2:295563/1‑143 (MQ=255)                                                                                                                                     atgcgcgATCAGGATCGCGAAGGGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGga    <  1:516972/140‑1 (MQ=255)                                                                                                                                      ttCGCGATCCTGATCGCGAAGTGCTGGTAAACGAGTAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAGACATTGGCGATGCGCTGAAAGAACCCGCACAGTCACAGGCGgag   >  2:626435/1‑143 (MQ=255)                                                                                                                                       tCGCGATTAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCCGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGACCCCGCACAGTCCCTGGCGgaga  >  1:326650/1‑143 (MQ=255)                                                                                                                                               |                                                                                                                                       CATTTGAATCTGCAAAACATCACCACCATGGCGGAACAGCATCAGGGGATTATTCAGGCGACGATCAATAACGAGCTGTATGAGATCCGCATGTTCCGCAGCCAGGTCGCGCCGCGCACACAAATTTTCATTATTCGCGATCAGGATCGCGAAGTGCTGGTAAACAAGAAACTCAAGCAGGCGCAGCGTCTGTATGAGAAAAACCAGCAGGGGCGGATGATCTTTATGAAAAACATTGGCGATGCGCTGAAAGAACCCGCACAGTCCCTGGCGGAGA  >  NC_000913/2314658‑2314934

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