Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP000730	2,018,612	G→A	R103* (CGA→TGA)	mgt ←	transglycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,018,612	0	G	A	100.0%	87.9 / NA	28	R103* (CGA→TGA)	mgt	transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/15);  total (13/15)

GAACGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAAT  >  CP000730/2018473‑2018754                                                                                                                                            |                                                                                                                                               gAACGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCtt                                                                                                                                        >  1:70963/1‑148 (MQ=255)              cATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATaaa                                                                                                                           <  2:165618/148‑1 (MQ=255)                        atTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGCTTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTa                                                                                                                 >  2:387675/1‑148 (MQ=255)                        atTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTa                                                                                                                 >  2:116463/1‑148 (MQ=255)                              gACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATAc                                                                                                           <  2:28842/148‑1 (MQ=255)                              gACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATAc                                                                                                           <  1:249682/148‑1 (MQ=255)                                cAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATActc                                                                                                         >  2:229536/1‑148 (MQ=255)                                          tgtAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTa                                                                                               <  2:298330/148‑1 (MQ=255)                                          tgtAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTa                                                                                               <  2:367589/148‑1 (MQ=255)                                                          ccTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAac                                                                               >  2:188681/1‑148 (MQ=255)                                                          ccTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAac                                                                               >  2:389325/1‑148 (MQ=255)                                                            ttGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGCACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAacta                                                                             >  1:300561/1‑148 (MQ=255)                                                                 caTCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACtttt                                                                        >  1:189712/1‑148 (MQ=255)                                                                     tctGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTAttt                                                                    <  1:8328/148‑1 (MQ=255)                                                                     tctGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTAttt                                                                    <  1:11719/148‑1 (MQ=255)                                                                     tctGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTAttt                                                                    <  1:266810/148‑1 (MQ=255)                                                                                 cGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCtt                                                        <  2:9127/148‑1 (MQ=255)                                                                                 cGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCtt                                                        >  2:346156/1‑148 (MQ=255)                                                                                                   aGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATcg                                      <  1:229536/148‑1 (MQ=255)                                                                                                   aGTTGTACCTTTCAAATCAAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATcg                                      <  1:316390/148‑1 (MQ=255)                                                                                                     ttGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCATTTATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATcgcg                                    <  1:189540/148‑1 (MQ=255)                                                                                                                 aaTCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATaaaaa                        <  1:235056/148‑1 (MQ=255)                                                                                                                  aTCGAATCCATGATGATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATaaaaaa                       >  2:217701/1‑148 (MQ=255)                                                                                                                           atgatgATTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATa              <  1:81966/148‑1 (MQ=255)                                                                                                                                gatTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAAt         >  1:262041/1‑148 (MQ=255)                                                                                                                                 atTGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATa        <  1:372404/148‑1 (MQ=255)                                                                                                                                  ttGTAGAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATaa       >  1:342315/1‑148 (MQ=255)                                                                                                                                       gAATCATTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAAt  >  2:258506/1‑148 (MQ=255)                                                                                                                                            |                                                                                                                                               GAACGATCATTATCATAAAAATAATTTTTGACAACTTGTTGTGTAATGGTACTACCACCTTGCACATCTCTGTCGCTAATCGTTGAAAATAAAGCTCTAGTTGTACCTTTCAAATCGAATCCATGATGATTGTAGAATCGTTCATCTTCCATTGAAATAAAGGCACCTTTAACATACTCTGGCATGTTATCAGCTGACACAAAACTACTTTTATTTTCAATTTTTCTTAGTTCATCCACATTATCGCGTGTAGATAAAAAATACATGATACCAATAAACAAT  >  CP000730/2018473‑2018754

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