Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,240,292	G→A	A19A (GCC→GCT)	mglB ←	methyl‑galactoside transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,240,292	0	G	A	100.0%	158.1 / NA	44	A19A (GCC→GCT)	mglB	methyl‑galactoside transporter subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/22);  total (22/22)

GGTCATTCTGAGAATCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTGAGATAAATGCTCGGCGGGGCAGTAGAGTTAATGAAGTGTTACTGAACGCC  >  NC_000913/2240147‑2240434                                                                                                                                                  |                                                                                                                                               ggTCATTCTGAGAATCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCAcc                                                                                                                                           <  1:441458/151‑1 (MQ=255)   tCATTCTGAGAATCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCg                                                                                                                                          >  1:386295/1‑150 (MQ=255)    cATTCTGAGAATCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGa                                                                                                                                         >  2:786779/1‑150 (MQ=255)            gaATCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCTCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTAGTTATAGATTGTTACACCAATGCGAGTATCAGTAGAGTGTGCAGCAGCACCGAATAACATGc                                                                                                                                >  1:758832/1‑151 (MQ=255)              aTCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTg                                                                                                                              >  1:375511/1‑151 (MQ=255)                       cagcagCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATcac                                                                                                                      <  2:246025/150‑1 (MQ=255)                            gcTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCAcagc                                                                                                                   >  1:47204/1‑148 (MQ=255)                                aaCATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAgg                                                                                                            >  2:340477/1‑151 (MQ=255)                                  cATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAggg                                                                                                           >  2:380285/1‑150 (MQ=255)                                      tGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAac                                                                                                      >  2:361413/1‑151 (MQ=255)                                         cgcgGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACAcc                                                                                                   <  2:800886/151‑1 (MQ=255)                                            ggCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACcttc                                                                                                <  2:43397/151‑1 (MQ=255)                                               tttCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACcttctt                                                                                              >  1:746018/1‑150 (MQ=255)                                                ttCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTAt                                                                                            >  1:59053/1‑151 (MQ=255)                                                   gCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACC‑TCTTATTCAt                                                                                        >  2:275450/1‑151 (MQ=255)                                                    cATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCAt                                                                                        >  1:428375/1‑151 (MQ=255)                                                    cATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCAt                                                                                        <  1:265464/151‑1 (MQ=255)                                                        ttGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGt                                                                                     >  1:792515/1‑150 (MQ=255)                                                         tGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATg                                                                                       <  2:317819/147‑1 (MQ=255)                                                               aTAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCgg                                                                             >  1:675026/1‑151 (MQ=255)                                                                   ccTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGtttt                                                                         <  2:572090/151‑1 (MQ=255)                                                                      tGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGttttt                                                                        <  2:762961/149‑1 (MQ=255)                                                                         gcACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTAt                                                                   <  1:389221/151‑1 (MQ=255)                                                                                aGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGt                                                            >  2:216422/1‑151 (MQ=255)                                                                                   cATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGt                                                         <  2:518984/151‑1 (MQ=255)                                                                                   cATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGt                                                         <  2:836936/151‑1 (MQ=255)                                                                                    aTAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTg                                                        >  1:184185/1‑151 (MQ=255)                                                                                      aaaGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCt                                                      >  1:43152/1‑151 (MQ=255)                                                                                       aaGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCtt                                                     <  1:853546/151‑1 (MQ=255)                                                                                        aGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTg                                                    >  1:256624/1‑151 (MQ=255)                                                                                           tATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTgag                                                  <  1:669471/150‑1 (MQ=255)                                                                                              cgtcgtACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTGAGATa                                               >  1:789060/1‑150 (MQ=255)                                                                                                        ataGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTGAGATAAATGCTcggcg                                    <  2:391749/151‑1 (MQ=255)                                                                                                        ataGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTGAGATAAATGCTcggc                                     <  1:402558/150‑1 (MQ=255)                                                                                                        ataGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTGAGATAAATGCTcggc                                     <  2:801688/150‑1 (MQ=255)                                                                                                                  acacCAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTGAGATAAATGCTCGGCGGGGCAGTag                           <  1:358975/150‑1 (MQ=255)                                                                                                                     ccAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTGAGATAAATGCTCGGCGGGGCAGTAGAGt                        <  2:157672/150‑1 (MQ=255)                                                                                                                      cAATGCGAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTGAGATAAATGCTCGGCGGGGCAGTAGAGtt                       >  1:515153/1‑150 (MQ=255)                                                                                                                           cGAGTATCAGCAGCGGGTGCAGCAGCTCCTAATGACATGCTGGCCATCACAGCAGAACGGGTTAACGCCTTCTTATTCATGGTATCTTCGGTTTTGCTTATGCATGGTAGTGCTTGAGATAAATGCTCGTCGGGGCAGTAGAGTTaag                    <  1:849729/148‑2 (MQ=255)                                                                                                                            gAGTATCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTGAGATAAATGCTCGGCGGGGCAGTAGAGTTAATGa                  >  2:274207/1‑149 (MQ=255)                                                                                                                                aTCAGCAGCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTGAGATAAATGCTCGGCGGGGCAGTAGGGTTAATGAAGTGtt            >  2:710505/1‑151 (MQ=255)                                                                                                                                  cagcagCGTGTGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTGAGATAAATGCTCGGCGGGGCAGTAGAGTTAATGAAGTGTTAc          <  2:219986/151‑1 (MQ=255)                                                                                                                                         gtgtGCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTGAGATAAATGCTCGGCGGGGCAGTAGAGTTAATGAAGTGTTACTGAACGc   <  2:32390/151‑1 (MQ=255)                                                                                                                                           gtgCAGCAGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTGAGATAAATGCTCGGCGGGGCAGTAGAGTTAATGAAGTGTTACTGAACGcc  <  2:159672/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                               GGTCATTCTGAGAATCATTCATCAGCAGCTGAACATCTGGCGCGGCTTTCGCATCTTGCTCAATAGCCTTGCGCACTACAGACATAAAGTTATCGTCGTACTTATAGATTGTTACACCAATGCGAGTATCAGCAGCGTGTGCAGCGGCACCGAATAACATGCTGGCCATCACAGCAGACAGGGTTAACACCTTCTTATTCATGGTATCTCCGGTTTTTCTTATGCAGGGTAGTGCTTGAGATAAATGCTCGGCGGGGCAGTAGAGTTAATGAAGTGTTACTGAACGCC  >  NC_000913/2240147‑2240434

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