Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,644,142 G→A Q55* (CAG→TAG)  ydfU ← Qin prophage, uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,644,1420GA100.0% 167.8 / NA 48Q55* (CAG→TAG) ydfUQin prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (25/23);  total (25/23)

ATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTGGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGG  >  NC_000913/1644008‑1644274
                                                                                                                                      |                                                                                                                                    
aTCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCg                                                                                                                                   <  2:1143373/138‑1 (MQ=255)
aTCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCg                                                                                                                                   >  1:1143373/1‑138 (MQ=255)
aTCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGa                                                                                                                                  >  1:681446/1‑139 (MQ=255)
  cGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAAc                                                                                                                                >  2:198268/1‑139 (MQ=255)
       ttAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAg                                                                                                                           >  2:507020/1‑139 (MQ=255)
       ttAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGAGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAg                                                                                                                           >  1:28323/1‑139 (MQ=255)
               cTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACt                                                                                                                   >  1:54825/1‑139 (MQ=255)
                            aCATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCgg                                                                                                               >  2:670233/1‑130 (MQ=255)
                            aCATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCgg                                                                                                               <  1:670233/130‑1 (MQ=255)
                            aCATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGAt                                                                                                      <  2:28323/139‑1 (MQ=255)
                                  cGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCATGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCa                                                                                                >  2:768893/1‑139 (MQ=255)
                                         gCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCtttt                                                                                         <  1:569539/139‑1 (MQ=255)
                                            cggcGGACAATGCGCAGACGCCGCGAGCGGCACGAATCGCCCGGTCGGTACTAAAAAACGGCAGCAACGTATAGACTTCTGCCGGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTgg                                                                                      >  1:786452/4‑139 (MQ=255)
                                             gTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTggt                                                                                     >  1:105340/1‑139 (MQ=255)
                                             gTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTggt                                                                                     >  2:1098442/1‑139 (MQ=255)
                                                      tGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTggttccggttcc                                                                            >  1:426458/1‑139 (MQ=255)
                                                        cgcCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTggttccggttcc                                                                            >  1:587018/1‑137 (MQ=255)
                                                        cgcCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTggttccggttcc                                                                            <  2:587018/137‑1 (MQ=255)
                                                         gcCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTccacc                                                                         <  1:507020/139‑1 (MQ=255)
                                                           cAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAg                                                                       <  2:54825/139‑1 (MQ=255)
                                                           cAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAg                                                                       <  1:1029829/139‑1 (MQ=255)
                                                              cgccgcCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGaa                                                                                                                                 <  1:868419/78‑1 (MQ=255)
                                                              cgccgcCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGaa                                                                                                                                 >  2:868419/1‑78 (MQ=255)
                                                                ccgccAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTc                                                                  <  1:1098442/139‑1 (MQ=255)
                                                                  gccAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGa                                                                <  2:681446/139‑1 (MQ=255)
                                                                                    ccGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTggtt                                                                                    <  1:346963/101‑1 (MQ=255)
                                                                                    ccGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGt                                              <  2:902514/139‑1 (MQ=255)
                                                                                    ccGTTCGTTACTAAAAAACGGCAGCAATGTTTAGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGCCTCCGGATGGCAGATTACGCATGCTTTTTggtt                                                                                    >  2:346963/1‑101 (MQ=255)
                                                                                       ttcgttACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCtt                                                                                           >  1:227164/1‑91 (MQ=255)
                                                                                       ttcgttACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCtt                                                                                           <  2:227164/91‑1 (MQ=255)
                                                                                       ttcgttACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTccacc                                                                         >  2:676387/1‑109 (MQ=255)
                                                                                       ttcgttACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTccacc                                                                         <  1:676387/109‑1 (MQ=255)
                                                                                       ttcgttACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAg                                                               >  1:1149422/1‑119 (MQ=255)
                                                                                       ttcgttACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAg                                                               <  2:1149422/119‑1 (MQ=255)
                                                                                              cTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAAc                                                                                                                       >  2:990013/1‑56 (MQ=255)
                                                                                              cTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAAc                                                                                                                       <  1:990013/56‑1 (MQ=255)
                                                                                                 aaaaaCGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGc                                 <  2:105340/139‑1 (MQ=255)
                                                                                                 aaaaaCGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGc                                 >  2:651169/1‑139 (MQ=255)
                                                                                                 aaaaaCGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGc                                 <  1:296945/139‑1 (MQ=255)
                                                                                                  aaaaCGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGcc                                <  1:616951/139‑1 (MQ=255)
                                                                                                         cagcaATGTTTTGTCTTCCGACAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTa                         >  1:693757/1‑139 (MQ=255)
                                                                                                                      tCTTCCGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGcc            >  1:988912/1‑139 (MQ=255)
                                                                                                                           cGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGa                                                                >  1:745424/1‑82 (MQ=255)
                                                                                                                           cGCCAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGa                                                                <  2:745424/82‑1 (MQ=255)
                                                                                                                              cAGCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCa    <  1:656252/139‑1 (MQ=255)
                                                                                                                                gCGGCTAGCGAACGGCAGGAATGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATg                                                        <  2:1122438/85‑1 (MQ=255)
                                                                                                                                gCGGCTAGCGAACGGCAGGAATGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATg                                                        >  1:1122438/1‑85 (MQ=255)
                                                                                                                                gCGGCTAGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAgg  >  2:377973/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                    
ATCGCCATTAGGCCACTGGCAGGATGTAACATGACGTAATAGCCAGTCGGACAATGCGCCAACGCCGCCAGCAGCACGAATCACCCGTTCGTTACTAAAAAACGGCAGCAATGTTTTGTCTTCCGCCAGCGGCTGGCGAACGGCAGGAACGACTCCGGATGGCAGATTACGCATGCTTTTTGGTTCCGGTTCCACCAGCACTCGAGGATTATGAAATATCTGTATGGATTCACGGCCCGGCTTAAGGACCACCAGCCCAAGCTCAGG  >  NC_000913/1644008‑1644274

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