Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,147,379	T→C	I54M (ATA→ATG)	alkA ←	DNA‑3‑methyladenine glycosylase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,147,379	0	T	C	100.0%	125.5 / NA	38	I54M (ATA→ATG)	alkA	DNA‑3‑methyladenine glycosylase 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/25);  total (13/25)

CCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCA  >  NC_000913/2147245‑2147511                                                                                                                                       |                                                                                                                                     ccTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATc                                                                                                                                 <  2:62239/140‑1 (MQ=255)          cccAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGc                                                                                                                        >  2:493847/1‑140 (MQ=255)          cccAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGc                                                                                                                        >  1:672935/1‑140 (MQ=255)          cccAACGCACCGTTAACAATTTGCGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGGATGGCGGGCCATCTCCCGAATAGc                                                                                                                        >  2:112732/1‑140 (MQ=255)                 cACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATTCGGAATAAGAGTcacc                                                                                                                 >  2:83721/1‑140 (MQ=255)                      ttAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCAcgcc                                                                                                            <  1:493847/140‑1 (MQ=255)                         aCAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCcgcg                                                                                                         >  1:276528/1‑140 (MQ=255)                           aaTTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGat                                                                                                       >  2:290058/1‑140 (MQ=255)                              ttGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATAtt                                                                                                    >  2:320229/1‑140 (MQ=255)                              ttGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATAtt                                                                                                    >  1:20394/1‑140 (MQ=255)                              ttGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATAtt                                                                                                    >  1:634130/1‑140 (MQ=255)                              ttGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATAtt                                                                                                    <  2:634130/140‑1 (MQ=255)                              ttGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAACAGTCACCACGCCGCGATAtt                                                                                                    >  2:669056/1‑140 (MQ=255)                                  gggTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGcc                                                                                                <  2:20394/140‑1 (MQ=255)                                   ggTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGccc                                                                                               >  1:949873/1‑140 (MQ=255)                                    gTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCa                                                                                              <  1:320229/140‑1 (MQ=255)                                        cATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGc                                                                                          >  2:197311/1‑140 (MQ=255)                                           tGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAg                                                                                       >  2:729690/1‑140 (MQ=255)                                           tGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAg                                                                                       <  2:404213/140‑1 (MQ=255)                                             cAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGAc                                                                                     >  2:824571/1‑140 (MQ=255)                                               gATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTa                                                                                   <  2:939727/140‑1 (MQ=255)                                                aTCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGcc                                                                                                >  1:153141/1‑126 (MQ=255)                                                aTCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGcc                                                                                                <  2:153141/126‑1 (MQ=255)                                                 tCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACg                                                                                 <  1:1015288/140‑1 (MQ=255)                                                 tCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACg                                                                                 <  2:64241/140‑1 (MQ=255)                                                   aaaCAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACggg                                                                               <  1:290058/140‑1 (MQ=255)                                                          cggcTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCataata                                                                        >  2:489415/1‑140 (MQ=255)                                                           ggcTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCAtaataa                                                                       <  2:498342/140‑1 (MQ=255)                                                           ggcTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCAtaataa                                                                       <  1:729690/140‑1 (MQ=255)                                                               cATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGt                                                                   >  1:149539/1‑140 (MQ=255)                                                                      gCCGGACCCTCTGCGGCAACAGGTTCTAAACCGGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGAc                                                            <  1:112732/140‑1 (MQ=255)                                                                      gCCAGACCCGCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCCCCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGAc                                                            <  2:714409/140‑1 (MQ=255)                                                                      gCCAGACACTCTGCGGCAACAGGTTCTAAACCTTCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGAc                                                            <  1:83721/140‑1 (MQ=255)                                                                      gCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGAc                                                            <  2:949873/140‑1 (MQ=255)                                                                      gCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGAc                                                            <  1:198040/140‑1 (MQ=255)                                                                      gCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGAc                                                            <  1:56217/140‑1 (MQ=255)                                                                       ccAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGAcc                                                           >  2:550241/1‑140 (MQ=255)                                                                              ctctGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCca                                                    <  2:672935/140‑1 (MQ=255)                                                                                      cAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTc                                            >  2:664582/1‑140 (MQ=255)                                                                                            gTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCa                                      <  1:550241/140‑1 (MQ=255)                                                                                                       gCACTTAAATTTATGTGCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCgag                           <  2:149539/140‑1 (MQ=255)                                                                                                                    tgtgCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTACTCACCGCACGGGCGGCGAGAAATCCCAACATc              <  2:276528/140‑1 (MQ=255)                                                                                                                     gtgCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATcc             >  2:75814/1‑140 (MQ=255)                                                                                                                     gtgCAGAGTATGGCGGGCCATATCCGGAATAGCAGTCACCACCCCGCGGTATTCGCCCCCCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTTCACACCGCTCACCGCCCGGGCCGCGAGAAATCCCAACAAcc             >  2:939445/1‑140 (MQ=255)                                                                                                                          gagTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAc        <  1:489415/140‑1 (MQ=255)                                                                                                                           agTATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGAcc       <  2:432824/140‑1 (MQ=255)                                                                                                                             tATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACg          >  1:648252/1‑135 (MQ=255)                                                                                                                             tATGGCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACg          <  2:648252/135‑1 (MQ=255)                                                                                                                                ggCGGGCCATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCa  >  1:996708/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                     CCTAACCTGCCCAACGCACCGTTAACAATTTGTGGGTTACATTGCAGATCAAACAGGCGGCTCATTTTCGCCAGACACTCTGCGGCAACAGGTTCTAAACCTGCACTTAAATTTATGTGCAGAGTATGGCGGGCTATATCCGGAATAGCAGTCACCACGCCGCGATATTCGCCCACCGCCAGACTACGGGCATAATAACTGTCCGCGACCGTTTCCACACTGCTCACCGCACGGGCGGCGAGAAATCCCAACATCCACGACCAGTCA  >  NC_000913/2147245‑2147511

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