Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP009273	3,019,631	G→T	M169I (ATG→ATT)	guaD →	guanine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,019,631	0	G	T	100.0%	108.0 / NA	38	M169I (ATG→ATT)	guaD	guanine deaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (13/25);  total (13/25)

ATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAAC  >  CP009273/3019488‑3019776                                                                                                                                                |                                                                                                                                                  ttGTCGGCGGTCTTCTTCAAGCAGCTTTTACGTAACCGAACCACCACGGTGCTGGGGGTTGGCACTGTTCATCCGCAATCTTTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGa                                                                                                                                              <  2:356357/148‑1 (MQ=255)  gtgcGGCGTTCGTCTTCCAGCAGGTTTGACGTAACGGAACGCCCACGGCGCGGGTGTTTGTCACTTGTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGAc                                                                                                                                             <  2:117144/146‑1 (MQ=255)    tCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCg                                                                                                                                           <  2:60168/149‑1 (MQ=255)        cGTTCTTCATTAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAAc                                                                                                                                       <  1:218117/149‑1 (MQ=255)                atcaAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGa                                                                                                                               >  2:684182/1‑149 (MQ=255)                 tcaAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTAATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTGTGGACCGCAACGCACCGGag                                                                                                                              >  1:103571/1‑148 (MQ=255)                    agcagcTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTAt                                                                                                                           >  2:540881/1‑149 (MQ=255)                          ttttACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGCGATTATGAACCGCAACGCACCGGATTATCTGCTc                                                                                                                     >  2:476034/1‑149 (MQ=255)                                   aCGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGcc                                                                                                            >  1:218865/1‑149 (MQ=255)                                                      ggTGTTGGGCACTGTTCATCCGCAATCTGTGTCTGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGGATGATTGCTGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTACTCGACACTGCCGAAAGCAGCTATCACccaa                                                                                         <  2:348866/149‑1 (MQ=255)                                                       ggggttGGCACTGTTCATCCCCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACcaaag                                                                                        <  2:11261/145‑1 (MQ=255)                                                        tgtTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACcaaagc                                                                                       <  2:218865/149‑1 (MQ=255)                                                           ttGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACcaaagcaaa                                                                                    <  1:585091/149‑1 (MQ=255)                                                               cACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAAc                                                                                <  2:672793/149‑1 (MQ=255)                                                                 cTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATCTCAATATACGTATGATTGCCGGTCAAGCGATTATGGACCGCAACGCACCGGATTATCTGCCCGAGACTACCCAAAACAGCCAGCCCCAAAGCCAACAACTg                                                                              >  2:694549/1‑149 (MQ=255)                                                                   ggtCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGAt                                                                            <  1:68429/147‑1 (MQ=255)                                                                   gTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGAt                                                                            <  2:560032/149‑1 (MQ=255)                                                                   gTTCATACGCGATCTGTTGATGCGATGTTTGAAGACGCCAGTCATATCAAAATGCGTATGATAGCCGGTAAGGTGCTTATGGACCGCAACGCAACGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGAt                                                                            <  2:586897/149‑1 (MQ=255)                                                                        cccGCAATCTGTTGATGCGCATTTTGAAACAGCCAGTCATATCACTATGCAAATGATTGCGGGTAAGGTGGTTATGGACCGCAACGCACCGGATTATGTTCTCGGCACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAAc                                                                       <  2:614105/148‑1 (MQ=255)                                                                          cGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGc                                                                     >  1:519082/1‑149 (MQ=255)                                                                            cAATCTGTTGATGCTCTGTTTGAAGCCGCCAGGCATATCACTATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCGCGACACTGCGGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTg                                                                   <  2:362862/149‑1 (MQ=255)                                                                                     gATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCCATCACCAAAGCAAAGAACTGCACGAACGCTGGAACaaaaa                                                          >  2:652777/1‑149 (MQ=255)                                                                                      aTGCGCTGTTTGAAGCCGCCAGCCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAAt                                                         <  1:283912/149‑1 (MQ=255)                                                                                        gcgcTGTTTGAAGCCGCCAGTCATATCAATATGCGGATGATTGCCGGGAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACCCTGCCGAAAGCAGCTATCACCAAAGCAAAGAACATATCGAACGCTGCCACAAAAAAgg                                                       >  2:610849/1‑149 (MQ=255)                                                                                          gcTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGgtc                                                     <  2:731740/149‑1 (MQ=255)                                                                                          gcTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGgtc                                                     <  2:130525/149‑1 (MQ=255)                                                                                             gTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTACCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATGGAACGCTGGCACAAAAATGgtcggc                                                  >  1:469474/1‑147 (MQ=255)                                                                                                 gAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCt                                              <  1:382834/149‑1 (MQ=255)                                                                                                 gAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCt                                              <  2:714243/149‑1 (MQ=255)                                                                                                      cgccAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATg                                         <  2:293456/149‑1 (MQ=255)                                                                                                       gccAGTCTTATCAATATGCGTATGATTGCCGGTAAGGTGAGTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGc                                        <  2:620715/149‑1 (MQ=255)                                                                                                              atatCAATATGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCAAAGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAGAACAGGTCGTCTGCTATATGCGATTACg                                 >  2:701408/1‑149 (MQ=255)                                                                                                                 tCAATATGCGTATGACTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCa                              <  1:68484/149‑1 (MQ=255)                                                                                                                 tCAATACGCGTATGATTGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCCGCTATCACCAAAGCAAAGACCTGAGCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCa                              <  1:305362/149‑1 (MQ=255)                                                                                                                      atGCGTATGATTGCCGCTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCtt                         <  1:576865/149‑1 (MQ=255)                                                                                                                                ttGCCGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGAcc               <  2:559207/149‑1 (MQ=255)                                                                                                                                   ccGGTAAGGTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCa            >  1:140376/1‑149 (MQ=255)                                                                                                                                           gTGATTATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGCCCTCATCTCCtta    >  2:197260/1‑147 (MQ=255)                                                                                                                                            tgattaTGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCCCCCACCGGCGGTTCCACCATATGCGAGTACGCCACGCTTCGCCCCGACCTCATCTCCTGaa   >  2:519694/1‑149 (MQ=255)                                                                                                                                             gattatGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAAc  <  2:467817/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  ATGTCGGCGTTCTTCATCAAGCAGCTTTTACGTAACGGAACCACCACGGCGCTGGTGTTTGGCACTGTTCATCCGCAATCTGTTGATGCGCTGTTTGAAGCCGCCAGTCATATCAATATGCGTATGATTGCCGGTAAGGTGATGATGGACCGCAACGCACCGGATTATCTGCTCGACACTGCCGAAAGCAGCTATCACCAAAGCAAAGAACTGATCGAACGCTGGCACAAAAATGGTCGTCTGCTATATGCGATTACGCCACGCTTCGCCCCGACCTCATCTCCTGAAC  >  CP009273/3019488‑3019776

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