Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,082,793	G→A	intergenic (+122/+84)	loiP → / ← speB	metalloprotease LoiP/agmatinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,082,793	0	G	A	100.0%	145.4 / NA	41	intergenic (+122/+84)	loiP/speB	metalloprotease LoiP/agmatinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (21/20);  total (21/20)

ATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATC  >  NC_000913/3082717‑3082870                                                                             |                                                                              aTCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGc                                                                           <  1:535525/81‑1 (MQ=255) aTCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGc                                                                           >  1:535569/1‑81 (MQ=255)   cTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGcc                                                                          >  2:535578/1‑80 (MQ=255)     gACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGaa                                                                       >  2:535563/1‑81 (MQ=255)      aCATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAAc                                                                      <  1:535526/81‑1 (MQ=255)        aTGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGt                                                                    <  2:535527/81‑1 (MQ=255)           aaTGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAgg                                                                 >  2:535567/1‑81 (MQ=255)            aTGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGt                                                                >  2:535581/1‑81 (MQ=255)             tGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTc                                                               <  2:535528/81‑1 (MQ=255)               gCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCg                                                              <  1:535529/80‑1 (MQ=255)                  caTTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATa                                                          <  1:535530/81‑1 (MQ=255)                    tttGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAg                                                        >  1:535584/1‑81 (MQ=255)                       gTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGc                                                      <  1:535531/80‑1 (MQ=255)                        ttACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGt                                                    >  2:535593/1‑81 (MQ=255)                            cTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCAc                                                <  1:535532/81‑1 (MQ=255)                               gtgCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgcc                                             >  2:535585/1‑81 (MQ=255)                               gtgCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgcc                                             >  2:535579/1‑81 (MQ=255)                                tgCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgcc                                             <  2:535533/80‑1 (MQ=255)                                 gcgcATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgccg                                            >  1:535571/1‑80 (MQ=255)                                  cgcATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCa                                          >  1:535575/1‑81 (MQ=255)                                  cgcATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCa                                          >  2:535572/1‑81 (MQ=255)                                  cgcATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTACCGAACGTAGGTCGGATAAGGCGTTCACGCCGCa                                          >  2:535561/1‑81 (MQ=255)                                    cATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATc                                        <  2:535534/81‑1 (MQ=255)                                     aTTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATcc                                       >  2:535587/1‑81 (MQ=255)                                     aTTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATcc                                       >  1:535586/1‑81 (MQ=255)                                         ccGGATGCGATGCTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGAc                                    <  1:535535/80‑1 (MQ=255)                                            gatgcgatgcTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTTGGATAAGGCGTTCACGCCGCATCCGACAtt                                 >  2:535577/1‑80 (MQ=255)                                              tgcgatgcTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAAt                              <  1:535536/81‑1 (MQ=255)                                                cgatgcTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATgg                            <  2:535538/81‑1 (MQ=255)                                                cgatgcTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATg                             <  2:535537/80‑1 (MQ=255)                                                  atgcTGGCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCa                          >  1:535576/1‑81 (MQ=255)                                                       ggCGCATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGttt                     >  1:535597/1‑81 (MQ=255)                                                        gcgcATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtttt                    <  1:535540/81‑1 (MQ=255)                                                        gcgcATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtt                      <  2:535539/79‑1 (MQ=255)                                                         cgcATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCTGCATCCGACATTAATGGCACGTTTTa                   >  1:535582/1‑81 (MQ=255)                                                         cgcATCTTATCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTa                   <  1:535541/81‑1 (MQ=255)                                                                 aTCCGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCg               <  2:535542/77‑1 (MQ=255)                                                                    cGGCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcat        <  1:535543/81‑1 (MQ=255)                                                                     ggCCTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatc       <  2:535544/81‑1 (MQ=255)                                                                        cTACGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcgca    >  2:535594/1‑81 (MQ=255)                                                                          aCGAGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCgcatcgcatc  >  2:535609/1‑81 (MQ=255)                                                                             |                                                                              ATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATC  >  NC_000913/3082717‑3082870

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