Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,208,003 A→G N168S (AAT→AGT)  ycfK → e14 prophage; protein StfP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,208,0030AG100.0% 155.9 / NA 44N168S (AAT→AGT) ycfKe14 prophage; protein StfP
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/26);  total (18/26)

CAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACGGGTAGTTATTACTATAACTGGATTGCTATTGGGTATTAAGA  >  NC_000913/1207885‑1208132
                                                                                                                      |                                                                                                                                 
cAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAGTAATGCCTGTGTAGCCGTTGTTGCAAGTTAAGACAATGTCa                                                                                                                     >  2:205945/1‑133 (MQ=255)
  aTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAgc                                                                                                                   >  1:425453/1‑133 (MQ=255)
   tCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAGTGCCTGTGTTGCCGTTTTTGCAAGTTATGACAATGTCAgcg                                                                                                                  >  1:327429/1‑133 (MQ=255)
   tCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAgcg                                                                                                                  >  1:55359/1‑133 (MQ=255)
   tCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAgcg                                                                                                                  >  1:379386/1‑133 (MQ=255)
     aGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCa                                                                                                                <  2:563906/133‑1 (MQ=255)
        aaTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGTAAGTTATGACAATGTCAGCGCACCt                                                                                                             >  2:118030/1‑133 (MQ=255)
              tCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCAcc                                                                                                              <  1:415364/126‑1 (MQ=255)
              tCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCAcc                                                                                                              >  2:415364/1‑126 (MQ=255)
              tCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAAt                                                                                                       >  1:4801/1‑133 (MQ=255)
                    ggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGt                                                                                                  >  1:316993/1‑132 (MQ=255)
                    ggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGt                                                                                                  <  2:316993/132‑1 (MQ=255)
                    ggggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTg                                                                                                 >  1:422037/1‑133 (MQ=255)
                      ggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGc                                                                                               <  2:497052/133‑1 (MQ=255)
                       gTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCa                                                                                              >  2:534522/1‑133 (MQ=255)
                          gAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACg                                                                                           <  2:4801/133‑1 (MQ=255)
                          gAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACg                                                                                           <  2:13949/133‑1 (MQ=255)
                            aTGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggg                                                                                         >  1:486294/1‑133 (MQ=255)
                             tGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACgggg                                                                                        <  1:368586/133‑1 (MQ=255)
                             tGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACgggg                                                                                        <  1:393432/133‑1 (MQ=255)
                               ccGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACgggggg                                                                                      >  2:569886/1‑133 (MQ=255)
                                cGCAAAAGGCTCTGCCTTTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAGTGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGa                                                                                     <  2:327429/133‑1 (MQ=255)
                                cGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGa                                                                                     >  1:118439/1‑133 (MQ=255)
                                cGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGa                                                                                     >  2:562032/1‑133 (MQ=255)
                                    aaaGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGAtata                                                                                 <  2:121374/133‑1 (MQ=255)
                                       gCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACaa                                                                              <  2:486294/133‑1 (MQ=255)
                                            tGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACt                                                                         <  1:569886/133‑1 (MQ=255)
                                              ccTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTTCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTc                                                                       <  2:422037/133‑1 (MQ=255)
                                                               aaTACAGTTCTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTgcg                                                      <  2:379386/133‑1 (MQ=255)
                                                                  aCAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAg                                                   <  1:562032/133‑1 (MQ=255)
                                                                         tATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTGTATGACGGTGCGCGGCTCAAAc                                            >  1:194735/1‑133 (MQ=255)
                                                                           ttCCCGGTTGCGTTTGATAATGTCTGTGTTACCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACgg                                          <  1:200902/133‑1 (MQ=255)
                                                                           ttCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACgg                                          <  1:587379/133‑1 (MQ=255)
                                                                            tCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACggg                                         <  1:534522/133‑1 (MQ=255)
                                                                                   tGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTGTTATTACGGTGCGCAGCTCAAACGGGTAGttat                                  >  2:410153/1‑133 (MQ=255)
                                                                                    gCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACGGGTAGttatt                                 <  2:118439/133‑1 (MQ=255)
                                                                                     cGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACGGGTAGttatta                                <  1:597786/133‑1 (MQ=255)
                                                                                          aataatGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACGGGTAGTTATTACtata                           >  1:307748/1‑133 (MQ=255)
                                                                                           ataatGCCTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACGGGTAGTTATTACTATaa                          <  2:307748/133‑1 (MQ=255)
                                                                                           ataatGCCTGCGTTGCCGTTGTTGTAAGCTATGACAATGTCGGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACGGGTAGTTATTACTATaa                          <  2:194735/133‑1 (MQ=255)
                                                                                                 ccTGTGTTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACGGGTAGTTATTACTATAACTGGAt                    <  2:425453/133‑1 (MQ=255)
                                                                                                     tgtTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACGGGTAGTTATTACTATAACTGGAt                    >  1:46451/1‑129 (MQ=255)
                                                                                                     tgtTGCCGTTGTTGCAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACGGGTAGTTATTACTATAACTGGAt                    <  2:46451/129‑1 (MQ=255)
                                                                                                                   cAAGTTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACGGGTAGTTATTACTATAACTGGATTGCTATTGGGTATTAAGa  <  2:55359/133‑1 (MQ=255)
                                                                                                                   cAAGTTATGACAATGTCACCGCTCCTATCATTGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACGGGTAGTTATTACTATAACTGGATTGCTATTGGGTATTAAGa  <  1:296187/133‑1 (MQ=255)
                                                                                                                      |                                                                                                                                 
CAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGCAGCTCAAACGGGTAGTTATTACTATAACTGGATTGCTATTGGGTATTAAGA  >  NC_000913/1207885‑1208132

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