Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,033,458	T→G	I9I (ATA→ATC)	yedR ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,033,458	0	T	G	95.0%	279.0 / 4.2	80	I9I (ATA→ATC)	yedR	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base G (35/41);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CCAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGAA  >  NC_000913/2033361‑2033553                                                                                                  |                                                                                                ccAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACtatgat                                                                                               >  2:2026982/1‑100 (MQ=255)  cAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTatgatg                                                                                              <  2:625400/100‑1 (MQ=255)   aGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGAtgt                                                                                             <  1:1475796/100‑1 (MQ=255)      ttGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGAt                                                                                          <  2:1380105/100‑1 (MQ=255)       tGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATa                                                                                         <  2:1600248/100‑1 (MQ=255)        gccgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATaa                                                                                        >  1:799981/1‑100 (MQ=255)        gccgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATaa                                                                                        >  2:1297980/1‑100 (MQ=255)         ccgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATaaa                                                                                       >  2:1370797/1‑100 (MQ=255)          cgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATaaaa                                                                                      >  2:1747471/1‑100 (MQ=255)          cgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATaaaa                                                                                      >  1:1106974/1‑100 (MQ=255)            ccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGt                                                                                    <  2:646910/100‑1 (MQ=255)            ccAGGGCTGCGAGAGAGTCTCGCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAg                                                                                     <  2:1907561/99‑1 (MQ=255)             cAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTc                                                                                   <  2:1939707/100‑1 (MQ=255)               gggCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTcac                                                                                 >  2:462760/1‑100 (MQ=255)               gggCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTcac                                                                                 <  1:1751430/100‑1 (MQ=255)                ggCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTcaca                                                                                >  1:132821/1‑100 (MQ=255)                 gCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACAt                                                                               <  1:124426/100‑1 (MQ=255)                   tGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACAttt                                                                             >  1:1142769/1‑100 (MQ=255)                    gCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACAtttt                                                                            >  2:1933372/1‑100 (MQ=255)                    gCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACAtttt                                                                            >  2:1392593/1‑100 (MQ=255)                     cGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACAttttt                                                                           >  1:618865/1‑100 (MQ=255)                      gagagaGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTc                                                                          <  2:1492047/100‑1 (MQ=255)                      gagagaGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTc                                                                          <  1:66004/100‑1 (MQ=255)                       agagagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTcc                                                                         <  1:1572034/100‑1 (MQ=255)                       agagagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTcc                                                                         <  1:1561874/100‑1 (MQ=255)                       agagagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTcc                                                                         <  2:1788778/100‑1 (MQ=255)                       agagagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTc                                                                          <  2:425278/99‑1 (MQ=255)                        gagagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCca                                                                        >  1:274802/1‑100 (MQ=255)                        gagagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCca                                                                        <  1:308154/100‑1 (MQ=255)                        gagagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCca                                                                        <  2:290586/100‑1 (MQ=255)                         agagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCcac                                                                       <  2:1474428/100‑1 (MQ=255)                          gagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGTTGTGATAAAAGTCACATTTTTCcaca                                                                      >  1:1103120/1‑100 (MQ=255)                            gTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACAtt                                                                    >  1:609615/1‑100 (MQ=255)                            gTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACAtt                                                                    >  1:1921484/1‑100 (MQ=255)                            gTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACAtt                                                                    >  2:440764/1‑100 (MQ=255)                               tctctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGa                                                                 >  2:335119/1‑100 (MQ=255)                                 tctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATg                                                               <  2:1166133/100‑1 (MQ=255)                                  ctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGa                                                              >  1:1390845/1‑100 (MQ=255)                                    tCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGc                                                            <  1:654062/100‑1 (MQ=255)                                     ccATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGcc                                                           <  2:1312593/100‑1 (MQ=255)                                      cATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCa                                                          <  2:1791340/100‑1 (MQ=255)                                           ttGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCaa                                                     <  2:989951/100‑1 (MQ=255)                                            tGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAg                                                    <  2:132820/100‑1 (MQ=255)                                            tGTCGAACCGCATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAACTCTCATTTTTCCACATTGGATGAGCCAGTCAAg                                                    <  2:1344239/100‑1 (MQ=255)                                             gTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGa                                                   >  1:183353/1‑100 (MQ=255)                                              tCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGAc                                                  <  1:182304/100‑1 (MQ=255)                                               cGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTTCACATTGGATGAGCCAGTCAAGACt                                                 >  2:250217/1‑100 (MQ=255)                                               cGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACt                                                 >  2:1317746/1‑100 (MQ=255)                                               cGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACt                                                 <  2:1906456/100‑1 (MQ=255)                                                gAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTc                                                <  2:722883/100‑1 (MQ=255)                                                 aaCCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCa                                               >  1:2038275/1‑100 (MQ=255)                                                  aCCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCaa                                              >  2:1394897/1‑100 (MQ=255)                                                     gTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCa                                           >  2:1425069/1‑100 (MQ=255)                                                     gTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCa                                           <  1:462761/100‑1 (MQ=255)                                                        ttCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAttt                                        >  1:1068628/1‑100 (MQ=255)                                                        ttCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAttt                                        <  2:1106969/100‑1 (MQ=255)                                                         tCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAtttt                                       <  2:1921466/100‑1 (MQ=255)                                                         tCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAtttt                                       >  2:197205/1‑100 (MQ=255)                                                         tCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAtttt                                       <  2:609613/100‑1 (MQ=255)                                                            aTTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGTTGAGCCAGTCAAGACTCAATCATTTTTTg                                    >  1:685662/1‑100 (MQ=255)                                                            aTTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTg                                    <  2:799978/100‑1 (MQ=255)                                                            aTTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTg                                    <  1:263326/100‑1 (MQ=255)                                                             ttttGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGa                                   >  2:243023/1‑100 (MQ=255)                                                                tGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACa                                >  1:722001/1‑100 (MQ=255)                                                                  agaTAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACaaa                              <  2:685660/100‑1 (MQ=255)                                                                  agaTAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACaaa                              >  1:892512/1‑100 (MQ=255)                                                                    aTAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATa                            >  1:1194891/1‑100 (MQ=255)                                                                          cGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATaaa                      >  1:469747/1‑100 (MQ=255)                                                                          cGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATaaa                      >  2:1990546/1‑100 (MQ=255)                                                                            gAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTTTATATATAGATAAAtt                    >  2:1015700/1‑100 (MQ=255)                                                                                 ttCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgcg               <  1:440765/100‑1 (MQ=255)                                                                                      aCTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACg          <  2:726024/100‑1 (MQ=255)                                                                                       cTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACgg         >  1:1574164/1‑100 (MQ=255)                                                                                         tAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCa       >  2:1505335/1‑100 (MQ=255)                                                                                         tAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCa       >  2:13101/1‑100 (MQ=255)                                                                                           gTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGt     <  2:618863/100‑1 (MQ=255)                                                                                            tACTTTGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTg    <  2:164928/100‑1 (MQ=255)                                                                                            tACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTg    <  1:1534988/100‑1 (MQ=255)                                                                                            tACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTg    <  1:1392602/100‑1 (MQ=255)                                                                                             aCTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGa   <  1:852083/100‑1 (MQ=255)                                                                                             aCTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGa   <  2:892509/100‑1 (MQ=255)                                                                                              cTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGaa  >  1:1108323/1‑100 (MQ=255)                                                                                                  |                                                                                                CCAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGAA  >  NC_000913/2033361‑2033553

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