Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,220,805	T→C	C438R (TGT→CGT)	patA →	putrescine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,220,805	0	T	C	100.0%	114.6 / NA	32	C438R (TGT→CGT)	patA	putrescine aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/13);  total (19/13)

GATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGC  >  NC_000913/3220669‑3220939                                                                                                                                         |                                                                                                                                       gtttgtttttGTTAAAAAAAAAATGGGTTATAATTTTGCCGGGGAGATGTCCCCCCGGCGCGTCCGGGGGGCGGAAACGCAAAAAACCCCAAAAACGACCCGAATTAAACCGCCACTGACACTGACCATTGAACAGcgtg                                                                                                                                     <  2:923495/133‑1 (MQ=255)                      aaTCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGgcgc                                                                                                                >  2:247296/1‑140 (MQ=255)                               taACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGgcgc                                                                                                       >  1:920850/1‑140 (MQ=255)                                aaCTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCt                                                                                                      >  2:87254/1‑140 (MQ=255)                                     cgcgggggagtTGTCCCCCCGGGGCGTGTTGGTGGCCGGACCGCTAAAAACCAAAAAACCGACCCGCATGGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTggcgg                                                                                                 <  2:232815/133‑1 (MQ=255)                                      gCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTggcggc                                                                                                >  1:125963/1‑140 (MQ=255)                                       ccAGGGAGATGTCCCCCCAGGGCTTATCGGGGGCGGGAACGCCAAATAACGCAAAAACAACCCCAATGGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGcc                                                                                               <  2:687174/140‑1 (MQ=255)                                                     cGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAgtgt                                                                                 >  1:189823/1‑140 (MQ=255)                                                           cgcgTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCgaaga                                                                           >  2:594777/1‑140 (MQ=255)                                                             cgTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGaagaag                                                                         >  2:597848/1‑140 (MQ=255)                                                                    gtggCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAata                                                                  >  2:384852/1‑140 (MQ=255)                                                                              cGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGa                                                        <  1:87254/140‑1 (MQ=255)                                                                                          ccAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGa                                                        <  1:821884/128‑1 (MQ=255)                                                                                          ccAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGa                                                        >  2:821884/1‑128 (MQ=255)                                                                                             aaaCGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTc                                         >  2:526749/1‑140 (MQ=255)                                                                                             aaaCGATCCGCATTGAACCGCCACTGACACTGACCATCGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTc                                         >  1:347667/1‑140 (MQ=255)                                                                                                   cccGAATTAAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATTGCGATACGACGTCGCCATc                                   <  2:115862/139‑1 (MQ=255)                                                                                                     cGCATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCg                                 <  1:1056300/140‑1 (MQ=255)                                                                                                      gCATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGa                                <  1:247296/140‑1 (MQ=255)                                                                                                      gCATTAAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGa                                <  2:982244/140‑1 (MQ=255)                                                                                                      ccATTGAACCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATAAAATCGGATGGCGATACGACGTCGCCATCCGa                                <  1:386560/139‑1 (MQ=255)                                                                                                             aCCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTc                                         >  2:973378/1‑124 (MQ=255)                                                                                                             aCCGCCACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTc                                         <  1:973378/124‑1 (MQ=255)                                                                                                                 ccACTGACACTGACCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATac                     <  1:594777/140‑1 (MQ=255)                                                                                                                         tctgacCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCCTTTTATACGACGCGTa             >  1:414711/2‑140 (MQ=255)                                                                                                                            gacCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGc          >  2:526703/1‑140 (MQ=255)                                                                                                                            gacCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGc          >  2:638623/1‑140 (MQ=255)                                                                                                                             acCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGcc         >  1:90942/1‑140 (MQ=255)                                                                                                                             acCATTGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGcc         >  2:962082/1‑140 (MQ=255)                                                                                                                              ccATTGAACAGCGTTAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGccc        >  1:968497/1‑140 (MQ=255)                                                                                                                                  tGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAg    >  2:992302/1‑140 (MQ=255)                                                                                                                                  tGAACAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCCTTTTATACGACGCGTACGCCCGCAg    <  2:414711/140‑1 (MQ=255)                                                                                                                                    aaCAGCGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGc  <  2:367968/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                       GATTGAGTTTGTTGATAACGAAATCGGCTATAACTTTGCCAGCGAGATGTTCCGCCAGCGCGTACTGGTGGCCGGAACGCTCAATAACGCCAAAACGATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGC  >  NC_000913/3220669‑3220939

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