Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,239,199	T→C	N31D (AAC→GAC)	mglA ←	methyl‑galactoside ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,239,199	0	T	C	100.0%	204.6 / NA	57	N31D (AAC→GAC)	mglA	methyl‑galactoside ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/32);  total (25/32)

ATACCAAACAGGCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTTATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCGGAGGACGGAGTCGTTGAGCTGACCATATAA  >  NC_000913/2239105‑2239293                                                                                               |                                                                                               aTACCAAACAGGCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCaa                                                                                           >  1:116240/1‑100 (MQ=255)  tACCAAACAGGCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAg                                                                                          >  1:1247009/1‑100 (MQ=255)    ccAAACAGGCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTg                                                                                        <  1:490784/100‑1 (MQ=255)    ccAAACAGGCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTg                                                                                        <  2:1004006/100‑1 (MQ=255)     cAAACAGGCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGc                                                                                       >  1:384545/1‑100 (MQ=255)     cAAACAGGCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGc                                                                                       >  1:164357/1‑100 (MQ=255)      aaaCAGGCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGcc                                                                                      <  1:925850/100‑1 (MQ=255)          aGGCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTaa                                                                                  <  2:587351/100‑1 (MQ=255)           ggCATTTTAATAATGTCGATTTTCCTGTACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAac                                                                                 <  2:384541/100‑1 (MQ=255)           ggCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAac                                                                                 <  2:164355/100‑1 (MQ=255)            gCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAaca                                                                                <  2:142041/100‑1 (MQ=255)            gCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAac                                                                                 <  2:1203056/99‑1 (MQ=255)                tttAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAg                                                                            <  1:811953/100‑1 (MQ=255)                 ttAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAgg                                                                           <  1:347181/100‑1 (MQ=255)                 ttAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAgg                                                                           <  2:529859/100‑1 (MQ=255)                  taataaTGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGa                                                                          >  2:619264/1‑100 (MQ=255)                   aataatGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGaa                                                                         <  2:70815/100‑1 (MQ=255)                      aatGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAgg                                                                      >  2:951817/1‑100 (MQ=255)                         gTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACt                                                                   >  1:1084972/1‑100 (MQ=255)                           cGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACttg                                                                 >  1:237300/1‑100 (MQ=255)                           cGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACttg                                                                 <  2:851066/100‑1 (MQ=255)                              ttttCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACttgttg                                                              <  1:726724/100‑1 (MQ=255)                                ttCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGAt                                                            <  1:647704/100‑1 (MQ=255)                                      cACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGc                                                       <  2:465914/99‑1 (MQ=255)                                       aCCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTc                                                     <  1:440793/100‑1 (MQ=255)                                        ccGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTGTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCt                                                      <  2:239142/98‑1 (MQ=255)                                          gTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCAtt                                                  <  2:1221911/100‑1 (MQ=255)                                            tttCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTc                                                <  1:1071151/100‑1 (MQ=255)                                            tttCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTc                                                <  2:1132551/100‑1 (MQ=255)                                             ttCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTcc                                               >  2:254396/1‑100 (MQ=255)                                             ttCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTGGTTGATACCGCTCATTTcc                                               >  2:106750/1‑100 (MQ=255)                                              tCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCca                                              <  1:503427/100‑1 (MQ=255)                                               cccccATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCcaa                                             <  2:1246998/100‑1 (MQ=255)                                                ccccCTTAATGCATGGATAGAATGTGGCCGGCCTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCcaac                                            <  1:711274/100‑1 (MQ=255)                                                 cccATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCcaaca                                           >  1:1240309/1‑100 (MQ=255)                                                 cccATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCcaaca                                           >  1:570493/1‑100 (MQ=255)                                                   cATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCcaacaan                                         <  2:11992/100‑2 (MQ=255)                                                   cATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAg                                         <  1:1043342/100‑1 (MQ=255)                                                         tGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTcc                                   >  1:1256105/1‑100 (MQ=255)                                                         tGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTcc                                   >  2:728413/1‑100 (MQ=255)                                                         tGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTcc                                   >  2:730215/1‑100 (MQ=255)                                                          gCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTccc                                  <  2:495002/100‑1 (MQ=255)                                                           cATGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTcccc                                 >  1:1226950/1‑100 (MQ=255)                                                            aTGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCg                                <  2:809743/100‑1 (MQ=255)                                                             tGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCACCAAGTATTCCCCgg                               >  1:743697/1‑100 (MQ=255)                                                             tGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCgg                               >  2:898102/1‑100 (MQ=255)                                                             tGGATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCgg                               >  1:812086/1‑100 (MQ=255)                                                              ggATAGAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCgga                              >  2:147664/1‑100 (MQ=255)                                                                   gAATGTGGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCGGAGGACg                         <  2:1134075/100‑1 (MQ=255)                                                                       gtgGCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCGGAGGACGGAGt                     <  2:679035/100‑1 (MQ=255)                                                                         ggCCGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCGGAGGACGGAGTCg                   >  2:1303148/1‑100 (MQ=255)                                                                            cGGACTTTTAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCGGAGGACGGAGTCGTTg                <  1:1038441/100‑1 (MQ=255)                                                                                  tttAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCGGAGGACGGAGTCGTTGAGCTGa          >  2:980616/1‑100 (MQ=255)                                                                                  tttAAATTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCGGAGGACGGAGTCGTTGAGCTGa          >  1:414434/1‑100 (MQ=255)                                                                                      aaTTAACGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCGGAGGACGGAGTCGTTGAGCTGACCat      <  1:951825/100‑1 (MQ=255)                                                                                          aaCGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCGGAGGACGGAGTCGTTGAGCTGACCAtata   >  2:419495/1‑99 (MQ=255)                                                                                          aaCGTCATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCGGAGGACGGAGTCGTTGAGCTGACCATATaa  >  2:419489/1‑100 (MQ=255)                                                                                               |                                                                                               ATACCAAACAGGCATTTTAATAATGTCGATTTTCCTGCACCGTTTTCCCCCATTAATGCATGGATAGAATGTGGCCGGACTTTTAAATTAACGTTATCAAGTGCCTTAACACCAGGAAAGGACTTGTTGATACCGCTCATTTCCAACAAGTATTCCCCGGAGGACGGAGTCGTTGAGCTGACCATATAA  >  NC_000913/2239105‑2239293

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