Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	475,971	C→G	10.0%	intergenic (+20/+11)	ybaY → / ← ybaZ	outer membrane lipoprotein/excision repair protein, alkyltransferase‑like protein ATL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	475,971	0	C	G	10.0%	124.2 / 5.8	49	intergenic (+20/+11)	ybaY/ybaZ	outer membrane lipoprotein/excision repair protein, alkyltransferase‑like protein ATL
Reads supporting (aligned to +/- strand):  ref base C (22/22);  new base G (3/2);  total (25/24)
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.99e-01

ATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACC  >  NC_000913/475825‑476117                                                                                                                                                   |                                                                                                                                                   aTTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCGACGACTACCGTACCTTCGACTGCACCATCTCCGGTGCATCCGTCTTCGGCCGTGCCCGCTCCGCCGCTATATTAAGAGCCGTTAACCCTCACCGccc                                                                                                                                                  >  1:583142/1‑149 (MQ=255) aTTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGCCTTCACCAACGCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGCCGTTAACCCTCTCCGccc                                                                                                                                                  >  1:993559/1‑149 (MQ=255)  ttGGTGCCGGTAGAGCAAAGCGCCGTTCCGGTTGAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTGCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCg                                                                                                                                                 <  2:147604/149‑1 (MQ=255)        ccGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAgg                                                                                                                                           <  1:163815/149‑1 (MQ=255)         cGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTCAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGGGCGGAGAggg                                                                                                                                          <  1:299277/149‑1 (MQ=255)             aCAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtag                                                                                                                                      <  1:897904/149‑1 (MQ=255)             aCAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtag                                                                                                                                      <  1:173249/149‑1 (MQ=255)              cAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGGCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTtagt                                                                                                                                     <  2:344517/149‑1 (MQ=255)                   aaCCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTc                                                                                                                                <  1:261516/149‑1 (MQ=255)                        ccgTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCGCCTACGCAATATTTAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCg                                                                                                                           >  2:776164/1‑149 (MQ=255)                          gTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGAt                                                                                                                         <  1:484198/149‑1 (MQ=255)                               ggTTCAGGCCAGCGGTGGCGCACCGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCCGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCAGCCCGGAGAGGGTTTGTAGTTCCAGCGATAACGc                                                                                                                    >  1:871434/1‑149 (MQ=255)                               ggTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCCCCAACTCAGGTGAATCCGTCTTCGGCAGGTCCAGCTCCTCCGCAATATTAAGCGCCGTTAACCCTCTCCGCCCGGAGCGGGTTAGTAGTTCCAGCGATAACGc                                                                                                                    >  1:192725/1‑149 (MQ=255)                               ggTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGc                                                                                                                    >  2:932101/1‑149 (MQ=255)                                     ggCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGTTATCGCTGGaac                                                                                                              <  1:381177/149‑2 (MQ=255)                                          gCGGTGGCGCAACGACTAGAGTACCTTCGACTTCCCCAAATCAGGTGAATCCGTCTCCGGCAGTTCCAGCGCCTACGCAATATTACGCGACGTTACCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCGGCAAGTCGAt                                                                                                         <  1:736670/149‑1 (MQ=255)                                           cGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCTAtt                                                                                                        <  1:432257/149‑1 (MQ=255)                                                         cTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGa                                                                                          >  1:108611/1‑149 (MQ=255)                                                            ccGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATTACGCTGCAAGTCGATTTGCCCGCTTCCCgctcc                                                                                       >  2:905911/1‑145 (MQ=255)                                                                   ttCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAATCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATcaca                                                                                >  1:832360/1‑149 (MQ=255)                                                                       aCTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCtt                                                                            <  2:457104/149‑1 (MQ=255)                                                                         ttCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCCGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGTGAGGGGTTGTAGTTACAGCGATATCGCTGCAAGTAGATTTGCCAGCGTCCCGATACCATCACACCTTcg                                                                          >  2:668764/1‑148 (MQ=255)                                                                          tCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTg                                                                         <  1:932101/149‑1 (MQ=255)                                                                           cACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGc                                                                        <  2:130314/149‑1 (MQ=255)                                                                            aCCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTGAGTAGTTCCCGCGATAACGCTGCAAGTCGATTTTCCCGCTTCCCGATACCATCACACctgctgac                                                                       >  1:74657/1‑147 (MQ=255)                                                                             ccAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCa                                                                      >  2:5658/1‑149 (MQ=255)                                                                              cAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAg                                                                     >  2:154427/1‑149 (MQ=255)                                                                                    aGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGc                                                               <  1:511582/149‑1 (MQ=255)                                                                                          aTCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCg                                                         <  1:623754/149‑1 (MQ=255)                                                                                              gTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGa                                                     >  1:69925/1‑149 (MQ=255)                                                                                                    ggCAGTTCCAGCTGCTACGCAATATTAAGCTACGTTAACCCTCTCCTCCCGGAGTGTGTTAGTAGTTCCAGCGATAACGCTGAAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGAAGCTGt                                               <  1:71323/149‑1 (MQ=255)                                                                                                    ggCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGt                                               <  1:430763/149‑1 (MQ=255)                                                                                                       aGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTaaa                                            <  1:400191/149‑1 (MQ=255)                                                                                                          tCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtcc                                         <  1:136237/149‑1 (MQ=255)                                                                                                            cAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccgg                                       <  2:22486/149‑1 (MQ=255)                                                                                                               cTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggtcc                                    >  2:845672/1‑149 (MQ=255)                                                                                                               cTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggtcc                                    >  1:37276/1‑149 (MQ=255)                                                                                                                 ccTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCTTCCCACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggtccgg                                  >  2:808365/1‑149 (MQ=255)                                                                                                                 ccTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggtccgg                                  >  1:925600/1‑149 (MQ=255)                                                                                                                 ccTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAAtccggtccgg                                  >  2:913501/1‑149 (MQ=255)                                                                                                                                gCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGc                   >  2:412211/1‑149 (MQ=255)                                                                                                                                      ctAACCCTCTCCGGGCGGAGGGGGTTAGTAGTTCCAGCGATACCGCTGCAAGTCGATTTGCCCGCTTCCCGATACCCTCACACCTTCTGCCAGTAATGCCTGTCGCTGCCGCTGTAACTCCGGTCCGGTTAGCGAAATTGTGCCGTGGc             >  1:133418/2‑149 (MQ=255)                                                                                                                                       tAACCCTCTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCCCACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAATTCCGGTCCGGTTATCGAAATTGTGCCGTGGCg            >  2:1004107/1‑149 (MQ=255)                                                                                                                                        aaCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCCATAACGCTGCATGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCAGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGa           <  2:558327/149‑1 (MQ=255)                                                                                                                                         aCCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCAGCCAGTACTGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGGGCCGTGGCGAt          <  2:333906/149‑1 (MQ=255)                                                                                                                                         aCCCACTCCGGGCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCGGCTACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGGTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGAt          >  1:187499/1‑149 (MQ=255)                                                                                                                                          cccTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTACCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGCCGCTGTATCTCCGGTCCGGTTAGCGAAAATGTGCCGTGGCGAtt         >  1:346750/1‑149 (MQ=255)                                                                                                                                            ctctCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCAGATCCCATCACACCTTCTGCCAGTAATCCCTGTCTCTGACGCTGAAAAACCCGTCCGGTTAGCGATCTTGTCCAGAGGCGCTTaa       >  1:573001/1‑149 (MQ=255)                                                                                                                                                 cGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCGCACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAaccacc  <  2:925600/149‑1 (MQ=255)                                                                                                                                                 cGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAaccacc  >  1:488577/1‑149 (MQ=255)                                                                                                                                                 cGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGACGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAaccacc  >  2:98695/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   ATTGGTGCCGGTACAGCAAACCGCCGTGCCGGTTCAGGCCAGCGGTGGCGCAACGACTACCGTACCTTCGACTTCACCAACTCAGGTGAATCCGTCTTCGGCAGTTCCAGCTCCTACGCAATATTAAGCGACGTTAACCCTCTCCGCCCGGAGAGGGTTAGTAGTTCCAGCGATAACGCTGCAAGTCGATTTGCCCGCTTCCCGATACCATCACACCTTCTGCCAGTAATGCCTGTCGCTGACGCTGTAAATCCGGTCCGGTTAGCGAAATTGTGCCGTGGCGATTAACCACC  >  NC_000913/475825‑476117

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