Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 287,971 A→G T395A (ACG→GCG)  yagH → CP4‑6 prophage, putative xylosidase/arabinosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913287,9710AG100.0% 114.0 / NA 32T395A (ACG→GCG) yagHCP4‑6 prophage, putative xylosidase/arabinosidase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/16);  total (16/16)

TCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCATCAAAG  >  NC_000913/287835‑288094
                                                                                                                                        |                                                                                                                           
tCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGcg                                                                                                                           <  1:146602/139‑1 (MQ=255)
            cacCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTtc                                                                                                               >  2:417429/1‑139 (MQ=255)
                      tcgctcACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGc                                                                                                     <  2:348709/139‑1 (MQ=255)
                                  cgcCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCaga                                                                                         >  1:487619/1‑139 (MQ=255)
                                     ccGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                                      <  1:483571/139‑1 (MQ=255)
                                     ccGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                                      <  1:417429/139‑1 (MQ=255)
                                     ccGGGCTTCTTACGGCTCTACGGCAACGACTCGCTCAATTCGACCTTCACCCCAACGACCGTGGCGCGCCCCTGGCAGGACTTCGCCATCCGCGCAGAAGCCCGGATGGATTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                                      <  1:541489/139‑1 (MQ=255)
                                        ggCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCgg                                                                                   >  2:130340/1‑139 (MQ=255)
                                              ttACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGa                                                                             <  2:503253/139‑1 (MQ=255)
                                                   gCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGt                                                                                                       <  1:426090/108‑1 (MQ=255)
                                                   gCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGt                                                                                                       >  2:426090/1‑108 (MQ=255)
                                                       tATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGAAGCACATCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGctact                                                                    >  1:258106/1‑139 (MQ=255)
                                                                               aCCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                                      <  1:143770/97‑1 (MQ=255)
                                                                               aCCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAgcg                                                                                      >  2:143770/1‑97 (MQ=255)
                                                                                ccTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTg                                           >  2:296557/1‑139 (MQ=255)
                                                                                     aCCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTCCTACAACAGCAAAAACTGGCGCTACTGCTTtg                                      >  1:34522/1‑139 (MQ=255)
                                                                                       ccAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTtgtg                                    <  2:487619/139‑1 (MQ=255)
                                                                                              aCCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGc                                                                        <  1:196971/96‑1 (MQ=255)
                                                                                              aCCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGc                                                                        >  2:196971/1‑96 (MQ=255)
                                                                                              aCCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACg                             >  1:430866/1‑139 (MQ=255)
                                                                                               ccGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGAc                                                                            <  1:214968/91‑1 (MQ=255)
                                                                                               ccGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGAc                                                                            >  2:214968/1‑91 (MQ=255)
                                                                                                   ggCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACgaggag                        >  1:261591/1‑139 (MQ=255)
                                                                                                     cgcgcCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTtg                                      <  1:140643/123‑1 (MQ=255)
                                                                                                     cgcgcCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTtg                                      <  1:173865/123‑1 (MQ=255)
                                                                                                     cgcgcCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTtg                                      >  2:173865/1‑123 (MQ=255)
                                                                                                     cgcgcCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTtg                                      >  2:140643/1‑123 (MQ=255)
                                                                                                      gcgcCGCTGGCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGa                     <  2:266779/139‑1 (MQ=255)
                                                                                                              ggCAGCACTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAg             <  1:296557/139‑1 (MQ=255)
                                                                                                               gcagcaCTTCGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGa            <  2:430866/139‑1 (MQ=255)
                                                                                                                       acGCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCATCaa    >  1:381624/2‑139 (MQ=255)
                                                                                                                         gCCTTCCGGGCAGAAGCGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCATCAAAg  >  1:568600/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                           
TCCGTTCGACGACACCCTCGGCTCGCTCACCGCGCGCCCGGGCTTCTTACGGCTCTATGGCAACGACTCGCTCAATTCGACCTTCACCCAATCGACCGTGGCGCGCCGCTGGCAGCACTTCGCCTTCCGGGCAGAAACGCGGATGGAGTTCTCGCCGGTGCACTTCCAGCAGAGCGCGGGGCTGACCTGCTACTACAACAGCAAAAACTGGAGCTACTGCTTTGTGGACTACGAGGAGGGACAGGGTAGAACCATCAAAG  >  NC_000913/287835‑288094

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