Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,409,177 G→A G91G (GGC→GGT vacJ ← predicted lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,409,1770GA100.0% 79.3 / NA 22G91G (GGC→GGTvacJpredicted lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/14);  total (8/14)

AACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCGCCCTGCAAGAAGTAGTTAACCATCACCGCAGGTTCTTCAAGGTTGCCAGTAAAGTTGCTCAAACCGTTACGC  >  REL606/2409110‑2409249
                                                                   |                                                                        
naCCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCACCCTGCaa                                                                  <  1:673345/75‑1 (MQ=255)
 aCCGCGCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCACCCTGCaag                                                                 <  2:620815/76‑1 (MQ=255)
       cATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCACCCTGCAAGAagtag                                                           >  1:338833/1‑76 (MQ=255)
          ccccAAAATGGTGTTTAGGAAAAAGAGGGAAAAGTGGACAATCCCCTGATAAGGGTGACCCTGCAAGTAGTAGTTa                                                        >  1:417935/1‑76 (MQ=255)
            ccAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCACCCTGCAAGAAGCAGTTAAc                                                      <  1:491642/76‑1 (MQ=255)
             cAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCACCCTGCAAGAAGTAGTTAAcc                                                     >  2:111028/1‑76 (MQ=255)
              aaaaTGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCACCCTGCAAGAAGTAGTTAACCa                                                    >  2:539390/1‑76 (MQ=255)
                aaTGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCACCCTGCAAGAAGTAGTTAACCATc                                                  >  2:501338/1‑76 (MQ=255)
                              aaaaGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCACCCTGCAAGAAGTAGTTAACCATCACCGCAGGttcttc                                    <  2:757253/76‑1 (MQ=255)
                                aaGCGGGTAAAGTGGACCTTCCCCTGATAAGGGTCACCCTGCAAGAAGTAGTTAACCATCACCGCAGGTTCTTCaa                                  >  2:711764/1‑76 (MQ=255)
                                 aGCGGGTAAAGTGGACCATCTCCTGATAAGGGTCACCCTGCAAGAAGTAGTTAACCATCACCGCAGGTTCTTCAAg                                 >  1:548569/1‑76 (MQ=255)
                                 aGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCACCCTGCAAGAAGTAGTTAACCATCACCGCAGGTTCTTCAAg                                 <  1:501338/76‑1 (MQ=255)
                                 aGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCACCCTGCAAGAAGTAGTTAACCATCACCGCAGGTTCTTCAAg                                 >  1:290448/1‑76 (MQ=255)
                                    gggTAAAGTGGACCATCCCCTGATAAGGGTCACCCTGCAAGAAGTAGTTAACCATCACCGCAGGTTCTTCAAGGtt                              <  1:539390/76‑1 (MQ=255)
                                    gggTAAAGTGGACCATCCCCTGATAAGGGTCACCCTGCAAGAAGTAGTTAACCATCACCGCAGGTTCTTCAAGGtt                              <  2:402769/76‑1 (MQ=255)
                                            nggNCNACCCCCTTNTACGGGGTGCCGTGCAANAAGTAGTTCACCGTCACCGCAGGTTCCTCAAGGTTGCCAGTaa                      <  1:664751/75‑1 (MQ=255)
                                               aCCATCCCCTGATAAGGGTCACCCTGCAAGTAGTAGTTAACCATCACCGCAGGTTCTTCAAGGTTGCCAGTAAAGt                   <  1:475162/76‑1 (MQ=255)
                                                      ccTGATAAGGGTCACCCTGCAAGAAGTAGTTAACCATCACCGCAGGTTCTTCAAGGTTGCCAGTAAAGTTGCTCaa            >  1:419598/1‑76 (MQ=255)
                                                         gATAAGGGTCACCCTGCAAGAAGTAGTTAACCATCACCGCAGGTTCTTCAAGGTTGCCAGTAAAGTTGCTCAAAcc         <  1:333765/76‑1 (MQ=255)
                                                            aaGGGTCACCCTGCAAGAAGTAGTTAACCATCACCGCAGGTTCTTCAAGGTTGCCAGTAAAGTTGCTCAAACCGtt      <  2:460610/76‑1 (MQ=255)
                                                             aGGGTCACCCTGCAAGAAGTAGTTAACCATCACCGCAGGTTCTTCAAGGTTGCCAGTAAAGTTGCTCAAACCGTTa     <  2:349033/76‑1 (MQ=255)
                                                             aGGGTCACCCTGCAAGAAGTAGTTAACCATCACCGCAGGTTCTTCAAGGTTGCCAGTAAAGTTGCTCAAACCGTTa     <  2:74558/76‑1 (MQ=255)
                                                               ggTCACCCTGCAAGAAGTAGTTAACCATCACCGCAGGTTCTTCAAGGTTGCCAGTAAAGTTGCACAAACCGTTAcg   <  2:6155/76‑1 (MQ=255)
                                                                gTCACCCTGCAAGAAGTAGTTAACCATCACCGCTGGTTCTTCAAGGTTGCCAGTAAAGTTGCTCAAACGGTTAcgc  <  1:517079/76‑1 (MQ=255)
                                                                   |                                                                        
AACCGCCCATCCCCAAAATGGTGTTCAGGAAAAAGCGGGTAAAGTGGACCATCCCCTGATAAGGGTCGCCCTGCAAGAAGTAGTTAACCATCACCGCAGGTTCTTCAAGGTTGCCAGTAAAGTTGCTCAAACCGTTACGC  >  REL606/2409110‑2409249

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