Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,649,041 G→A *73* (TAG→TAA ydfC → uncharacterized protein, Qin prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,649,0410GA100.0% 94.7 / NA 30*73* (TAG→TAAydfCuncharacterized protein, Qin prophage
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/14);  total (16/14)

CTGAAATGGGCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCAT  >  NC_000913/1648916‑1649175
                                                                                                                             |                                                                                                                                      
cTGAAATGGGCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATcaaca                                                                                                                           <  1:81759/139‑1 (MQ=255)
  gAAATGGGCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACaaa                                                                                                                         >  1:175344/1‑139 (MQ=255)
     aTGGGCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCt                                                                                                                      <  1:573395/139‑1 (MQ=255)
       gggCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTaa                                                                                                                    >  2:317392/1‑139 (MQ=255)
        ggCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAg                                                                                                                   >  2:145289/1‑139 (MQ=255)
             aTGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAg                                                                                                              >  2:522180/1‑139 (MQ=255)
                                ggCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCAc                                                                                           <  1:145289/139‑1 (MQ=255)
                                            gTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGa                                                                               <  1:317392/139‑1 (MQ=255)
                                                     ggAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGtgcgtgcg                                                                      <  1:522180/139‑1 (MQ=255)
                                                     ggAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGtgcgtgcg                                                                      <  1:329570/139‑1 (MQ=255)
                                                            gATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGaca                                                               <  2:175344/139‑1 (MQ=255)
                                                                      atatcgAAAAAACAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACATAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCa                                                     >  1:586945/6‑139 (MQ=255)
                                                                      aTCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATcaac                                                                                                                            <  1:567852/68‑1 (MQ=255)
                                                                      aTCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATcaac                                                                                                                            >  2:567852/1‑68 (MQ=255)
                                                                      aTCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGaca                                                               >  2:739547/1‑129 (MQ=255)
                                                                      aTCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGaca                                                               <  1:739547/129‑1 (MQ=255)
                                                                      aTCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCa                                                     >  2:851559/1‑139 (MQ=255)
                                                                      aTCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCa                                                     >  2:948385/1‑139 (MQ=255)
                                                                       tCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATcaaca                                                                                                                           <  2:143183/68‑1 (MQ=255)
                                                                       tCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATcaaca                                                                                                                           >  1:143183/1‑68 (MQ=255)
                                                                            aaaaaGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGa                                               >  2:417466/1‑139 (MQ=255)
                                                                                     cGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGt                                      >  1:348489/1‑139 (MQ=255)
                                                                                                  tGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTa                                                                                                                     <  1:880550/47‑1 (MQ=255)
                                                                                                  tGCCAGAATGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTa                                                                                                                     >  2:880550/1‑47 (MQ=255)
                                                                                                        aaTGATTACATACCCGCGTTAATAGTTAATCAACAAAGCTAAGGtt                                                                                                                <  1:594836/46‑1 (MQ=255)
                                                                                                        aaTGATTACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGtt                                                                                                                >  2:594836/1‑46 (MQ=255)
                                                                                                             ttACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGa                                               >  2:249086/1‑106 (MQ=255)
                                                                                                             ttACATACCCGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGa                                               <  1:249086/106‑1 (MQ=255)
                                                                                                                     ccGCGTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGa      >  2:388148/1‑139 (MQ=255)
                                                                                                                         gTTAATAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGGCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCAt  <  2:292485/139‑1 (MQ=255)
                                                                                                                             |                                                                                                                                      
CTGAAATGGGCTAATGAAGCCAGGGCTAAGCAGGCTGAAAATCTGTCCGCTCTGGAGCAGGATATCACCAATCTGGAAAAAGCAGCGGCATTAAGCATTGCCAGAATGATTACATACCCGCGTTAGTAGCTAATCAACAAAGCTAAGGTTAGTAATTAAGGAGTTCTCCACGGGTGAGGTGGAGTGCGTGCGCCGGACACGGGTGCGCATCCGGAACTGACAGTTTACTGAAAGGATATTTCCCTGAAAAGTCAGACCAT  >  NC_000913/1648916‑1649175

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