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% 66.8 / NA 22*274* (TAG→TAAyfjQCP4‑57 prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/11);  total (11/11)

GGGCGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAGACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCT  >  NC_000913/2769354‑2769596
                                                                                                                                    |                                                                                                              
gggCGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGttt                                                                                                          >  1:10673/1‑139 (MQ=255)
           gTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                               <  2:110890/139‑1 (MQ=255)
                                gACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcaccac                                                                          <  1:563213/139‑1 (MQ=255)
                                                     tgcgtgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAatat                                                                                    <  1:240158/108‑1 (MQ=255)
                                                     tgcgtgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAatat                                                                                    >  2:240158/1‑108 (MQ=255)
                                                                                 tAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGacac                                                                             >  1:512020/1‑87 (MQ=255)
                                                                                 tAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGacac                                                                             <  2:512020/87‑1 (MQ=255)
                                                                                    ccGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCtttt                      >  1:792642/1‑139 (MQ=255)
                                                                                        gcgcTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAAcc                  <  2:792642/139‑1 (MQ=255)
                                                                                            tgtgGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                               >  1:459543/1‑58 (MQ=255)
                                                                                            tgtgGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                               <  2:459543/58‑1 (MQ=255)
                                                                                            tgtgGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGt                                                                      >  1:691441/1‑83 (MQ=255)
                                                                                            tgtgGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGt                                                                      <  2:691441/83‑1 (MQ=255)
                                                                                                   gATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                               >  1:330113/1‑51 (MQ=255)
                                                                                                   gATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                               <  2:330113/51‑1 (MQ=255)
                                                                                                        cAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCt  >  2:294579/1‑139 (MQ=255)
                                                                                                          gAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTg                                                                                                   <  1:674714/40‑1 (MQ=255)
                                                                                                          gAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTg                                                                                                   >  2:674714/1‑40 (MQ=25)
                                                                                                                   cTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCAt                                                                <  1:348861/66‑1 (MQ=255)
                                                                                                                   cTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCAt                                                                >  2:348861/1‑66 (MQ=255)
                                                                                                                               cAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGa                                                                                <  2:61102/38‑1 (MQ=255)
                                                                                                                               cAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGa                                                                                >  1:61102/1‑38 (MQ=255)
                                                                                                                                    |                                                                                                              
GGGCGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAGACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCT  >  NC_000913/2769354‑2769596

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