Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,580,128	A→G	F221F (TTT→TTC)	ydeM ←	putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,580,128	0	A	G	96.0%	249.0 / 0.3	75	F221F (TTT→TTC)	ydeM	putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (40/32);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

CCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTC  >  NC_000913/1580032‑1580223                                                                                                 |                                                                                                cccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAAc                                                                                              >  1:574608/1‑100 (MQ=255) cccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAAc                                                                                              >  1:94605/1‑100 (MQ=255) cccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAAc                                                                                              <  1:1188084/100‑1 (MQ=255)   cATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACtt                                                                                            <  1:802391/100‑1 (MQ=255)    aTTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTg                                                                                           <  1:142853/100‑1 (MQ=255)     ttCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGc                                                                                          <  1:946684/100‑1 (MQ=255)       ccccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCa                                                                                        <  2:1662313/100‑1 (MQ=255)        cccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCAt                                                                                       >  1:844149/1‑100 (MQ=255)        cccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCAt                                                                                       >  1:1753485/1‑100 (MQ=255)          cAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTGTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATaa                                                                                     >  2:1432552/1‑100 (MQ=255)            aaaaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGc                                                                                   <  1:691071/100‑1 (MQ=255)            aaaaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGc                                                                                   >  1:1222518/1‑100 (MQ=255)                aCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTg                                                                               >  2:1296022/1‑100 (MQ=255)                  ggCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTggg                                                                             <  2:1782576/100‑1 (MQ=255)                   gCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGa                                                                            <  2:844146/100‑1 (MQ=255)                   gCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGa                                                                            <  2:432905/100‑1 (MQ=255)                    cTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAg                                                                           <  1:1607943/100‑1 (MQ=255)                       aCAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGca                                                                        >  1:597405/1‑100 (MQ=255)                        cAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGcac                                                                       <  1:1296029/100‑1 (MQ=255)                        cAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGcac                                                                       <  2:224571/100‑1 (MQ=255)                        cAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTG‑TTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGcaca                                                                      <  1:1310851/100‑1 (MQ=255)                          aaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAg                                                                     <  2:94604/100‑1 (MQ=255)                          aaGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAg                                                                     >  1:1211781/1‑100 (MQ=255)                             cTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGaaa                                                                  >  1:1256779/1‑100 (MQ=255)                              tttCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGaaaa                                                                 >  1:406165/1‑100 (MQ=255)                               ttCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGaaaaa                                                                <  2:597403/100‑1 (MQ=255)                               ttCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGaaaaa                                                                >  2:135578/1‑100 (MQ=255)                                    aCTGACGGATGAAAATTTCACCCACATTGTTTTTAATCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                                                           <  1:1432561/100‑1 (MQ=255)                                    aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                                                           <  1:1869289/100‑1 (MQ=255)                                    aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                                                           <  1:111884/100‑1 (MQ=255)                                    aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                                                           <  1:1675245/100‑1 (MQ=255)                                     ctaaCGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                                                           <  2:1800209/96‑1 (MQ=255)                                      tGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGa                                                         >  1:510455/1‑100 (MQ=255)                                       gACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGa                                                         <  1:132761/99‑1 (MQ=255)                                       gACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGAt                                                        >  1:511771/1‑100 (MQ=255)                                       gACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGAt                                                        >  1:805254/1‑100 (MQ=255)                                              gAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGaa                                                  >  1:440158/1‑99 (MQ=255)                                                aaaTTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAAtgt                                               <  2:1104780/100‑1 (MQ=255)                                                   tttCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGtt                                            <  2:406164/100‑1 (MQ=255)                                                    ttCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtc                                           >  1:80079/1‑100 (MQ=255)                                                     tCACCCACATTGTTTTTAACCCATTGCATAAAAAGGGTTAACATGAACGTGCAATAAGACGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtct                                          <  2:440157/100‑1 (MQ=255)                                                       aCCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCa                                        >  1:537280/1‑100 (MQ=255)                                                        cccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCAc                                       >  1:1100438/1‑100 (MQ=255)                                                        cccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCAc                                       >  1:393406/1‑100 (MQ=255)                                                         ccACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACt                                      >  1:277806/1‑100 (MQ=255)                                                           acaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTAt                                    <  2:1874468/100‑1 (MQ=255)                                                             attgTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGa                                  >  1:1579089/4‑100 (MQ=255)                                                             attgTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGa                                  >  1:434984/4‑100 (MQ=255)                                                               tgTTTTTAACCCATTGGATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGAcc                                <  2:556210/99‑1 (MQ=255)                                                               cGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGAcc                                <  2:1222512/100‑1 (MQ=255)                                                                gTTTTTAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCa                               >  2:1823763/1‑100 (MQ=255)                                                                  ttttAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACt                             >  2:409530/1‑100 (MQ=255)                                                                   tttAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTg                            <  2:511770/100‑1 (MQ=255)                                                                   tttAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTg                            <  2:1753471/100‑1 (MQ=255)                                                                    ttAACCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGa                           <  2:823800/100‑1 (MQ=255)                                                                      aaCCCATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGaaa                         >  2:1224684/1‑100 (MQ=255)                                                                        cccATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATc                       <  1:1957585/100‑1 (MQ=255)                                                                          cATTGCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCaa                     >  1:1546229/1‑100 (MQ=255)                                                                              gCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATAtt                 >  2:1947016/1‑100 (MQ=255)                                                                              gCATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATAtt                 >  2:818605/1‑100 (MQ=255)                                                                               cATAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTc                >  2:430471/1‑100 (MQ=255)                                                                                aTAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCGATATTCg               >  1:403923/1‑100 (MQ=255)                                                                                aTAAAAATGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCg               <  2:277805/100‑1 (MQ=255)                                                                                   aaaaTGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCg            >  2:951807/1‑100 (MQ=255)                                                                                   aaaaTGGTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCg            >  2:329854/1‑100 (MQ=255)                                                                                      aTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTcc         <  1:2016937/100‑1 (MQ=255)                                                                                         gTTGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCgg      >  2:269512/1‑100 (MQ=255)                                                                                          ttGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGt     >  1:912347/1‑100 (MQ=255)                                                                                          ttGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGt     >  1:912920/1‑100 (MQ=255)                                                                                           tGACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGtt    <  2:662963/100‑1 (MQ=255)                                                                                            gACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGttt   >  2:770126/1‑100 (MQ=255)                                                                                            gACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGttt   >  1:843371/1‑100 (MQ=255)                                                                                            gACATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGttt   >  1:1458365/1‑100 (MQ=255)                                                                                             aCATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTc  >  2:1658732/1‑100 (MQ=255)                                                                                             aCATGAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTc  >  1:1001344/1‑100 (MQ=255)                                                                                                 |                                                                                                CCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTC  >  NC_000913/1580032‑1580223

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: ‑