Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,765,666 T→C T37T (ACA→ACG)  yfjM ← CP4‑57 prophage, uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,765,6660TC100.0% 204.1 / NA 53T37T (ACA→ACG) yfjMCP4‑57 prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (28/25);  total (28/25)

TCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCC  >  NC_000913/2765533‑2765793
                                                                                                                                     |                                                                                                                               
tCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTc                                                                                                                            >  1:300771/1‑139 (MQ=255)
    atCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACg                                                                                                                        >  2:37933/1‑139 (MQ=255)
          aTGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAAt                                                                                                                  >  2:283416/1‑139 (MQ=255)
                         tGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTctgct                                                                                                   >  1:420354/1‑139 (MQ=255)
                          gTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTctgctg                                                                                                  >  1:229901/1‑139 (MQ=255)
                             aaCATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCGCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCg                                                                                               <  1:455828/139‑1 (MQ=255)
                             aaCATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCg                                                                                               <  1:551096/139‑1 (MQ=255)
                                   gcgTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGaataat                                                                                         >  1:263538/1‑139 (MQ=255)
                                    cgTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAAtt                                                                                        >  2:390202/1‑139 (MQ=255)
                                     gtgtTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTa                                                                                       >  1:396590/1‑139 (MQ=255)
                                                              aTGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGt                                                               <  2:444763/138‑1 (MQ=255)
                                                              aTGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGt                                                               >  1:444763/1‑138 (MQ=255)
                                                              aTGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTg                                                              <  2:396590/139‑1 (MQ=255)
                                                                     gCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTctg                                                                                                     <  1:277168/93‑1 (MQ=255)
                                                                     gCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTctg                                                                                                     >  2:277168/1‑93 (MQ=255)
                                                                     gCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGaat                                                                                            <  2:826813/102‑1 (MQ=255)
                                                                     gCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGaat                                                                                            >  1:826813/1‑102 (MQ=255)
                                                                                         ttATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGAt                                   <  1:390202/139‑1 (MQ=255)
                                                                                          tataGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTgg                                                                >  2:171203/1‑109 (MQ=255)
                                                                                          tataGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTgg                                                                <  1:171203/109‑1 (MQ=255)
                                                                                           ataGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATgg                                 >  1:893436/1‑139 (MQ=255)
                                                                                                tcatcaTGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCg                                                                                               <  1:97414/72‑1 (MQ=255)
                                                                                                tcatcaTGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCg                                                                                               >  2:97414/1‑72 (MQ=255)
                                                                                                  atcatGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACt                                                                                     >  1:206447/1‑80 (MQ=255)
                                                                                                  atcatGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACt                                                                                     <  2:206447/80‑1 (MQ=255)
                                                                                                  atcatGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGCGTTGCTGCTTGGTg                                                              >  2:224707/1‑103 (MQ=255)
                                                                                                  atcatGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTg                                                              >  2:535595/1‑103 (MQ=255)
                                                                                                  atcatGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTg                                                              <  1:224707/103‑1 (MQ=255)
                                                                                                  atcatGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGAcatcat                          >  1:540427/1‑139 (MQ=255)
                                                                                                  atcatGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGAcatcat                          >  1:22107/1‑139 (MQ=255)
                                                                                                  atcatGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAACTATCGAGTTGCTGCTTGGTg                                                              <  1:535595/103‑1 (MQ=255)
                                                                                                  atcatGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGACTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGAt                                                          >  1:335301/1‑107 (MQ=255)
                                                                                                  atcatGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGACTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGAt                                                          <  2:335301/107‑1 (MQ=255)
                                                                                                   tcatGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTa                                                                                                                      >  1:710308/1‑46 (MQ=255)
                                                                                                   tcatGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTa                                                                                                                      <  2:710308/46‑1 (MQ=255)
                                                                                                    catGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAAc                                      >  2:905607/1‑125 (MQ=255)
                                                                                                    catGATGGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAAc                                      <  1:905607/125‑1 (MQ=255)
                                                                                                       gatgGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTtgctgc                                                                    <  1:722985/92‑1 (MQ=255)
                                                                                                       gatgGCTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTtgctgc                                                                    >  2:722985/1‑92 (MQ=255)
                                                                                                            cTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAaga                                                  <  2:307012/105‑1 (MQ=255)
                                                                                                            cTTGTTTTGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAaga                                                  >  1:307012/1‑105 (MQ=255)
                                                                                                                   tGGGCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACt         <  1:37933/139‑1 (MQ=255)
                                                                                                                     ggCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACt                                                                                     >  2:351903/1‑61 (MQ=255)
                                                                                                                     ggCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACt                                                                                     <  1:351903/61‑1 (MQ=255)
                                                                                                                      gCTCTCTGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcc      <  2:300771/139‑1 (MQ=255)
                                                                                                                          tctGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCg                                                                                               <  2:309636/46‑1 (MQ=255)
                                                                                                                          tctGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCg                                                                                               >  1:309636/1‑46 (MQ=255)
                                                                                                                          tctGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGaataa                                                                                          <  2:266508/51‑1 (MQ=255)
                                                                                                                          tctGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGaataa                                                                                          >  1:266508/1‑51 (MQ=255)
                                                                                                                          tctGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGaataa                                                                                          >  1:319749/1‑51 (MQ=255)
                                                                                                                          tctGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGaataa                                                                                          <  2:319749/51‑1 (MQ=255)
                                                                                                                          tctGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcccccc  <  2:420354/139‑1 (MQ=255)
                                                                                                                          tctGCTCGTTCCGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcccccc  <  1:283416/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                               
TCATATCTCGATGCTCTTTGAGCAATGTCAACATCGCGTGTTCATGGCTTTCTATATTGTTGATGCCTTGCCCATCCGGACCCCACTCCTTATAGCTCATCATGATGGCTTGTTTTGGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCC  >  NC_000913/2765533‑2765793

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: ‑