Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,782,883	G→A	intergenic (‑157/+182)	ypjA ← / ← pinH	adhesin‑like autotransporter/pseudogene, invertase resolvase family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,782,883	0	G	A	98.7%	263.9 / ‑4.2	76	intergenic (‑157/+182)	ypjA/pinH	adhesin‑like autotransporter/pseudogene, invertase resolvase family
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (43/32);  total (43/33)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

CAACATCAGCCCACAGGCACCAGACCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATGATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGT  >  NC_000913/2782786‑2782977                                                                                                  |                                                                                               cAACATCAGCCCACAGGCACCAGACCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAaataat                                                                                              <  1:985223/100‑1 (MQ=255)  aaCATCAGCCCACAGGCACCAGACCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAataatg                                                                                             >  2:489932/1‑100 (MQ=255)  aaCATCAGCCCACAGGCACCAGACCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAataatg                                                                                             <  1:300075/100‑1 (MQ=255)     aTCAGCCCACAGGCACCAGACCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGata                                                                                          >  1:885696/1‑100 (MQ=255)         gCCCACAGGCACCAGACCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTcc                                                                                      <  1:582974/100‑1 (MQ=255)           ccACAGGCACCAGACCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGc                                                                                    >  2:1078739/1‑100 (MQ=255)           ccACAGGCACCAGACCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGc                                                                                    <  2:1121258/100‑1 (MQ=255)           ccACAGGCACCAGACCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGc                                                                                    >  2:1369439/1‑100 (MQ=255)            cacaGGCACCAGACCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCa                                                                                   <  1:1298042/100‑1 (MQ=255)             acaGGCACCAGACCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCAt                                                                                  >  2:1240777/1‑100 (MQ=255)               aGGCACCAGACCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATaa                                                                                >  1:1557000/1‑100 (MQ=255)                ggCACCAGACCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAg                                                                               >  2:1390629/1‑100 (MQ=255)                 gCACCAGTCCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGa                                                                              <  2:122234/100‑1 (MQ=255)                    ccagaccagGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATg                                                                           <  2:1256437/100‑1 (MQ=255)                      agaccagGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAg                                                                         <  1:990539/100‑1 (MQ=255)                        accagGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGc                                                                       <  2:100344/100‑1 (MQ=255)                              ggATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCAttttt                                                                 <  2:774238/100‑1 (MQ=255)                               gATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCCTGATGAAATAATGATATTCCGCATAAGAATGAGGCAtttttt                                                                >  2:1106154/1‑100 (MQ=255)                                aTTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTa                                                               <  1:947664/100‑1 (MQ=255)                                 tccATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTaa                                                              >  2:406072/3‑100 (MQ=255)                                  tCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTaan                                                             <  1:531918/100‑2 (MQ=255)                                     tCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGc                                                          >  1:437952/1‑100 (MQ=255)                                      ccTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCa                                                         <  1:355098/100‑1 (MQ=255)                                       cTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAg                                                        >  1:1510783/1‑100 (MQ=255)                                       cTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAg                                                        <  2:156620/100‑1 (MQ=255)                                        tGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGt                                                       >  1:1413366/1‑100 (MQ=255)                                         gAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTg                                                      >  1:1162674/1‑100 (MQ=255)                                            gagaCAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTgcgc                                                   <  2:152293/100‑1 (MQ=255)                                             agaCAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCt                                                  <  1:1403372/100‑1 (MQ=255)                                              gaCAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTg                                                 >  2:286398/1‑100 (MQ=255)                                              gaCAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTg                                                 >  1:478847/1‑100 (MQ=255)                                                   cgcAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGCATGAGGCATTTTTTAAACGCAGTGCGCTGAAgtg                                            >  2:780762/1‑100 (MQ=255)                                                   cgcAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAgtg                                            >  2:1123859/1‑100 (MQ=255)                                                   cgcAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAgtg                                            >  1:276638/1‑100 (MQ=255)                                                    gcAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAgtgt                                           >  2:595106/1‑100 (MQ=255)                                                      aaGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTgg                                         <  2:1426488/100‑1 (MQ=255)                                                        gtgtATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGtt                                       <  2:610751/100‑1 (MQ=255)                                                          gtATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATGATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTgg                                     <  1:88181/100‑1 (MQ=255)                                                          gtATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTgg                                     >  1:1532361/1‑100 (MQ=255)                                                          gtATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTgg                                     <  1:1240785/100‑1 (MQ=255)                                                              tgtgTTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATaa                                 >  1:1252123/1‑100 (MQ=255)                                                               gtgtTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATaaa                                >  2:922077/1‑100 (MQ=255)                                                               gtgtTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATaaa                                >  2:151497/1‑100 (MQ=255)                                                                 gtTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATaaaaa                              >  2:107488/1‑100 (MQ=255)                                                                 gtTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATaaaaa                              >  1:325445/1‑100 (MQ=255)                                                                 gtTCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATaaaaa                              <  1:1362550/100‑1 (MQ=255)                                                                  ttCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGGTaaaaaa                             >  2:1189504/1‑100 (MQ=255)                                                                  ttCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATaaaaaa                             >  1:791117/1‑100 (MQ=255)                                                                   tCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAg                            >  1:472276/1‑100 (MQ=255)                                                                   tCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAg                            >  1:100099/1‑100 (MQ=255)                                                                   tCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAg                            <  1:595109/100‑1 (MQ=255)                                                                   tCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAg                            <  1:971291/100‑1 (MQ=255)                                                                   tCACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTAGTTGGATAAAAAAg                            <  2:181788/100‑1 (MQ=255)                                                                    cACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGt                           >  2:694617/1‑100 (MQ=255)                                                                    cACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGt                           <  1:774721/100‑1 (MQ=255)                                                                    cACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGt                           >  1:762004/1‑100 (MQ=255)                                                                    cACCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGAGTGGTTGGATAAAAAAg                            >  2:827404/1‑99 (MQ=255)                                                                     aCCGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTc                          <  2:885689/100‑1 (MQ=255)                                                                      ccGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCa                         <  1:193480/100‑1 (MQ=255)                                                                      ccGCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCa                         <  1:1189511/100‑1 (MQ=255)                                                                        gCTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAAt                       >  2:873717/1‑100 (MQ=255)                                                                         cTCATCAAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATc                      >  2:1528843/1‑100 (MQ=255)                                                                              caAAGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTc                 <  2:1510773/100‑1 (MQ=255)                                                                               aaaGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTATTTTTTTTTTTCTGAAGTGTGGTTGGATAAAAAAGTTATTCCCTTc                 >  2:1386692/1‑99 (MQ=9)                                                                               aaaGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCa                >  1:1165743/1‑100 (MQ=255)                                                                               aaaGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCa                >  1:463032/1‑100 (MQ=255)                                                                               aaaGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCa                >  1:585242/1‑100 (MQ=255)                                                                                aaGACATCATGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAg               >  2:1164744/1‑100 (MQ=255)                                                                                  atcttcatGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGa             <  2:262172/96‑1 (MQ=255)                                                                                    catcatGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGaaa           <  1:855791/100‑1 (MQ=255)                                                                                      tcatGATGAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGaagat         >  1:262173/1‑97 (MQ=255)                                                                                        atgatgAAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACg       <  1:725410/100‑1 (MQ=255)                                                                                         tgatgaAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACgg      >  2:151446/1‑100 (MQ=255)                                                                                          gatgaAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACggg     >  2:1436806/1‑100 (MQ=255)                                                                                          gatgaAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACggg     <  1:342726/100‑1 (MQ=255)                                                                                           atgaAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACggg     >  2:1528676/1‑99 (MQ=255)                                                                                           atgaAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGc    >  2:476913/1‑100 (MQ=255)                                                                                            tgaAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCg   <  1:1350106/100‑1 (MQ=255)                                                                                             gaAATAATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGt  <  1:1369448/100‑1 (MQ=255)                                                                                                  |                                                                                               CAACATCAGCCCACAGGCACCAGACCAGGGGATTCATCCTGAAGAGACAGCGCAAGTGTATTGTGTTCACCGCTCATCAAAGACATCATGATGAAATGATGATATTCCGCATAAGAATGAGGCATTTTTTAAACGCAGTGCGCTGAAGTGTGGTTGGATAAAAAAGTCAATCCATTCAGGAAATACGGGCGT  >  NC_000913/2782786‑2782977

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