Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,988,933	T→C	V73A (GTT→GCT)	ftnA →	ferritin iron storage protein (cytoplasmic)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,988,933	0	T	C	100.0%	217.2 / NA	59	V73A (GTT→GCT)	ftnA	ferritin iron storage protein (cytoplasmic)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/34);  total (25/34)

AGCGCCTGGTGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGCTCATGCTGCAATGACCAA  >  NC_000913/1988797‑1989047                                                                                                                                         |                                                                                                                   aGCGCCTGGTGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTg                                                                                                                  <  1:1264172/139‑1 (MQ=255)   cgcCTGGTGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGaa                                                                                                                <  1:75789/139‑1 (MQ=255)          tGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGt                                                                                                         <  2:546341/139‑1 (MQ=255)             aGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTg                                                                                                      <  2:528413/139‑1 (MQ=255)               cTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCt                                                                                                    <  1:740054/139‑1 (MQ=255)                   cATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAat                                                                                                >  2:951295/1‑139 (MQ=255)                    aTACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAata                                                                                               >  2:499564/1‑139 (MQ=255)                         ttCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCt                                                                                          <  1:1096183/139‑1 (MQ=255)                         ttCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCt                                                                                          <  2:554201/139‑1 (MQ=255)                           cGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGt                                                                                                         >  2:85032/1‑122 (MQ=255)                           cGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGt                                                                                                         <  1:85032/122‑1 (MQ=255)                            gAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCAc                                                                                       >  1:1089967/1‑139 (MQ=255)                            gAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCAc                                                                                       >  1:1257325/1‑139 (MQ=255)                               ggTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTtg                                                                                    <  1:249517/139‑1 (MQ=255)                               ggTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTtg                                                                                    <  1:806095/139‑1 (MQ=255)                                    tgcCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGaa                                                                               >  2:898907/1‑139 (MQ=255)                                    tgcCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGaa                                                                               <  1:499564/139‑1 (MQ=255)                                              cTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGAAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAgg                                                                     <  2:305734/139‑1 (MQ=255)                                               tGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGa                                                                    <  2:999007/139‑1 (MQ=255)                                               tGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGa                                                                    <  2:1257325/139‑1 (MQ=255)                                               tGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGa                                                                    <  1:951295/139‑1 (MQ=255)                                                      tCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCtat                                                             >  2:218093/1‑139 (MQ=255)                                                                      aTGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACg                                                      >  2:465332/1‑130 (MQ=255)                                                                      aTGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACg                                                      <  1:465332/130‑1 (MQ=255)                                                                      aTGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTaa                                             >  2:1205307/1‑139 (MQ=255)                                                                      aTGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTaa                                             >  2:435249/1‑139 (MQ=255)                                                                          cGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCAc                                         >  1:743247/1‑139 (MQ=255)                                                                                ttcAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGaa                                   <  1:1208273/137‑1 (MQ=255)                                                                                    gCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGa                                                                    <  2:971609/102‑1 (MQ=255)                                                                                    gCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGa                                                                    >  1:971609/1‑102 (MQ=255)                                                                                     cGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGt                                                                                                         >  2:602501/1‑64 (MQ=255)                                                                                     cGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGt                                                                                                         <  1:602501/64‑1 (MQ=255)                                                                                     cGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAgg                                                                     >  1:639188/1‑100 (MQ=255)                                                                                     cGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAgg                                                                     <  2:639188/100‑1 (MQ=255)                                                                                     cGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTa                              <  1:898907/139‑1 (MQ=255)                                                                                     cGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTa                              <  2:743247/139‑1 (MQ=255)                                                                                       tCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTc                                                                                         <  1:1172033/78‑1 (MQ=255)                                                                                       tCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTc                                                                                         >  2:1172033/1‑78 (MQ=255)                                                                                       tCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAAc                            >  2:1057252/1‑139 (MQ=255)                                                                                                    ctgactgaTACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGCTCAtg               <  2:1089967/139‑1 (MQ=255)                                                                                                     tgactgaTACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAattatt                                                                          <  1:616501/79‑1 (MQ=255)                                                                                                     tgactgaTACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAattatt                                                                          >  2:616501/1‑79 (MQ=255)                                                                                                       actgaTACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGGTGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCAccc                                       >  2:767067/1‑112 (MQ=255)                                                                                                       actgaTACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCAccc                                       <  1:767067/112‑1 (MQ=255)                                                                                                        ctgaTACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGCTCAtgctgc           <  1:435249/139‑1 (MQ=255)                                                                                                         tgaTACCGGCAATTTACCACGTAATAATACCGCTGAATCTCCGTTTGCTGAATATTGCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAATAATAAATCACCCAGAAGATTAACGAACTGGCTCATGCTGCa          <  2:220298/139‑1 (MQ=255)                                                                                                           aTACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGaa                                                                   >  2:175203/1‑80 (MQ=255)                                                                                                           aTACCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGaa                                                                   <  1:175203/80‑1 (MQ=255)                                                                                                             aCCGGCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGCTCATGCTGCAATGa      <  1:218093/139‑1 (MQ=255)                                                                                                                 gCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCt                                                                                          <  2:90534/51‑1 (MQ=255)                                                                                                                 gCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCt                                                                                          >  1:90534/1‑51 (MQ=255)                                                                                                                 gCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCAccc                                       >  1:219048/1‑102 (MQ=255)                                                                                                                 gCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCAccc                                       <  2:219048/102‑1 (MQ=255)                                                                                                                 gCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGCTCATGCTGCAATGACCaa  >  2:39736/1‑139 (MQ=255)                                                                                                                 gCAATTTACCGCGTATTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGCTCATGCTGCAATGACCaa  >  2:297632/1‑139 (MQ=255)                                                                                                                               aTTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGaa                                                                   >  1:314715/1‑60 (MQ=255)                                                                                                                               aTTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGaa                                                                   <  2:314715/60‑1 (MQ=255)                                                                                                                               aTTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAAt                                            <  1:35620/83‑1 (MQ=255)                                                                                                                               aTTAATACCGCTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAAt                                            >  2:35620/1‑83 (MQ=255)                                                                                                                                         |                                                                                                                   AGCGCCTGGTGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGCTCATGCTGCAATGACCAA  >  NC_000913/1988797‑1989047

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