Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,411,637	Δ1 bp	coding (453/651 nt)	yfbT ←	sugar phosphatas

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,411,634	0	C	.	98.1%	235.3 / ‑1.4	52	coding (456/651 nt)	yfbT	sugar phosphatas
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (26/25);  total (26/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 = 9.81e-01

GCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGCCC  >  NC_000913/2411536‑2411729                                                                                                   |                                                                                                gCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGcc                                                                                                >  1:819646/1‑100 (MQ=255)   cGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGccca                                                                                              <  2:742449/100‑2 (MQ=255)     gAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGcccagc                                                                                            <  1:1294897/100‑4 (MQ=255)       gCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCcagcag                                                                                         <  1:379609/100‑1 (MQ=255)         gTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCt                                                                                       >  2:271758/1‑100 (MQ=255)              ccGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTgcgcg                                                                                  >  1:1043109/1‑100 (MQ=255)                gCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGcc                                                                                <  1:1233332/100‑1 (MQ=255)                gCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGcc                                                                                <  2:500861/100‑1 (MQ=255)                 cAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCt                                                                               <  2:813749/100‑1 (MQ=255)                   atgacatgacAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTa                                                                              >  2:204737/1‑99 (MQ=255)                       catgacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGa                                                                         >  1:420704/1‑100 (MQ=255)                          gacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATAc                                                                      >  1:358906/1‑100 (MQ=255)                          gacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATAc                                                                      >  1:359338/1‑100 (MQ=255)                          gacaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATAc                                                                      <  1:1444255/100‑1 (MQ=255)                           acaACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACg                                                                     <  1:1402852/100‑1 (MQ=255)                            caccccGCCGCCAGCCCAGAAATCACGCCACCGGGCGCATCTTCCACCACCAAACACTCCTGCGCCGCAAGCCCCAGCAGCTGCGCGCATAACAAATAc                                                                      <  2:237147/96‑1 (MQ=255)                              aCCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCAt                                                                  >  1:353189/1‑100 (MQ=255)                                 cgccgcCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAg                                                               >  1:1473515/1‑100 (MQ=255)                                  gccgccAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAgg                                                              <  1:781198/100‑1 (MQ=255)                                  gccgccAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCACCAgg                                                              >  2:59939/1‑100 (MQ=255)                                     gccaggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTc                                                           >  1:764350/1‑100 (MQ=255)                                       caggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTg                                                         <  2:1165380/100‑1 (MQ=255)                                       caggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTACGGCGCAAG‑CCCAACAGCTGCACGCCTAACAGCTACGCATCAGGTTccg                                                         >  2:161568/1‑98 (MQ=255)                                        aggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTgg                                                        <  2:770386/100‑1 (MQ=255)                                        aggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTgg                                                        <  1:982417/100‑1 (MQ=255)                                         ggccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGt                                                       >  1:1087093/1‑100 (MQ=255)                                           ccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGttt                                                     <  2:805520/100‑1 (MQ=255)                                           ccagAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGttt                                                     >  2:992451/1‑100 (MQ=255)                                             agAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTc                                                   <  1:1279033/100‑1 (MQ=255)                                                   cccGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGGTCTGGTTTTCCGCGCt                                             >  2:106975/3‑100 (MQ=255)                                                   cACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCt                                             >  2:233820/1‑100 (MQ=255)                                                     cGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTc                                           <  2:1137228/100‑1 (MQ=255)                                                     cGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTc                                           <  2:1137412/100‑1 (MQ=255)                                                     cGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTc                                           <  1:1180487/100‑1 (MQ=255)                                                     cGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTc                                           <  1:387345/100‑1 (MQ=255)                                                       ccAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCAc                                         <  1:1466509/100‑1 (MQ=255)                                                       ccAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCAc                                         >  1:1193524/1‑100 (MQ=255)                                                        cAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACt                                        >  2:1556241/1‑100 (MQ=255)                                                                gCATCTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGc                                >  1:211222/1‑100 (MQ=255)                                                                    cTTCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCATCg                               <  2:1477373/97‑1 (MQ=255)                                                                      tCCACCACCACACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTAc                          <  1:493688/100‑1 (MQ=255)                                                                       ccaccaccaCACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACa                         <  2:1208193/100‑1 (MQ=255)                                                                         accaccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACaaa                       >  1:769033/1‑100 (MQ=255)                                                                             ccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACAcc                   <  1:1556251/100‑1 (MQ=255)                                                                             ccacACACTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACAcc                   >  2:598711/1‑100 (MQ=255)                                                                                cacacTCCTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACctct                >  1:520264/1‑100 (MQ=255)                                                                                      ccTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCg          >  1:786998/1‑100 (MQ=255)                                                                                       cTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCgg         >  1:15142/1‑100 (MQ=255)                                                                                       cTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCgg         >  1:286512/1‑100 (MQ=255)                                                                                       cTGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCgg         >  1:341766/1‑100 (MQ=255)                                                                                        tGCGGCGCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCggg        <  2:761191/100‑1 (MQ=255)                                                                                         gcggcgCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTTGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGa       <  2:755320/100‑1 (MQ=255)                                                                                           ggcgCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAg     <  1:1357447/100‑1 (MQ=255)                                                                                           ggcgCAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAg     >  2:1286895/1‑100 (MQ=255)                                                                                              gcAAG‑CCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGccc  <  1:481190/100‑1 (MQ=255)                                                                                                   |                                                                                                GCCGGAGCGTTAACCGCAATGACATGACAACCCGCCGCCAGGCCAGAAAGCACGCCAGCGGGAGCATCTTCCACCACCACACACTCCTGCGGCGCAAGCCCCAGCAGCTGCGCGCCTAACAGATACGCATCAGGTTCTGGTTTTCCGCGCTTCACTCGCTCAGCGGTTACAAACACCTCTGGTGCGGGAAGCCC  >  NC_000913/2411536‑2411729

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