Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 286,882 G→T E32* (GAG→TAG)  yagH → CP4‑6 prophage; putative xylosidase/arabinosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913286,8820GT100.0% 147.0 / NA 42E32* (GAG→TAG) yagHCP4‑6 prophage; putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (17/25);  total (17/25)

CGCGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCT  >  NC_000913/286756‑287018
                                                                                                                              |                                                                                                                                        
cgcgACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCggg                                                                                                                             >  2:124580/1‑140 (MQ=255)
    aCAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTg                                                                                                                         <  2:134366/140‑1 (MQ=255)
      aGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTgcg                                                                                                                       <  1:266891/140‑1 (MQ=255)
      aGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTgcg                                                                                                                       <  1:516210/140‑1 (MQ=255)
                  ttAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTc                                                                                                           >  1:931896/1‑140 (MQ=255)
                          aGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGAcc                                                                                                   >  2:79661/1‑140 (MQ=255)
                                     aaaTCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTgg                                                                                        >  2:419796/1‑140 (MQ=255)
                                      aaTCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGt                                                                                       >  2:688897/1‑140 (MQ=255)
                                        tCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCg                                                                                     <  2:183396/140‑1 (MQ=255)
                                            tAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGTATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTgg                                                                                 >  2:704594/1‑140 (MQ=255)
                                            tAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTgg                                                                                 >  2:389767/1‑140 (MQ=255)
                                             aaCCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGt                                                                                <  1:124580/140‑1 (MQ=255)
                                              aCCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTg                                                                                  >  1:200554/1‑137 (MQ=255)
                                              aCCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTg                                                                                  <  2:200554/137‑1 (MQ=255)
                                              aCCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTc                                                                               <  1:688897/140‑1 (MQ=255)
                                                ccGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAg                                                                             <  1:596217/140‑1 (MQ=255)
                                                  gATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCa                                                                           <  1:537816/140‑1 (MQ=255)
                                                  gATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCa                                                                           <  1:533430/140‑1 (MQ=255)
                                                   aTACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCAc                                                                          >  1:93124/1‑140 (MQ=255)
                                                    tACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCAcc                                                                         <  1:340550/140‑1 (MQ=255)
                                                              cTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGt                                                                                <  1:301887/123‑1 (MQ=255)
                                                              cTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGt                                                                                >  2:301887/1‑123 (MQ=255)
                                                              cTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGAcc                                                               >  1:872173/1‑140 (MQ=255)
                                                                tCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACcgc                                                             <  1:533205/140‑1 (MQ=255)
                                                                    cccGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCgtgt                                                         <  2:625372/140‑1 (MQ=255)
                                                                        gACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGAt                                                     >  2:462820/1‑140 (MQ=255)
                                                                          cccGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGc                                                   <  1:419796/140‑1 (MQ=255)
                                                                           ccGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCt                                                  <  1:840650/140‑1 (MQ=255)
                                                                             gTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCCCCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCAGTGCCGTGACCTGAAAGACTGGTCGCTGGTCAGCACCCCGTTGGACCGTGTGTCGATGCTgg                                                <  2:701930/140‑1 (MQ=255)
                                                                                 cTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACAt                                            <  2:93124/140‑1 (MQ=255)
                                                                                     gccgccAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACAt                                            <  1:325593/136‑1 (MQ=255)
                                                                                     gccgccAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACAt                                            >  2:325593/1‑136 (MQ=255)
                                                                                      ccgccAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAgg                                       <  2:931896/140‑1 (MQ=255)
                                                                                        gccAGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGc                                     <  1:737805/140‑1 (MQ=255)
                                                                                           aGGGCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAAc                                  <  1:462820/140‑1 (MQ=255)
                                                                                            gggCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAAcc                                 >  2:224426/1‑140 (MQ=255)
                                                                                            gggCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAAcc                                 <  1:224426/140‑1 (MQ=255)
                                                                                              gCGAGGACTACTACATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCg                               <  1:167696/140‑1 (MQ=255)
                                                                                                        ctacATCGCCACCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCcggcg                     >  2:252092/1‑140 (MQ=255)
                                                                                                                  aCCTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGgc           >  2:770066/1‑140 (MQ=255)
                                                                                                                    cTCGACCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGgcgc         <  1:591917/140‑1 (MQ=255)
                                                                                                                        aCCTTCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTg     >  2:924508/1‑140 (MQ=255)
                                                                                                                           ttCTAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCt  <  2:872173/140‑1 (MQ=255)
                                                                                                                              |                                                                                                                                        
CGCGACAGCCGCATCCCATTAACCGGAGGCAATATGGAAATCACTAACCCGATACTCACCGGCTTCAACCCGGACCCGTCCCTGTGCCGCCAGGGCGAGGACTACTACATCGCCACCTCGACCTTCGAGTGGTTCCCGGGCGTGCGCATCTACCACTCCCGTGACCTGAAAAACTGGTCGCTGGTCAGCACCCCGTTGGACCGCGTGTCGATGCTGGACATGAAGGGCAACCCGGACTCCGGCGGCATCTGGGCGCCGTGCCT  >  NC_000913/286756‑287018

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