Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,769,486 G→A *274* (TAG→TAA)  yfjQ → CP4‑57 prophage, uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,769,4860GA100.0% 179.9 / NA 57*274* (TAG→TAA) yfjQCP4‑57 prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/27);  total (30/27)

CGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAGACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGT  >  NC_000913/2769357‑2769617
                                                                                                                                 |                                                                                                                                   
cGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCAt                                                                                                                            >  2:236194/1‑139 (MQ=255)
 ggACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATg                                                                                                                           <  1:1343578/139‑1 (MQ=255)
  gACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGt                                                                                                                          >  1:1438433/1‑139 (MQ=255)
      cAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                                      <  2:1438433/139‑1 (MQ=255)
       aGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCa                                                                                                                     <  2:1497985/139‑1 (MQ=255)
       aGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCa                                                                                                                     >  1:749729/1‑139 (MQ=255)
                            ggACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                                    <  2:1492410/119‑1 (MQ=255)
                            ggACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                                    >  1:1492410/1‑119 (MQ=255)
                            ggACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcacca                                                                                                >  2:1460193/1‑139 (MQ=255)
                             gACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcaccac                                                                                               >  1:1566482/1‑139 (MQ=255)
                             gACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcacca                                                                                                <  2:1618842/138‑1 (MQ=255)
                             gACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcacca                                                                                                >  1:1618842/1‑138 (MQ=255)
                             gACGAACTCATAATCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGCTGCCACGTTGTTAATATCGGAcaccac                                                                                               <  2:1007207/139‑1 (MQ=255)
                                             tgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTAt                                                                               >  1:1212514/1‑139 (MQ=255)
                                              gccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                                      >  1:1065781/1‑99 (MQ=255)
                                              cccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                                      <  2:1065781/98‑1 (MQ=255)
                                                        gCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGtgtg                                                                    <  2:380750/139‑1 (MQ=255)
                                                            tGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                                      >  1:298468/1‑85 (MQ=255)
                                                            tGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                                      <  2:298468/85‑1 (MQ=255)
                                                            tGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTc                                                                >  2:110529/1‑139 (MQ=255)
                                                                  ggATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATccc                                                          <  2:749729/139‑1 (MQ=255)
                                                                  ggATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATccc                                                          <  1:236194/139‑1 (MQ=255)
                                                                   gATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcacca                                                                                                >  1:892774/1‑100 (MQ=255)
                                                                   gATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcacca                                                                                                <  2:892774/100‑1 (MQ=255)
                                                                     tgtgAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCa                                                                                      <  1:742473/108‑1 (MQ=255)
                                                                     tgtgAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCa                                                                                      >  2:742473/1‑108 (MQ=255)
                                                                      gtgAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcacca                                                                                                >  2:506055/1‑97 (MQ=255)
                                                                      gtgAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcacca                                                                                                <  1:506055/97‑1 (MQ=255)
                                                                                 ccGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCtttt                                           >  2:1544192/1‑139 (MQ=255)
                                                                                    tGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                                      >  1:483755/1‑61 (MQ=255)
                                                                                    tGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                                      <  2:483755/61‑1 (MQ=255)
                                                                                    tGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAAc                                        >  1:1441002/1‑139 (MQ=255)
                                                                                      cgcTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACgttgt                                                                                                               >  2:559711/1‑66 (MQ=255)
                                                                                      cgcTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACgttgt                                                                                                               <  1:559711/66‑1 (MQ=255)
                                                                                      cgcTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAAc                                        >  1:1640615/1‑137 (MQ=255)
                                                                                      cgcTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAAc                                        <  2:1640615/137‑1 (MQ=255)
                                                                                          gtgGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGt                                  <  1:110529/139‑1 (MQ=255)
                                                                                              gTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                                   >  2:1415149/1‑54 (MQ=255)
                                                                                              gTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                                   <  1:1415149/54‑1 (MQ=255)
                                                                                               tGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACAt                             <  1:1544192/139‑1 (MQ=255)
                                                                                                gATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                                      <  1:473443/49‑1 (MQ=255)
                                                                                                gATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                                      >  2:473443/1‑49 (MQ=255)
                                                                                                 aTGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGt                                                                                           <  2:953875/75‑1 (MQ=255)
                                                                                                 aTGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGt                                                                                           >  1:953875/1‑75 (MQ=255)
                                                                                                   ggCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAatat                                                                                                         >  1:513111/1‑59 (MQ=255)
                                                                                                   ggCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAatat                                                                                                         <  2:513111/59‑1 (MQ=255)
                                                                                                   ggCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCt                                                                                 >  1:147707/1‑83 (MQ=255)
                                                                                                   ggCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCt                                                                                 <  2:147707/83‑1 (MQ=255)
                                                                                                      aGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTa                                          <  2:1269398/119‑1 (MQ=255)
                                                                                                      aGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTa                                          >  1:1269398/1‑119 (MQ=255)
                                                                                                              tGCTCACGCAACTGCAGTAAACGTTTCATGTTGCGAc                                                                                                                    <  1:781752/37‑1 (MQ=17)
                                                                                                              tGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                                    >  2:781752/1‑37 (MQ=255)
                                                                                                                    cGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCg                                                                                        >  2:226065/1‑59 (MQ=255)
                                                                                                                    cGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCg                                                                                        <  1:226065/59‑1 (MQ=255)
                                                                                                                        aCTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATa    >  1:760027/1‑139 (MQ=255)
                                                                                                                          tGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGt  >  2:622032/1‑139 (MQ=255)
                                                                                                                          tGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGt  >  2:924202/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                   
CGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAGACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTGCAGAAATAAAAAATAGT  >  NC_000913/2769357‑2769617

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