Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,984,411 C→T *394* (TAG→TAA yqeF ← short chain acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,984,4110CT100.0% 79.5 / NA 27*394* (TAG→TAAyqeFshort chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/11);  total (16/11)

GTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGC  >  NC_000913/2984298‑2984541
                                                                                                                 |                                                                                                                                  
gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                                <  1:111510/134‑1 (MQ=255)
gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGt                                                                                                                >  2:111510/1‑134 (MQ=255)
gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATg                                                                                                           >  1:548353/1‑139 (MQ=255)
        aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACAccc                                                                                                   >  1:471031/1‑139 (MQ=255)
          tGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                 >  1:174539/1‑139 (MQ=255)
          tGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTg                                                                                                 >  1:422767/1‑139 (MQ=255)
             aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                              <  1:583027/139‑1 (MQ=255)
                      ggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                         <  2:626129/119‑1 (MQ=255)
                      ggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                         >  1:626129/1‑119 (MQ=255)
                                                    gccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                       <  1:140118/139‑1 (MQ=255)
                                                          attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                 >  2:228767/1‑139 (MQ=255)
                                                                   cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGc                                                                                               <  2:93642/84‑1 (MQ=255)
                                                                   cAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGc                                                                                               >  1:93642/1‑84 (MQ=255)
                                                                    aGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCcac                                                                                                       <  1:199234/75‑1 (MQ=255)
                                                                    aGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCcac                                                                                                       >  2:199234/1‑75 (MQ=255)
                                                                      aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTc                                                                                                                            <  2:96/52‑1 (MQ=255)
                                                                      aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTc                                                                                                                            >  1:96/1‑52 (MQ=255)
                                                                                               tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                         >  2:565327/1‑46 (MQ=255)
                                                                                               tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGcc                                                                                                         <  1:565327/46‑1 (MQ=255)
                                                                                               tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                              >  1:207405/1‑57 (MQ=255)
                                                                                               tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGcc                                                                                              <  2:207405/57‑1 (MQ=255)
                                                                                               tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg                                                           <  2:579811/92‑1 (MQ=255)
                                                                                               tAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCaa            >  1:382022/1‑139 (MQ=255)
                                                                                               tAAGGAGGAAATGCAGAGTTATTCGTCACGATCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg                                                           >  1:579811/1‑92 (MQ=255)
                                                                                                 aggaggAAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc          >  1:369567/1‑139 (MQ=255)
                                                                                                      gaAATGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg     <  2:422767/139‑1 (MQ=255)
                                                                                                         aTGCAGAGTTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc  >  1:666003/1‑139 (MQ=255)
                                                                                                                 |                                                                                                                                  
GTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGC  >  NC_000913/2984298‑2984541

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