Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,296,381:1	+GC	100%	intergenic (+587/+55)	gltP → / ← yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,380	1	.	C	92.1%	104.5 / 0.6	39	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base C (20/15);  minor base T (3/0);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01
*	NC_000913	4,296,380	2	.	G	100.0%	110.0 / ‑6.7	39	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base G (23/14);  minor base C (0/1);  total (23/16)

TTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACC  >  NC_000913/4296347‑4296470                                     ||                                                                                           ttAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                            <  1:2111813/100‑1 (MQ=255)  tAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAccaacg                           <  2:286270/100‑5 (MQ=11)  tAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                           <  2:2534666/100‑1 (MQ=255)   aGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                          <  1:1884044/100‑1 (MQ=255)   aGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                          <  2:2116478/100‑1 (MQ=255)    gCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTATGTTagacggtg                           >  2:699214/1‑91 (MQ=1)    gCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         >  2:17226/1‑100 (MQ=255)    gCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         >  2:955761/1‑100 (MQ=255)    gCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         <  1:464330/100‑1 (MQ=255)    gCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         <  1:1645233/100‑1 (MQ=255)    gCGTGATTTTGCAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         >  2:1965369/1‑100 (MQ=14)     cGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                        >  2:1878298/1‑100 (MQ=255)      gTGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                       >  2:326832/1‑100 (MQ=255)       tGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         <  2:425299/97‑1 (MQ=255)       tGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                      >  1:2474464/1‑100 (MQ=255)       tGATTTTGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                      >  1:588772/1‑100 (MQ=255)          ttttGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                   >  1:1472016/1‑100 (MQ=255)          ttttGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                   >  1:2074730/1‑100 (MQ=255)           tttGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                  >  2:762170/1‑100 (MQ=255)            ttGTAGGTCGGATAAGGCGTT‑‑TATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTATGTTAGccagtggaata                      >  2:1130856/1‑86 (MQ=1)            ttGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                 <  1:797855/100‑1 (MQ=255)            ttGGAGGTCCGCTAAGTCGTT‑‑TCCCCCGCATCCTCCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                 <  2:791899/100‑1 (MQ=14)             tGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                >  1:804474/1‑100 (MQ=255)             tGTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCGTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                >  2:19769/1‑100 (MQ=17)              gTAGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGTCTATCTTAACCGTTGGTTAATTTTTCa               >  2:455776/1‑100 (MQ=17)                aGGTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGACTTATCAGACCTATCTTAACCGTTGGTTAATTTTTCaaa             >  2:387629/1‑100 (MQ=17)                aGGTAGGATAAGGCGTTTCTG‑‑GCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCTACGGTTaagataggc               <  1:2205754/98‑8 (MQ=11)                 ggTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTCGCCGCGTCTTATCAGGCCTATCTTCACCGTTGGTTAATCTTTCaa              <  2:619701/98‑1 (MQ=14)                  gTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct           >  1:453462/1‑100 (MQ=255)                  gTCGGATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct           >  1:1871226/1‑100 (MQ=255)                    cGGATAAGGCGTT‑‑TATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct         >  1:517065/1‑100 (MQ=14)                     ggATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt        >  2:1460006/1‑100 (MQ=255)                      gATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc       >  2:2157271/1‑100 (MQ=255)                      gATAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc       <  2:1103614/100‑1 (MQ=255)                       aTAAGGCGTT‑‑TACGCCTCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca      <  2:459653/100‑1 (MQ=18)                       aTAAGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca      <  1:1768527/100‑1 (MQ=255)                          aGGCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac   >  1:656499/1‑100 (MQ=255)                           ggCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc  <  2:685707/100‑1 (MQ=255)                           ggCGTT‑‑TACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc  >  2:893370/1‑100 (MQ=255)                                     ||                                                                                           TTAGCGTGATTTTGTAGGTCGGATAAGGCGTT‑‑TA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACC  >  NC_000913/4296347‑4296470

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