Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 277,996 A→G G269G (GGT→GGC afuC ← CP4‑6 prophage; putative ferric transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913277,9960AG100.0% 127.9 / NA 35G269G (GGT→GGCafuCCP4‑6 prophage; putative ferric transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/19);  total (16/19)

CTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGAGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCACCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCGGGAACAGGTTGGCATCGCCCATAAAGCTCGCCATAAAGCGGGAGGCGGGCTGGCGGTAAAGATCCTGC  >  NC_000913/277849‑278140
                                                                                                                                                   |                                                                                                                                                
cTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGccc                                                                                                                                               >  1:98323/1‑151 (MQ=255)
   cAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCtgt                                                                                                                                            <  1:265317/151‑1 (MQ=255)
   cAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCtgt                                                                                                                                            >  1:333350/1‑151 (MQ=255)
       aataTCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTAc                                                                                                                                        >  2:534676/1‑151 (MQ=255)
                       gCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGgcg                                                                                                                        >  1:283822/1‑151 (MQ=255)
                       gCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGgc                                                                                                                         <  2:97286/150‑1 (MQ=255)
                            tCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGc                                                                                                                    >  2:36591/1‑150 (MQ=255)
                               aCCGTCACTTCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGc                                                                                                                    >  1:177288/1‑147 (MQ=255)
                                      cTTCATACTGCGGCCCCATATATGCGACATGGCGGGTCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACCCCGACCATCCCATCGCCCTGTGTACCAAAGTGCAGCGGGCGCTGCAGAGGATa                                                                                                            <  2:610357/148‑1 (MQ=255)
                                        tCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCg                                                                                                        <  1:493971/150‑1 (MQ=255)
                                          aTACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAg                                                                                                     <  2:133143/151‑1 (MQ=255)
                                              tGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGatat                                                                                                 >  2:130724/1‑151 (MQ=255)
                                                  gCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAAc                                                                                             <  1:213766/151‑1 (MQ=255)
                                                        tataGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATcc                                                                                       >  2:50380/1‑151 (MQ=255)
                                                          taGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGt                                                                                     >  2:98569/1‑151 (MQ=255)
                                                              cGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTATCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCt                                                                                 <  1:226646/151‑1 (MQ=255)
                                                               gACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTg                                                                                <  2:338362/151‑1 (MQ=255)
                                                                  aTGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAg                                                                             <  2:65459/151‑1 (MQ=255)
                                                                  aTGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAg                                                                             <  1:176271/151‑1 (MQ=255)
                                                                  aTGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAg                                                                             <  1:406005/151‑1 (MQ=255)
                                                                        gATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGcc                                                                       <  1:98569/151‑1 (MQ=255)
                                                                         aTCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCg                                                                      >  1:615280/1‑151 (MQ=255)
                                                                            aCGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATGTCAACGTATCCGTCGCTGAAGGTTGCCGGGa                                                                   <  1:572949/151‑1 (MQ=255)
                                                                              gCAGCGCTGGCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCGGGAAc                                                                 <  1:8124/151‑1 (MQ=255)
                                                                                      ggCTCTCTTCGCCGCGATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCGGGAACAGGTTg                                                           <  2:572740/149‑1 (MQ=255)
                                                                                                   gcgATCGCTGGGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCGGGAACAGGTTGGCATCGCCCATaa                                              >  1:73038/1‑149 (MQ=255)
                                                                                                                  gATCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCGGGAACAGGTTGGCATCGCCCATAAAGCTCGCCATAAAGCg                              <  2:121813/150‑1 (MQ=255)
                                                                                                                    tCGCTTCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCGGGAACAGGTTGGCATCGCCCATAAAGCTCGCCATAAAGCggg                            >  2:382106/1‑150 (MQ=255)
                                                                                                                         tCCGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCGGGAACAGGTTGGCATCGCCCATAAAGCTCGCCATAAAGCGGGAGGCgg                      >  1:206044/1‑151 (MQ=255)
                                                                                                                          ccGGGCGCACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCGGGAACAGGTTGGCATCGCCCATAAAGCTCGCCATAAAGCGGGAGGCggg                     >  2:174494/1‑151 (MQ=255)
                                                                                                                              gcgcACACCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCGGGAACAGGTTGGCATCGCCCATAAAGCTCGCCATAAAGCGGGAGGCGGGCTgg                 >  2:64184/1‑151 (MQ=255)
                                                                                                                                   cacCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCGGGAACAGGTTGGCATCGCCCATAAAGCTCGCCATAAAGCGGGAGGCGGGCTGGCGGTa            <  2:192151/151‑1 (MQ=255)
                                                                                                                                   cacCGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCGGGAACAGGTTGGCATCGCCCATAAAGCTCGCCATAAAGCGGGAGGCGGGCTGGCGGTa            <  2:631390/151‑1 (MQ=255)
                                                                                                                                      cGACCATCCCTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCGGGAACAGGTTGGCATCGCCCATAAAGCTCGCCATAAAGCGGGAGGCGGGCTGGCGGTaa           >  1:357170/1‑149 (MQ=255)
                                                                                                                                             cccTTCGCCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCGGGAACAGGTTGGCATCGCCCATAAAGCTCGCCATAAAGCGGGAGGCGGGCTGGCGGTAAAGATCCTGc  <  1:130727/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CTGCAATAATATCTCCTGCCCGTGCCATTCCACCGTCACTTCATACTGCGGCCCCATATAGGCGACATGGCGGATCACGCAGCGCTGGCTCTCTTCGCCGCGATCGCTGAGCGTGATCGCTTCCGGGCGCACACCGACCATCCCTTCACCCTGTGTACCAAAGTGCAGCGGGCGCGGCAGATGATAGCCGTAGATATCAACGTATCCGTCGCTGAAGGTTGCCGGGAACAGGTTGGCATCGCCCATAAAGCTCGCCATAAAGCGGGAGGCGGGCTGGCGGTAAAGATCCTGC  >  NC_000913/277849‑278140

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