Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,205,246	C→T	pseudogene (843/1134 nt)	aaaE →	pseudogene, phage terminase protein A family, e14 prophage,Phage or Prophage Related

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,205,246	0	C	T	100.0%	106.9 / NA	34	pseudogene (843/1134 nt)	aaaE	pseudogene, phage terminase protein A family, e14 prophage,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/16);  total (18/16)

GGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATT  >  NC_000913/1205114‑1205378                                                                                                                                     |                                                                                                                                     ggATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCAtt                                                                                                                                <  2:188857/139‑1 (MQ=255) ggATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCAtt                                                                                                                                <  1:636537/139‑1 (MQ=255)  gATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTc                                                                                                                               >  2:721120/1‑139 (MQ=255)      aTCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGaa                                                                                                                           <  2:4995/139‑1 (MQ=255)       tCGCTACATCCTCGAAGAGGCCCAAGCGGCGAACAAAATCAGCCCGGCCAGTGAGTCACCCATCGACCCCTCCGGGGCGACCAGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAAc                                                                                                                          >  2:80099/1‑139 (MQ=255)          cTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTac                                                                                                                       >  2:398537/1‑139 (MQ=255)           tACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTaca                                                                                                                      <  1:721120/139‑1 (MQ=255)                  tCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACAt                                                                                                               >  1:568501/1‑139 (MQ=255)                          gCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGaa                                                                                                                           >  1:111820/1‑119 (MQ=255)                          gCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGaa                                                                                                                           <  2:111820/119‑1 (MQ=255)                                  ggcgAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGaaa                                                                                               <  1:370847/139‑1 (MQ=255)                                             tCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAgc                                                                                    <  2:568501/139‑1 (MQ=255)                                              cAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAgcg                                                                                   <  1:469607/139‑1 (MQ=255)                                                     gTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTg                                                                            >  1:735661/1‑139 (MQ=255)                                                              tCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGgc                                                                   <  2:735661/139‑1 (MQ=255)                                                                aCCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGgcgc                                                                 >  2:735881/1‑139 (MQ=255)                                                                      cGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTtcat                                                           >  2:459658/1‑139 (MQ=255)                                                                                   gCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATc                                              <  1:346262/139‑1 (MQ=255)                                                                                       ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGtccga                                                                                                         <  1:46790/76‑1 (MQ=255)                                                                                       ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGtccga                                                                                                         >  2:46790/1‑76 (MQ=255)                                                                                           gCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATg                                      <  1:81438/139‑1 (MQ=255)                                                                                                cacaTGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGCCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTg                                 <  1:425220/139‑1 (MQ=255)                                                                                                         ttGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCggg                                                                     <  2:462954/94‑1 (MQ=255)                                                                                                         ttGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCggg                                                                     >  1:462954/1‑94 (MQ=255)                                                                                                         ttGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCgg                        >  1:259614/1‑139 (MQ=255)                                                                                                         ttGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCgg                        >  1:197433/1‑139 (MQ=255)                                                                                                                 gaAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGTCAACGTgg                >  2:749122/1‑139 (MQ=255)                                                                                                                         gAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCaaa        >  2:528310/1‑139 (MQ=255)                                                                                                                         gAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCaaa        >  2:7789/1‑139 (MQ=255)                                                                                                                               cGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGAcc                                   >  2:411603/1‑106 (MQ=255)                                                                                                                               cGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGAcc                                   <  1:411603/106‑1 (MQ=255)                                                                                                                               cGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCAtt  >  2:238233/1‑139 (MQ=255)                                                                                                                               cGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCAtt  <  1:238233/139‑1 (MQ=255)                                                                                                                               cGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCAtt  >  2:534938/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                     GGATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATT  >  NC_000913/1205114‑1205378

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