Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,420,793 (G)6→7 coding (110/423 nt) ydaT → Rac prophage, uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,420,7871.G100.0% 174.5 / NA 49Q35R (CAG→CGG) ydaTRac prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (28/21);  total (28/21)

ACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTG  >  NC_000913/1420656‑1420917
                                                                                                                                    |                                                                                                                                  
aCCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑ggggggg                                                                                                                             >  1:213873/1‑138 (MQ=255)
          tCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTg                                                                                                                    >  2:1247727/1‑139 (MQ=255)
          tCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGACAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTg                                                                                                                    >  2:994645/1‑139 (MQ=255)
           cAAGAGGCAAATGATACATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGc                                                                                                                   >  1:1078808/1‑139 (MQ=255)
              gagGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCACTTACTCAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGcccc                                                                                                                <  2:1078808/139‑1 (MQ=255)
                 gCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTgg                                                                                                             >  2:250672/1‑139 (MQ=255)
                  cAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGt                                                                                                            >  2:1215429/1‑139 (MQ=255)
                    aaTGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTc                                                                                                          <  2:213873/139‑1 (MQ=255)
                         ttCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCTAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCACTGGTTCCAGGc                                                                                                     >  1:1520549/1‑139 (MQ=255)
                          tCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCa                                                                                                    <  1:1247727/139‑1 (MQ=255)
                                   tCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGa                                                                                           <  1:722712/139‑1 (MQ=255)
                                              cacaTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcagcag                                                                                >  1:741430/1‑139 (MQ=255)
                                               acaTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGTATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGa                                                                               <  1:994645/139‑1 (MQ=255)
                                                caTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGaa                                                                              <  1:250672/139‑1 (MQ=255)
                                                caTCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGaa                                                                              <  2:1160643/139‑1 (MQ=255)
                                                  tCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAata                                                                            >  2:181151/1‑139 (MQ=255)
                                                   cGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGTAAg                                                                                                  <  1:453605/116‑1 (MQ=255)
                                                   cGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGTAAg                                                                                                  >  2:453605/1‑116 (MQ=255)
                                                           tgttgtTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcag                                                                                   <  2:137322/123‑1 (MQ=255)
                                                           tgttgtTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATcag                                                                                   >  1:137322/1‑123 (MQ=255)
                                                                tttGCCCTAGCAGCCGAAAAAGGCCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCgttgg                                                              >  1:1572678/1‑139 (MQ=255)
                                                                     ccTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGa                                                                                                >  2:13506/1‑100 (MQ=255)
                                                                     ccTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGa                                                                                                <  1:13506/100‑1 (MQ=255)
                                                                     ccTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGaa                                                         <  2:741430/139‑1 (MQ=255)
                                                                       tAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAAtt                                                                                                                     >  2:1346452/1‑77 (MQ=255)
                                                                       tAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAAtt                                                                                                                     <  1:1346452/77‑1 (MQ=255)
                                                                                       gcaggcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGGAAGGTGAAACGAAAACGCa                                       >  1:1008966/1‑139 (MQ=255)
                                                                                       gcaggcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCa                                       >  1:1271339/1‑139 (MQ=255)
                                                                                                ggTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCa                                                                                                    <  1:564507/69‑1 (MQ=255)
                                                                                                ggTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCa                                                                                                    >  2:564507/1‑69 (MQ=255)
                                                                                                     ccAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTc                                                                                                          <  2:1106237/58‑1 (MQ=255)
                                                                                                     ccAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTc                                                                                                          <  2:217979/58‑1 (MQ=255)
                                                                                                     ccAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTc                                                                                                          >  1:1106237/1‑58 (MQ=255)
                                                                                                     ccAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTc                                                                                                          >  1:217979/1‑58 (MQ=255)
                                                                                                     ccAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTc                         >  1:384619/1‑139 (MQ=255)
                                                                                                         tGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATGGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGaa                     >  2:1029659/1‑139 (MQ=255)
                                                                                                         tGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGACGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAAt                                                                                      >  2:250287/1‑74 (MQ=255)
                                                                                                         tGCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGACGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAAt                                                                                      <  1:250287/74‑1 (MQ=255)
                                                                                                          gCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGaaaaaa                            <  2:464327/131‑1 (MQ=255)
                                                                                                          gCAATTACTGAAGAATATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGaaaaaa                            >  1:464327/1‑131 (MQ=255)
                                                                                                                    aagaaTATCTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAg          >  2:492739/1‑139 (MQ=255)
                                                                                                                            cTGCGCCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTg  >  1:1437546/1‑139 (MQ=255)
                                                                                                                              gcgcCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGc                                                                                                     <  2:1026331/38‑1 (MQ=255)
                                                                                                                              gcgcCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGc                                                                                                     >  1:1026331/1‑38 (MQ=37)
                                                                                                                              gcgcCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGAc                                                                                               <  1:530713/44‑1 (MQ=255)
                                                                                                                              gcgcCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGAc                                                                                               >  2:530713/1‑44 (MQ=38)
                                                                                                                              gcgcCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCa                                       <  2:422336/100‑1 (MQ=255)
                                                                                                                              gcgcCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCa                                       >  1:422336/1‑100 (MQ=255)
                                                                                                                              gcgcCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCa                        >  1:704105/1‑115 (MQ=255)
                                                                                                                              gcgcCAGGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCa                        <  2:704105/115‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                  
ACCTGAAACATCAAGAGGCAAATGATTCATGAAAATCAAGCATGAGCACATCGAATCAGTGTTGTTTGCCCTAGCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCA‑GGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCAGCAATTCTG  >  NC_000913/1420656‑1420917

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