Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,988,905	A→G	T64A (ACC→GCC)	ftnA →	ferritin iron storage protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,988,905	0	A	G	100.0%	176.9 / NA	50	T64A (ACC→GCC)	ftnA	ferritin iron storage protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/23);  total (27/23)

TACTCTTCACTGCTTTATCAGCAAATGAGCGCCTGGTGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGCTCATGCTGCAAT  >  NC_000913/1988770‑1989041                                                                                                                                        |                                                                                                                                         tACTCTTCACTGCTTTATCAGCAAATGAGCGCCTGGTGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCgg                                                                                                                                      >  2:393177/1‑140 (MQ=255) tACTCTTCACTGCTTTATCAGCAAATGAGCGCCTGGTGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCgg                                                                                                                                      <  1:586267/140‑1 (MQ=255)      ttCACTGCTTTATCAGCAAATGAGCGCCTGGTGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAAtt                                                                                                                                 >  2:320282/1‑140 (MQ=255)           tGCTTTATCAGCAAATGAGCGCCTGGTGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACcg                                                                                                                            >  1:161782/1‑140 (MQ=255)               ttATCAGCAAATGAGCGCCTGGTGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTa                                                                                                                        <  2:761355/140‑1 (MQ=255)                  tCAGCAAATGAGCGCCTGGTGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTa                                                                                                                     >  1:484356/1‑140 (MQ=255)                                ccTGGTGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAAtct                                                                                                       >  1:792873/1‑140 (MQ=255)                                ccTGGTGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAAtct                                                                                                       >  1:358899/1‑140 (MQ=255)                                       cAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTg                                                                                                <  1:746034/140‑1 (MQ=255)                                       cAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTg                                                                                                <  1:320282/140‑1 (MQ=255)                                         gCTATCATACCTTCGAAGGTGTTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCt                                                                                              <  1:740387/140‑1 (MQ=255)                                           tATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGa                                                                                            <  1:695989/140‑1 (MQ=255)                                            aTCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGaa                                                                                           <  2:613262/140‑1 (MQ=255)                                              cATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGaa                                                                                           >  1:134631/1‑138 (MQ=255)                                              cATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGaa                                                                                           <  2:134631/138‑1 (MQ=255)                                               aTACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAatat                                                                                        >  1:738900/1‑140 (MQ=255)                                                   cTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCt                                                                                    <  1:106334/140‑1 (MQ=255)                                                   cTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCt                                                                                    <  2:494441/140‑1 (MQ=255)                                                     tCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCa                                                                                  >  2:456111/1‑140 (MQ=255)                                                      cGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCAc                                                                                 >  2:759615/1‑140 (MQ=255)                                                              ctgcCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGaa                                                                         <  1:393177/140‑1 (MQ=255)                                                              ctgcCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGaa                                                                         >  1:957077/1‑140 (MQ=255)                                                              ctgcCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGaa                                                                         >  1:641751/1‑140 (MQ=255)                                                               tgcCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAat                                                                        <  2:266077/140‑1 (MQ=255)                                                                  cgcgTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAattat                                                                     <  2:136238/140‑1 (MQ=255)                                                                      ttCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTc                                                                   <  1:798400/138‑1 (MQ=255)                                                                      ttCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTc                                                                   >  2:798400/1‑138 (MQ=255)                                                                      ttCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCa                                                                 >  2:867275/1‑140 (MQ=255)                                                                        ccTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAgg                                                               <  2:358899/140‑1 (MQ=255)                                                                           gcgcCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGaaa                                                            >  1:307025/1‑140 (MQ=255)                                                                           gcgcCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGaaa                                                            >  1:307026/1‑140 (MQ=255)                                                                                 tCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGa                                                              <  2:973399/132‑1 (MQ=255)                                                                                 tCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGa                                                              >  1:973399/1‑132 (MQ=255)                                                                                 tCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCtata                                                      >  2:676969/1‑140 (MQ=255)                                                                                 tCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCtata                                                      >  2:702158/1‑140 (MQ=255)                                                                                 tCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCtata                                                      <  1:759615/140‑1 (MQ=255)                                                                                           gAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACa                                            >  2:693124/1‑140 (MQ=255)                                                                                           gAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACa                                            <  2:792873/140‑1 (MQ=255)                                                                                                        atatGCAGCGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAg                               <  2:484356/140‑1 (MQ=255)                                                                                                                cGTCTGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTaa                       >  1:945646/1‑140 (MQ=255)                                                                                                                    tGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGaa                   >  1:299671/1‑140 (MQ=255)                                                                                                                    tGTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGaa                   <  2:641751/140‑1 (MQ=255)                                                                                                                     gTTTGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAAc                  >  1:792793/1‑140 (MQ=255)                                                                                                                        tGATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTgg               >  2:526527/1‑140 (MQ=255)                                                                                                                         gATTACCTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGc              <  2:259455/140‑1 (MQ=255)                                                                                                                              ccTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGCTCAtg         <  1:867275/140‑1 (MQ=255)                                                                                                                              ccTGACTGATGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGCTCAtg         >  1:474338/1‑140 (MQ=255)                                                                                                                                  actgaTGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGCTCAtgctgc     <  1:150656/140‑1 (MQ=255)                                                                                                                                   ctgaTGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGCTCATGCTGCa    >  2:378822/1‑140 (MQ=255)                                                                                                                                     gaTGCCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGCTCATGCTGCAAt  >  1:98241/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                         TACTCTTCACTGCTTTATCAGCAAATGAGCGCCTGGTGCAGCTATCATACCTTCGAAGGTGCTGCCGCGTTCCTGCGCCGTCACGCCCAGGAAGAGATGACGCATATGCAGCGTCTGTTTGATTACCTGACTGATACCGGCAATTTACCGCGTATTAATACCGTTGAATCTCCGTTTGCTGAATATTCCTCACTTGATGAATTATTCCAGGAAACCTATAAACACGAACAATTAATCACCCAGAAAATTAACGAACTGGCTCATGCTGCAAT  >  NC_000913/1988770‑1989041

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