Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 640,862 G→A A370A (GCG→GCA ahpF → alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913640,8620GA100.0% 44.8 / NA 14A370A (GCG→GCAahpFalkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/9);  total (5/9)

TGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGGCGTGGAAGCGGCAATTGACCTGGCGGGTATCGTTGAGCACGTAACGCTGCTGGAATTTGCGCCAGAAATGAAAGCCGACCAGGTTCTGCAGGACAAACTGCGCAGCCTGAAA  >  NC_000913/640728‑640964
                                                                                                                                      |                                                                                                      
tGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGGCGTGGAAGCAGCaa                                                                                                    <  2:200133/139‑1 (MQ=255)
     cAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGGCGTGGAAGCAGCAATTGAc                                                                                               <  1:182191/139‑1 (MQ=255)
             aCGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGGCGTGGAAGCAGCAATTGACCTGGCggg                                                                                       <  2:448316/139‑1 (MQ=255)
               gTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGGCGTGGAAGCAGCAATTGACCTGGCGGGTa                                                                                     <  1:185225/139‑1 (MQ=255)
                                        ccAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGGCGTGGAAGCAGCAATTGACCTGGCGGGTATCGTTGAGCACGTAACGCTGCTGGa                                                            >  2:3065/1‑139 (MQ=255)
                                                        cTGCCCGCACTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGGCGTGGAAGCAGCAATTGACCTGGCGGGTATCGTTGAGCACGTAACGCTGCTGGAATTTGCGCCAGAAATg                                            >  2:308077/1‑139 (MQ=255)
                                                                aCTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGGCGTGGAAGCAGCAATTGACCTGGCGGGTATCGtt                                                                                <  1:406852/95‑1 (MQ=255)
                                                                aCTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGGCGTGGAAGCAGCAATTGACCTGGCGGGTATCGtt                                                                                >  2:406852/1‑95 (MQ=255)
                                                                aCTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGGCGTGGAAGCAGCAATTGACCTGGCGGGTATCGTTGAGCACGTAACGCTGCTGGAATTTGCACCAGAAATGAAAGCCGa                                    <  2:313710/139‑1 (MQ=255)
                                                                       cGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGGCGTGGAAGCAGCAATTGACCTGGCGGGTATCGTTGAGCACGTAACGCTGCTGGAATTTGCGCCAGAAATGAAAGCCGACCAGGtt                             <  2:157140/139‑1 (MQ=255)
                                                                                        gTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGGCGTGGAAGCAGCAATTGACCTGGCGGGTATCGTTGAGCACGTAACGCTGCTGGAATTTGCGCCAGAAATGAAAGCCGACCAGGTTCTGCAGGACAAACTgcg            <  1:159241/139‑1 (MQ=255)
                                                                                                  aGCGGTTATCGGCGGCGGTAACTCCGGCGTGGAAGCAGCAATTGACCTGGCGGGTATCGTTGAGCACGTAACGCTGCTGGAATTTGCGCCAGAAATGAAAGCCGACCAGGTTCTGCAGGACAAACTGCGCAGCCTGaaa  >  2:76008/1‑139 (MQ=255)
                                                                                                                                  aaGCAGCAATTGACCTGGCGGGTATCGTTGAGCACGTAACGCTGCTGGAATTTGCGCCAGAAATGAAAGc                                       <  1:113991/70‑1 (MQ=255)
                                                                                                                                  aaGCAGCAATTGACCTGGCGGGTATCGTTGAGCACGTAACGCTGCTGGAATTTGCGCCAGAAATGAAAGc                                       >  2:113991/1‑70 (MQ=255)
                                                                                                                                      |                                                                                                      
TGGCGCAACATGAACGTTCCGGGCGAAGATCAGTATCGCACCAAAGGCGTGACCTACTGCCCGCACTGCGACGGCCCGCTGTTTAAAGGTAAACGCGTAGCGGTTATCGGCGGCGGTAACTCCGGCGTGGAAGCGGCAATTGACCTGGCGGGTATCGTTGAGCACGTAACGCTGCTGGAATTTGCGCCAGAAATGAAAGCCGACCAGGTTCTGCAGGACAAACTGCGCAGCCTGAAA  >  NC_000913/640728‑640964

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