Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,182,453	C→A	intergenic (‑19/‑97)	ygiD ← / → zupT	4,5‑DOPA dioxygenase extradiol/heavy metal divalent cation transporter ZupT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,182,453	0	C	A	100.0%	192.7 / NA	54	intergenic (‑19/‑97)	ygiD/zupT	4,5‑DOPA dioxygenase extradiol/heavy metal divalent cation transporter ZupT
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (32/22);  total (32/22)

AAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGA  >  NC_000913/3182376‑3182529                                                                              |                                                                             aaaaaaaaTGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCt                                                                           <  2:591021/81‑1 (MQ=255)    aaCAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGt                                                                        >  1:591085/1‑81 (MQ=255)     aCAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGtt                                                                       <  1:591022/81‑1 (MQ=255)       aaTGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTcg                                                                     >  2:591094/1‑81 (MQ=255)       aaTGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATACTCACAATAAGCTTGTTcg                                                                     >  2:591075/1‑81 (MQ=255)         tGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTcggc                                                                   <  1:591023/81‑1 (MQ=255)         tGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTcggc                                                                   >  1:591099/1‑81 (MQ=255)           cTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTcggcgg                                                                 <  2:591024/81‑1 (MQ=255)           cTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTcggcgg                                                                 <  1:591025/81‑1 (MQ=255)           cTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTcggcg                                                                  >  2:591077/1‑80 (MQ=255)              gCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAg                                                               <  1:591026/80‑1 (MQ=255)               cATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGt                                                              >  2:591070/1‑80 (MQ=255)               cATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTa                                                             <  1:591028/81‑1 (MQ=255)               cATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTa                                                             <  2:591027/81‑1 (MQ=255)                aTACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTaa                                                            <  1:591029/81‑1 (MQ=255)                  aCGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAg                                                           >  2:591071/1‑80 (MQ=255)                  aCGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAg                                                           >  2:591086/1‑80 (MQ=255)                         gAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCgg                                                    <  1:591030/80‑1 (MQ=255)                         gAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGa                                                   >  2:591079/1‑81 (MQ=255)                              cATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAAc                                               >  1:591074/1‑80 (MQ=255)                               atgatgATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACaa                                             >  1:591064/1‑81 (MQ=255)                                tgatgaTATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAat                                            <  2:591031/81‑1 (MQ=255)                                  atgatATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAatga                                          >  2:591109/1‑81 (MQ=255)                                    gatATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAatgat                                         >  2:591107/1‑80 (MQ=255)                                     atatCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAatgatg                                        <  1:591032/80‑1 (MQ=255)                                      tatCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAatgatgat                                      >  2:591097/1‑81 (MQ=255)                                          cTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATgatgatgat                                   <  2:591033/80‑1 (MQ=255)                                          cTTAACTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGatc                                  <  1:591034/81‑1 (MQ=255)                                             aaCTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGatcat                                >  2:591095/1‑80 (MQ=255)                                             aaCTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGatcat                                >  2:591106/1‑80 (MQ=255)                                              aCTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGatcatc                               <  2:591036/80‑1 (MQ=255)                                              aCTAAAGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGAtcatca                              <  2:591035/81‑1 (MQ=255)                                                  aaGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGtt                           <  1:591037/80‑1 (MQ=255)                                                   aGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTAt                         >  1:591091/1‑81 (MQ=255)                                                   aGGTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTAt                         <  1:591038/81‑1 (MQ=255)                                                    ggTGTCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTAtt                        >  2:591076/1‑81 (MQ=255)                                                       gtCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTAtttt                      >  1:591101/1‑80 (MQ=255)                                                       gtCATTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTAtttt                      >  1:591092/1‑80 (MQ=255)                                                          aTTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGAc                   >  1:591110/1‑80 (MQ=255)                                                          aTTTTGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGAc                   <  2:591039/80‑1 (MQ=255)                                                            tttGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGAt                >  1:591089/1‑81 (MQ=255)                                                             ttGATATCCTCTCAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATc               >  2:591088/1‑81 (MQ=255)                                                             ttGATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATc               >  2:591105/1‑81 (MQ=255)                                                               gATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCt              <  2:591040/80‑1 (MQ=255)                                                               gATATCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTg             >  2:591087/1‑81 (MQ=255)                                                                atatCCTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGc            <  1:591041/81‑1 (MQ=255)                                                                     cTCACAATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTg         >  1:591081/1‑79 (MQ=255)                                                                       cacaATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAgg     >  1:591119/1‑81 (MQ=255)                                                                       cacaATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAgg     >  2:591121/1‑81 (MQ=255)                                                                         caATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGt    >  2:591113/1‑80 (MQ=255)                                                                         caATAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGt    >  1:591112/1‑80 (MQ=255)                                                                          aaTAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTg   >  1:591090/1‑80 (MQ=255)                                                                          aaTAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGa  <  1:591042/81‑1 (MQ=255)                                                                           aTAAGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTg   >  1:591103/1‑79 (MQ=255)                                                                              |                                                                             AAAAACAATGCTGGCATACGTGTTGAAGACATGATGATATCCTTAACTAAAGGTGTCATTTTGATATCCTCACAATACGCTTGTTCGGCGGAGTAAGAACCCGGATAACAATGATGATGATCATCAGTTATTTTGACGATCTGCCTGAAGGTGA  >  NC_000913/3182376‑3182529

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