Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,421,856 C→T S246S (AGC→AGT ydaU → Rac prophage; DUF1376 domain‑containing protein YdaU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,421,8560CT100.0% 147.4 / NA 43S246S (AGC→AGTydaURac prophage; DUF1376 domain‑containing protein YdaU
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/25);  total (18/25)

GCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGCATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAAT  >  NC_000913/1421718‑1421991
                                                                                                                                          |                                                                                                                                       
gCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTCACCTGAGCAGGTGATTCGACACAGTAt                                                                                                                                       <  2:972733/140‑1 (MQ=255)
   ttGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCg                                                                                                                                     <  2:506446/140‑1 (MQ=255)
    tGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGc                                                                                                                                    <  2:207031/140‑1 (MQ=255)
        gTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAAt                                                                                                                                <  1:938883/140‑1 (MQ=255)
         tGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATg                                                                                                                               <  1:399901/140‑1 (MQ=255)
          gAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGa                                                                                                                              <  1:250414/140‑1 (MQ=255)
           aGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGc                                                                                                                                    <  1:706628/133‑1 (MQ=255)
           aGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGc                                                                                                                                    >  2:706628/1‑133 (MQ=255)
           aGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGaa                                                                                                                             >  2:542170/1‑140 (MQ=255)
              gggTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATAc                                                                                                                          <  2:844726/140‑1 (MQ=255)
                gTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCa                                                                                                                        >  2:314654/1‑140 (MQ=255)
                     aTTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCg                                                                                                                                     <  1:196343/122‑1 (MQ=255)
                     aTTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCg                                                                                                                                     >  2:196343/1‑122 (MQ=255)
                        tCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGt                                                                                                                <  1:241942/140‑1 (MQ=255)
                          gCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTc                                                                                                              <  1:139029/140‑1 (MQ=255)
                          gCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTc                                                                                                              <  1:180336/140‑1 (MQ=255)
                               gCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGcc                                                                                                         <  1:277575/140‑1 (MQ=255)
                                  ttGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGaa                                                                                                      <  2:389752/140‑1 (MQ=255)
                                   tGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGaaa                                                                                                     <  1:697616/140‑1 (MQ=255)
                                   tGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGaaa                                                                                                     <  1:923789/140‑1 (MQ=255)
                                                    aCGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGAc                                                                                    <  1:145634/140‑1 (MQ=255)
                                                    aCGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGAc                                                                                    <  1:570730/140‑1 (MQ=255)
                                                    aCGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGAc                                                                                    <  1:352680/140‑1 (MQ=255)
                                                             ggCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCg                                                                           >  1:751/1‑140 (MQ=255)
                                                                  tAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAgg                                                                      >  1:626801/1‑140 (MQ=255)
                                                                  tAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAgg                                                                      >  1:28089/1‑140 (MQ=255)
                                                                  tAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAgg                                                                      >  2:2992/1‑140 (MQ=255)
                                                                  tAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAgg                                                                      <  2:844363/140‑1 (MQ=255)
                                                                  tAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAgg                                                                      >  1:36348/1‑140 (MQ=255)
                                                                        ggCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAaca                                                                <  1:320404/140‑1 (MQ=255)
                                                                        ggCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAaca                                                                >  2:683165/1‑140 (MQ=255)
                                                                            aTACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGt                                                            >  2:357188/1‑140 (MQ=255)
                                                                                  gAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCTCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGtt                                                      >  2:869515/1‑140 (MQ=255)
                                                                                  gAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGtt                                                      >  2:435312/1‑140 (MQ=255)
                                                                                  gAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGtt                                                      >  2:333210/1‑140 (MQ=255)
                                                                                    gCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTAc                                                    <  1:435312/140‑1 (MQ=255)
                                                                                    gCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTAc                                                    <  2:852042/140‑1 (MQ=255)
                                                                                              ttCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTAGCGGATAGTGc                                          >  2:613173/1‑140 (MQ=255)
                                                                                                ccGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGa                                        <  1:542170/140‑1 (MQ=255)
                                                                                                             ttttCACCTGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCtt                           <  1:357188/140‑1 (MQ=255)
                                                                                                                     tGAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGt                   >  2:993549/1‑140 (MQ=255)
                                                                                                                      gAGCAGGTGATTCGACACAGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTa                  >  2:21483/1‑140 (MQ=255)
                                                                                                                                      acaGTATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAat  >  1:691109/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                       
GCTTTGGCGTGAGTGGGTTGAATTTCGCCAGGCATTGCGAAAACCGATTCGAACGGAGCAGGGCGCTAACGGGGCGATACGGGAGCTGGAAAAATTCCGCCAGCAGGGTTTTTCACCTGAGCAGGTGATTCGACACAGCATCGCCAATGAATACCAGGGCTTGTTCGCGCCGAAAGGTGTTCGACCTGAGACGTTACTCCGACAGGTTAACACCGTCTCGTTACCGGATAGTGCGATCCCGCCAGGCTTCAGGGGGTAACTGACCATGAAAAAT  >  NC_000913/1421718‑1421991

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