Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,916,484	T→A	7.4%	M76K (ATG→AAG)	yebS →	inner membrane PqiA domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,916,484	0	T	A	7.4%	173.2 / 3.7	54	M76K (ATG→AAG)	yebS	inner membrane PqiA domain protein
Reads supporting (aligned to +/- strand):  ref base T (24/26);  new base A (2/2);  total (26/28)
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.88e-01

TTACCAACGTTGCCCGCAATGTGACATGCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTC  >  NC_000913/1916338‑1916626                                                                                                                                                   |                                                                                                                                               ttACCAACGTTGCCCGCAATGTGACATGCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATg                                                                                                                                               >  5:311104/1‑148 (MQ=255)  tACCAACGTTGCCCGCAATGTGACATGCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTATCGCGCCTGGCGGCAATGGCCTTCACCAAGc                                                                                                                                              >  7:216287/1‑148 (MQ=255)  tACCAACGTTGCCCGCAATGTGACATGCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTATCGCGCCTGGCGGCAATGGCCTTCACCAAGc                                                                                                                                              >  7:222357/1‑148 (MQ=255)    ccAACGTTGCCCGCAATGTGACATGCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGTAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTg                                                                                                                                            >  4:242063/1‑148 (MQ=255)          ttGCCCGCAATGTGACATGCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCtt                                                                                                                                      <  6:234508/148‑1 (MQ=255)          ttGCCCGCAATGTGACATGCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCtt                                                                                                                                      <  6:210837/148‑1 (MQ=255)           tGCCCGCAATGTGACATGCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTa                                                                                                                                     >  8:40231/1‑148 (MQ=255)               cGCAATGTGACATGCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGcc                                                                                                                                 >  3:326609/1‑148 (MQ=255)                  aaTGTGACATGCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGtt                                                                                                                              <  2:287483/148‑1 (MQ=255)                       gACATGCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCt                                                                                                                         <  6:311104/148‑1 (MQ=255)                         cATGCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTgg                                                                                                                       >  1:57303/1‑148 (MQ=255)                          aTGCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTggg                                                                                                                      >  2:172453/1‑148 (MQ=255)                            gCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTg                                                                                                                            <  6:301113/140‑1 (MQ=255)                            gCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTg                                                                                                                            >  5:301113/1‑140 (MQ=255)                            gCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGc                                                                                                                    >  1:233749/1‑148 (MQ=255)                                tttAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAAc                                                                                                                >  4:70522/1‑148 (MQ=255)                                    gCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGc                                                                                                             <  1:172453/147‑1 (MQ=255)                                          ccGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCa                                                                                                      >  6:362660/1‑148 (MQ=255)                                                 aaaTTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTgg                                                                                               >  3:88889/1‑148 (MQ=255)                                                  aaTTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTATCGCGCCTGGCGGCAATGGCCTTCACCAAGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGc                                                                                              <  8:216287/148‑1 (MQ=255)                                                  aaTTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTATCGCGCCTGGCGGCAATGGCCTTCACCAAGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGc                                                                                              <  8:222357/148‑1 (MQ=255)                                                    ttCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTg                                                                                            <  2:57303/148‑1 (MQ=255)                                                      ctcATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGt                                                                                           <  2:233749/147‑1 (MQ=255)                                                      ctcATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGGCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGt                                                                                           <  5:362660/147‑1 (MQ=255)                                                         atcaAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGtt                                                                                          <  4:88889/145‑1 (MQ=255)                                                           caAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGtt                                                                                          <  6:343011/143‑1 (MQ=255)                                                           caAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGtt                                                                                          >  5:343011/1‑143 (MQ=255)                                                               gTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCg                                                                                 <  7:40231/148‑1 (MQ=255)                                                                    tATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGc                                                                                                       <  6:97524/121‑1 (MQ=255)                                                                    tATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGc                                                                                                       >  5:97524/1‑121 (MQ=255)                                                                      ttGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGAc                                                                          >  1:288299/1‑148 (MQ=255)                                                                        gTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCa                                                                                                      <  5:346358/118‑1 (MQ=255)                                                                        gTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCa                                                                                                      >  6:346358/1‑118 (MQ=255)                                                                         tCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCCGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGc                                                                        >  3:309259/1‑147 (MQ=255)                                                                           cgcgCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTa                                                                               >  4:87360/1‑138 (MQ=255)                                                                           cgcgCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTa                                                                               <  3:87360/138‑1 (MQ=255)                                                                               cTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTGGGCATTCGTATTGACGCCAACGtt                                                                 <  3:70522/148‑1 (MQ=255)                                                                                      gCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCtt                                                                                                                                      <  5:163892/72‑1 (MQ=255)                                                                                      gCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCtt                                                                                                                                      >  6:163892/1‑72 (MQ=255)                                                                                        aaaaaTTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGt                                                                                >  2:374464/1‑124 (MQ=255)                                                                                        aaaaaTTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGt                                                                                <  1:374464/124‑1 (MQ=255)                                                                                           aaTTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATc                                                     <  3:242063/148‑1 (MQ=255)                                                                                               cGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGa                                                                           <  8:98761/122‑1 (MQ=255)                                                                                               cGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGa                                                                           >  7:98761/1‑122 (MQ=255)                                                                                                 tGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCCGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCaa                                               <  4:309259/148‑1 (MQ=255)                                                                                                    cGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAAt                                             >  8:231347/1‑147 (MQ=255)                                                                                                     gggCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGa                                           <  4:326609/148‑1 (MQ=255)                                                                                                                 tCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGgcg                               <  7:231347/148‑1 (MQ=255)                                                                                                                         gcgcCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTTGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCAGCTGGCAAATGACCAAACAGGGCGATGCGATa                       <  2:288299/148‑1 (MQ=255)                                                                                                                                             tCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCAtt                                                                                   <  6:162373/68‑1 (MQ=255)                                                                                                                                             tCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCAtt                                                                                   >  5:162373/1‑68 (MQ=255)                                                                                                                                             tCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGgcg                               >  2:372806/1‑120 (MQ=255)                                                                                                                                             tCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGgcg                               <  1:372806/120‑1 (MQ=255)                                                                                                                                              cACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTc  >  5:241674/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                               TTACCAACGTTGCCCGCAATGTGACATGCTGTTTAGCCTGCCCGAGATAAATTCTCATCAAAGTGCCTATTGTCCGCGCTGTCAGGCAAAAATTCGTGACGGGCGCGACTGGTCGCTAACGCGCCTGGCGGCAATGGCCTTCACCATGCTGTTGCTTATGCCGTTTGCCTGGGGCGAACCGCTGTTGCATATCTGGCTGTTAGGCATTCGTATTGACGCCAACGTTATGCAAGGCATCTGGCAAATGACCAAACAGGGCGATGCGATAACGGGGTCGATGGTCTTTTTC  >  NC_000913/1916338‑1916626

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