Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,220,860	T→C	V456A (GTC→GCC)	patA →	putrescine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,220,860	0	T	C	100.0%	147.8 / NA	41	V456A (GTC→GCC)	patA	putrescine aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/19);  total (22/19)

ATAACGCCAAAACG‑ATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACA  >  NC_000913/3220752‑3220995                                                                                                              |                                                                                                                                        aataccgcAAAACGAAGCCGCATTAAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCg                                                                                                           <  1:746369/133‑1 (MQ=255)       ccAAAACG‑ATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATACGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATacgac                                                                                                    >  2:17868/1‑140 (MQ=255)        cAAAACG‑ATCCGCATTGAACCGCCACTGACACTTACCACTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATacgacg                                                                                                   >  2:364536/1‑140 (MQ=255)                    gaattacaCCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGa                                                                                        <  1:467867/133‑1 (MQ=255)                     cATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGAt                                                                                       <  1:1006627/140‑1 (MQ=255)                     cATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGACCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGAt                                                                                       >  2:1006627/1‑140 (MQ=255)                      aTTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATc                                                                                      >  2:324181/1‑140 (MQ=255)                              gCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATACGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTtata                                                                              <  1:17868/140‑1 (MQ=255)                                     acactgacCATTGAACAGTGTGAACTGGCGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGAcgcg                                                                       >  2:736033/1‑140 (MQ=255)                                           acCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGcc                                                                 >  2:780101/1‑140 (MQ=255)                                           acCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGcc                                                                 >  1:1070864/1‑140 (MQ=255)                                            ccATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGccc                                                                >  1:139783/1‑140 (MQ=255)                                             cATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCg                                                               <  1:902245/140‑1 (MQ=255)                                               ttGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCa                                                             >  2:873404/1‑140 (MQ=255)                                                tGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAg                                                            >  2:757079/1‑140 (MQ=255)                                                  aaCAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGc                                                          <  1:780101/140‑1 (MQ=255)                                                       tgtgAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGcatcat                                                     >  2:43387/1‑140 (MQ=255)                                                         tgAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATcc                                                   >  1:1060675/1‑140 (MQ=255)                                                            aCTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCt                                                <  2:1060675/140‑1 (MQ=255)                                                                 tGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCt                                                >  2:643154/1‑135 (MQ=255)                                                                 tGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCt                                                <  1:643154/135‑1 (MQ=255)                                                                 tGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGc                                           >  2:1050037/1‑140 (MQ=255)                                                                 tGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGc                                           >  1:329234/1‑140 (MQ=255)                                                                 tGATCAAAGCGGCGCGTAAGGCGCTCGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGc                                           >  2:1085683/1‑140 (MQ=255)                                                                   aTCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGt                                         >  2:948994/1‑140 (MQ=255)                                                                   aTCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGt                                         >  1:863920/1‑140 (MQ=255)                                                                   aTCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGt                                         <  1:873404/140‑1 (MQ=255)                                                                        aGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACa                                    >  2:1148145/1‑140 (MQ=255)                                                                          cggcgCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCt                                                >  2:367881/1‑126 (MQ=255)                                                                          cggcgCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCt                                                <  1:367881/126‑1 (MQ=255)                                                                           ggcgCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCa                                 <  2:139783/140‑1 (MQ=255)                                                                             cgcgTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGa                                                                                        >  2:596267/1‑83 (MQ=255)                                                                             cgcgTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGa                                                                                        <  1:596267/83‑1 (MQ=255)                                                                             cgcgTAAGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACg                               >  1:1081082/1‑140 (MQ=255)                                                                                   aGGCGCTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCa                         <  1:757079/140‑1 (MQ=255)                                                                                      cgcTGGCGGCCATGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTg                      <  2:236616/140‑1 (MQ=255)                                                                                                  tGCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGtt          >  2:587713/1‑140 (MQ=255)                                                                                                   gCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGCTACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCCTCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGttt         <  2:1081082/140‑1 (MQ=255)                                                                                                   gCGAGTAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGttt         <  2:863920/140‑1 (MQ=255)                                                                                                        tAAGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATATGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGca    <  1:932270/140‑1 (MQ=255)                                                                                                          aGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGcaca  <  2:1146291/140‑1 (MQ=255)                                                                                                          aGTGCCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGcaca  <  1:948994/140‑1 (MQ=255)                                                                                                              |                                                                                                                                        ATAACGCCAAAACG‑ATCCGCATTGAACCGCCACTGACACTGACCATTGAACAGTGTGAACTGGTGATCAAAGCGGCGCGTAAGGCGCTGGCGGCCATGCGAGTAAGTGTCGAAGAAGCGTAATATCAAATCGGATGGCGATACGACGTCGCCATCCGATCTTTTTATACGACGCGTACGCCCGCAGGCATCATCCGCTCCGGCGTTAACAGCACGCTTTCACTGCCGTCATCGGTTTCCGCACA  >  NC_000913/3220752‑3220995

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