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_000913	2,259,134	0	T	A	100.0%	91.1 / NA	26	N55Y (AAC→TAC)	psuK	putative pseudouridine kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (13/13);  total (13/13)

CCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTTACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAATTTTATTTTACCTGGATTTGA  >  NC_000913/2259058‑2259209                                                                             |                                                                            ccAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAAcc                                                                          >  1:355145/1‑80 (MQ=255)  cAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAAcc                                                                          >  1:355146/1‑79 (MQ=255)        gattgGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACg                                                                  >  1:355143/1‑81 (MQ=255)         attgGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACg                                                                  <  1:355118/80‑1 (MQ=255)         attgGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACg                                                                  <  2:355119/80‑1 (MQ=255)             gTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAgg                                                             >  2:355142/1‑81 (MQ=255)               ttGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGtt                                                           >  2:355151/1‑81 (MQ=255)                tGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGtt                                                           <  1:355120/80‑1 (MQ=255)                 gCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGttt                                                          <  2:355121/80‑1 (MQ=255)                           cGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATAtt                                               <  2:355122/81‑1 (MQ=255)                            gATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTg                                              <  2:355123/81‑1 (MQ=255)                                 aCCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGccc                                         >  1:355144/1‑81 (MQ=255)                                             aCTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAgg                             <  1:355124/81‑1 (MQ=255)                                             aCTGCCTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAg                              <  2:355125/80‑1 (MQ=255)                                                 ccTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGTTTTTGTGCAATCTTGCTCCCTACTCCACCAGGCGt                          >  1:355154/1‑80 (MQ=255)                                                  cTACGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTaa                        <  1:355126/81‑1 (MQ=255)                                                    aCGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAat                      >  1:355152/1‑81 (MQ=255)                                                    aCGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAat                      >  1:355156/1‑81 (MQ=255)                                                     cGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTaaa                       <  2:355128/79‑1 (MQ=255)                                                     cGGCGCTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAat                      <  2:355127/80‑1 (MQ=255)                                                         gcTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAattttat                 >  1:355150/1‑81 (MQ=255)                                                         gcTCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAattttat                 >  2:355153/1‑81 (MQ=255)                                                           tCAGTAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAattttattt               >  2:355158/1‑81 (MQ=255)                                                               tAGCCAGGCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAATTTTATTTTAcc           >  2:355165/1‑81 (MQ=255)                                                                      gCTTTGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAATTTTATTTTACCTGGAttt    <  1:355129/81‑1 (MQ=255)                                                                          tGTAACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAATTTTATTTTACCTGGATTTGa  <  2:355130/79‑1 (MQ=255)                                                                             |                                                                            CCAGATTGATTGGTTTGCGTTAGCAGCGATTGACCATAAAAATCACTGCCTACGGCGCTCAGTAGCCAGGCTTTGTTACCCAGCAACGCCAGGTTTTGTGCAATATTGCGCCCTACTCCACCAGGCGTAAATTTTATTTTACCTGGATTTGA  >  NC_000913/2259058‑2259209

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