Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,259,134 T→A N55Y (AAC→TAC)  psuK ← putative pseudouridine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,259,1340TA100.0% 61.7 / NA 17N55Y (AAC→TAC) psuKputative pseudouridine kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (7/10);  total (7/10)

CCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTTACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAATTTTATTTTACCTGGATTTGAA  >  NC_000913/2259058‑2259210
                                                                            |                                                                            
ccAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAAccc                                                                          <  1:402858/81‑1 (MQ=255)
ccAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAAccc                                                                          <  1:402857/81‑1 (MQ=255)
      tgattgGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAAc                                                                    >  2:402886/1‑81 (MQ=255)
        attgGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGc                                                                  <  1:402859/81‑1 (MQ=255)
        attgGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGc                                                                  <  1:402860/81‑1 (MQ=255)
             tttGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGt                                                             >  2:857057/1‑81 (MQ=255)
                 cGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGtttt                                                          <  2:402861/80‑1 (MQ=255)
                  gTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTtgt                                                        >  2:402885/1‑81 (MQ=255)
                      gcagcGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCaa                                                    >  1:402879/1‑81 (MQ=255)
                            aTTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTgc                                              <  2:402862/81‑1 (MQ=255)
                            aTTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTgc                                              <  2:402863/81‑1 (MQ=255)
                                    tAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACt                                      <  1:402864/81‑1 (MQ=255)
                                       aaaTCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTcca                                   >  1:402902/1‑81 (MQ=255)
                                            aCTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAgg                              >  1:402883/1‑81 (MQ=255)
                                            aCTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAg                               <  2:402865/80‑1 (MQ=255)
                                                                 ccAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAATTTTATTTTACCTgg         <  1:402866/81‑1 (MQ=255)
                                                                        ttGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAATTTTATTTTACCTGGATTTGaa  >  2:402887/1‑81 (MQ=255)
                                                                            |                                                                            
CCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTTACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAATTTTATTTTACCTGGATTTGAA  >  NC_000913/2259058‑2259210

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