Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,582,309	T→C	Q72Q (CAA→CAG)	ydeN ←	putative sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,582,309	0	T	C	96.6%	199.2 / NA	60	Q72Q (CAA→CAG)	ydeN	putative sulfatase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  major base C (27/30);  minor base A (2/0);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGTTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCTTTCAGCTTTACATCATCTGCCGCATGA  >  NC_000913/1582174‑1582444                                                                                                                                        |                                                                                                                                        aTCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTc                                                                                                                                     >  2:378325/1‑140 (MQ=255)   ccATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCa                                                                                                                                   <  1:244422/140‑1 (MQ=255)      ttAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATaa                                                                                                                                >  1:164760/1‑140 (MQ=255)       tAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAAc                                                                                                                               >  2:977682/1‑140 (MQ=255)         aTGAAAGGAGCGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCa                                                                                                                             >  2:151639/1‑140 (MQ=255)             aaGGAGCGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGat                                                                                                                         <  2:975127/140‑1 (MQ=255)             aaGGAGCGTCGGCGTAGATTTTCGTCCAGCTTCAATGTCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGATTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGat                                                                                                                         <  1:492064/140‑1 (MQ=255)              aGGAGCGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGatc                                                                                                                        >  1:720933/1‑140 (MQ=255)               ggAGCGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGTTGTCCATAACCAAGatca                                                                                                                       >  2:69762/1‑140 (MQ=255)                  gCGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATcc                                                                                                                    <  1:670612/140‑1 (MQ=255)                  gCGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATcc                                                                                                                    <  2:731384/140‑1 (MQ=255)                  gCGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATcc                                                                                                                    <  2:838769/140‑1 (MQ=255)                   cGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATTCCTGATCGAAAGGAAGATGTCGATAACCAAGATCATCCa                                                                                                                   >  2:297514/1‑140 (MQ=255)                   cGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATACCTTATCCAAAGGAAGATGTCCATAACCAAGATCATCCa                                                                                                                   >  2:370569/1‑140 (MQ=255)                     tCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATg                                                                                                                 <  2:921416/140‑1 (MQ=255)                     tCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATg                                                                                                                 <  1:703798/140‑1 (MQ=255)                          gTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAg                                                                                                            <  2:546291/140‑1 (MQ=255)                          gTTGATTTTTGTGCAGCTTCAATGGCTTTACCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAg                                                                                                            <  2:885034/140‑1 (MQ=255)                               tttttGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGa                                                                                                       <  1:692461/140‑1 (MQ=255)                               tttttGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGa                                                                                                       <  1:796248/140‑1 (MQ=255)                                          cTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGGAGCTGTCCATAAACAAAATCATCCAAGGTCCGTACGATAATATTTGGc                                                                                            >  2:340149/1‑140 (MQ=255)                                          cTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGc                                                                                            <  2:550676/140‑1 (MQ=255)                                            tCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCtt                                                                                          >  2:117116/1‑140 (MQ=255)                                            tCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCtt                                                                                          <  1:253189/140‑1 (MQ=255)                                              aaTGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTc                                                                                        <  1:151639/140‑1 (MQ=255)                                                 ggCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCtt                                                                                     <  2:45927/140‑1 (MQ=255)                                                 ggCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCtt                                                                                     <  1:378325/140‑1 (MQ=255)                                                   cTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGGTTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTg                                                                                   <  1:636619/140‑1 (MQ=255)                                                   cTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTg                                                                                   <  1:578063/140‑1 (MQ=255)                                                    tttATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTgg                                                                                  >  2:166134/1‑140 (MQ=255)                                                     ttATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGt                                                                                 <  2:137439/140‑1 (MQ=255)                                                     ttATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGt                                                                                 >  1:899345/1‑140 (MQ=255)                                                            tcCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATAGTCAGTAGGATAATATATGGCCTTCCTTTGGTACTGTAt                                                                          >  2:377047/1‑140 (MQ=255)                                                                  ttttGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTc                                                                    >  2:196042/1‑140 (MQ=255)                                                                  ttttGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTc                                                                    <  2:686449/140‑1 (MQ=255)                                                                         ggTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTaa                                                             >  1:569301/1‑140 (MQ=255)                                                                            aTCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGt                                                          >  2:815214/1‑140 (MQ=255)                                                                              cGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCt                                                        >  1:891338/1‑140 (MQ=255)                                                                                   aCTTCACGATTTTCCATTGTTTATGGGTCCAAAGATCTCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGaa                                                   <  1:340149/140‑1 (MQ=255)                                                                                     ttCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATACTTGGCTTTCCTTCGGTTTTGTATTCGGTCGAGGTAAAGTTTGAGAAAg                                                 >  2:491143/1‑140 (MQ=255)                                                                                          cagtttCCATTGTTTTTGGGTCGAAAGATTCCTTATCTAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACg                                            <  1:466091/137‑1 (MQ=255)                                                                                            ttttCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACgtt                                          >  1:744690/1‑140 (MQ=255)                                                                                               tCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACgtttgt                                       >  1:940919/1‑140 (MQ=255)                                                                                               tCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACgtttgt                                       >  2:722541/1‑140 (MQ=255)                                                                                                        tttGGGTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACAGTATTCTGGCGGCGTAAAGTCTGAGTAAGCTGCGTTTGTCTTGGtctg                              >  2:466091/1‑137 (MQ=255)                                                                                                           gggTCAAAAGATCCCTTATCAATAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCttt                           <  1:496191/140‑1 (MQ=255)                                                                                                           gggTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCttt                           <  1:35677/140‑1 (MQ=255)                                                                                                           gggTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCttt                           <  1:977682/140‑1 (MQ=255)                                                                                                            ggTCAAAAGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCTTTc                          <  2:744690/140‑1 (MQ=255)                                                                                                                aaaaGATCCCTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCTTTCAGCt                      >  1:242335/1‑140 (MQ=255)                                                                                                                 aaaGATCCCTTATCAGAAGGAACCTGTCCATAACCAACATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCTTTCAGCtt                     <  1:491143/140‑1 (MQ=255)                                                                                                                       cccTTATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGGAAAGTCTGGGAAGGCAACGTTTGTTTTGGTGTCTTTCAGTTTTAcatc               >  2:496191/1‑140 (MQ=255)                                                                                                                          ttATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCTTTc                          <  1:956034/126‑1 (MQ=255)                                                                                                                          ttATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCTTTc                          >  2:956034/1‑126 (MQ=255)                                                                                                                           tATCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAACTCTGAGAAAGCACCGTTTGTTTTGGATGCTTTCAGCTTTACATCATCt           >  2:444486/1‑140 (MQ=255)                                                                                                                            aTCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCTTTCAGCTTTACATCATCTg          >  2:141167/1‑140 (MQ=255)                                                                                                                            aTCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCTTTCAGCTTTACATCATCTg          >  1:58580/1‑140 (MQ=255)                                                                                                                            aTCAAAAGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCTTTCAGCTTTACATCATCTg          >  1:669867/1‑140 (MQ=255)                                                                                                                                 aaGGAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCTTTCAGCTTTACATCATCTGCCGCa     <  2:707239/140‑1 (MQ=255)                                                                                                                                    gAAGCTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCTTTCAGCTTTACATCATCTGCCGCATGa  >  2:216799/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                        ATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAGCTTCAATGGCTTTATCTATCCCTATTTTGTAGGTATCGACAACTTCACGATTTTCCATTGTTTTTGGGTCAAAAGATCCCTTATCAAAAGGAAGTTGTCCATAACCAAGATCATCCATGGTCAGTACGATAATATTTGGCTTTCCTTTGGTACTGTATTCTGTCGGCGTAAAGTCTGAGAAAGCAACGTTTGTTTTGGTTGCTTTCAGCTTTACATCATCTGCCGCATGA  >  NC_000913/1582174‑1582444

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