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%	157.6 / NA	49	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 (26/23);  total (26/23)

GATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCC  >  NC_000913/1205115‑1205380                                                                                                                                    |                                                                                                                                       gATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTc                                                                                                                                 >  1:1439384/1‑139 (MQ=255)     aTCGCTACATCCTCGAAGAGGCCATAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGaa                                                                                                                             <  2:800343/139‑1 (MQ=255)           aCATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTacac                                                                                                                       >  1:1481730/1‑139 (MQ=255)            cATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTaca                                                                                                                        <  2:521504/137‑1 (MQ=255)            cATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTaca                                                                                                                        >  1:521504/1‑137 (MQ=255)                   gAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGt                                                                                                               <  2:913494/139‑1 (MQ=255)                         gCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCGCATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAAc                                                                                                                   >  1:924567/1‑129 (MQ=255)                         gCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCaa                                                                                                                    <  2:1340572/128‑1 (MQ=255)                         gCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCaa                                                                                                                    >  1:1340572/1‑128 (MQ=255)                         gCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAAc                                                                                                                   <  2:924567/129‑1 (MQ=255)                         gCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAAc                                                                                                                   >  1:820671/1‑129 (MQ=255)                         gCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAAc                                                                                                                   <  2:820671/129‑1 (MQ=255)                         gCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGtccgatc                                                                                                         >  1:755273/1‑139 (MQ=255)                              aGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTc                                                                                                                                 <  2:349717/110‑1 (MQ=255)                              aGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTc                                                                                                                                 >  1:349717/1‑110 (MQ=255)                                 ggcgAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTATCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGaaa                                                                                                 <  1:801516/139‑1 (MQ=255)                                             cAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAgcg                                                                                     <  2:1439384/139‑1 (MQ=255)                                                               aCCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTGCACCAACATGTCCGATCCGATGaaa                                                                                                 >  2:4514/1‑109 (MQ=255)                                                               aCCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTGCACCAACATGTCCGATCCGATGaaa                                                                                                 <  1:4514/109‑1 (MQ=255)                                                                      gACCCCTTCTGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTCAATCGGGGCGCTTtcatc                                                            >  1:1548472/1‑139 (MQ=255)                                                                         cccTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACGCCAACATGt                                                                                                               >  2:979906/1‑85 (MQ=255)                                                                         cccTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGt                                                                                                               <  1:979906/85‑1 (MQ=255)                                                                                     aCCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATcaca                                             >  2:986488/1‑139 (MQ=255)                                                                                      ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCcat                                            >  1:1040175/1‑139 (MQ=255)                                                                                      ccGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCcat                                            >  2:1161222/1‑139 (MQ=255)                                                                                       cGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCcatc                                           <  1:831291/139‑1 (MQ=255)                                                                                       cGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCcatc                                           <  2:1481730/139‑1 (MQ=255)                                                                                          gCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCcatca                                          <  1:848781/137‑1 (MQ=255)                                                                                          gCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCcatca                                          >  2:848781/1‑137 (MQ=255)                                                                                               cacaTGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTg                                   <  2:870712/139‑1 (MQ=255)                                                                                                acaTGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTgg                                  >  2:531335/1‑139 (MQ=255)                                                                                                   tGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTa                                                                                                                          >  2:1058431/1‑48 (MQ=255)                                                                                                   tGACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTa                                                                                                                          <  1:1058431/48‑1 (MQ=255)                                                                                                    gACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCgggg                                                                      <  2:318388/99‑1 (MQ=255)                                                                                                    gACCTTGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCgggg                                                                      >  1:318388/1‑99 (MQ=255)                                                                                                        ttGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCgggg                                                                      >  2:1286007/1‑95 (MQ=255)                                                                                                        ttGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCgggg                                                                      <  1:1286007/95‑1 (MQ=255)                                                                                                         tGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGtccgat                                                                                                          <  2:1523838/58‑1 (MQ=255)                                                                                                         tGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGtccgat                                                                                                          >  1:1523838/1‑58 (MQ=255)                                                                                                         tGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGaa                                                                                                  >  1:1047787/1‑66 (MQ=255)                                                                                                         tGCTGATGAAGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGAACCGATGaa                                                                                                  <  2:1047787/66‑1 (MQ=255)                                                                                                                 aaGACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGt                 <  1:1161222/139‑1 (MQ=255)                                                                                                                   gACCTGAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCg               >  2:1183665/1‑139 (MQ=255)                                                                                                                        gAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTtc                                                               >  2:72187/1‑86 (MQ=255)                                                                                                                        gAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTtc                                                               <  1:72187/86‑1 (MQ=255)                                                                                                                        gAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCaaa          >  2:1542293/1‑139 (MQ=255)                                                                                                                        gAACCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCaaa          >  2:479863/1‑139 (MQ=255)                                                                                                                         aaCCCCGTCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCaaaa         <  2:1612128/139‑1 (MQ=255)                                                                                                                                tCACTATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTcc  >  1:262977/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                       GATTATCGCTACATCCTCGAAGAGGCCAAAGCGGCGAACAAAATCAGCCCGGTCAGTGAGTCACCCATCGACCCCTTCGGGGCGACCGGGCTGTCACATGACCTTGCTGATGAAGACCTGAACCCCGTCACCATCATTCAGAACTACACCAACATGTCCGATCCGATGAAAGAGCTGGAAGCGGCGATTGAATCGGGGCGCTTTCATCATGACGGCAATCCCATCATGACCTGGTGTATCGGCAACGTGGTCGGCAAAACCATTCC  >  NC_000913/1205115‑1205380

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