Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,995,314	C→T	pseudogene (447/447 nt)	ygeN ←	pseudogene, orgB family, part of T3SS PAI ETT2 remnant

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,995,314	0	C	T	100.0%	76.7 / NA	26	pseudogene (447/447 nt)	ygeN	pseudogene, orgB family, part of T3SS PAI ETT2 remnant
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/13);  total (13/13)

AAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATT  >  NC_000913/2995245‑2995445                                                                      |                                                                                                                                    aaaaaaTCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCt                                                                                                                           <  2:609825/80‑1 (MQ=255) aaaaaaTCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCt                                                                                                                           >  1:609825/1‑80 (MQ=255)       tCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTaaa                                                          <  2:670248/139‑1 (MQ=255)           gagaTGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCaaaa                                                                                                            <  2:905570/85‑1 (MQ=255)           gagaTGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCaaaa                                                                                                            >  1:905570/1‑85 (MQ=255)                           tACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCa                                      >  1:40710/1‑139 (MQ=255)                            aCAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAg                                     >  2:392537/1‑139 (MQ=255)                                  aGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCg                               <  2:40710/139‑1 (MQ=255)                                        tatTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGaa                                                                                                       >  1:383731/1‑61 (MQ=255)                                        tatTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGaa                                                                                                       <  2:383731/61‑1 (MQ=255)                                        tatTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttttt                                                                  >  1:819735/1‑98 (MQ=255)                                        tatTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttttt                                                                  <  2:819735/98‑1 (MQ=255)                                             aaaTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCa                                                                >  2:40447/1‑95 (MQ=255)                                             aaaTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCa                                                                <  1:40447/95‑1 (MQ=255)                                              aaTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCaa                                                               >  1:572525/1‑95 (MQ=255)                                              aaTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCaa                                                               <  2:572525/95‑1 (MQ=255)                                               aTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttttt                                                                  <  1:665941/91‑1 (MQ=255)                                               aTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACAtttttt                                                                  >  2:665941/1‑91 (MQ=255)                                                  aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCg                               >  1:872852/1‑123 (MQ=255)                                                  aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCg                               <  2:872852/123‑1 (MQ=255)                                                  aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAc                          <  1:67870/128‑1 (MQ=255)                                                  aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAc                          >  2:67870/1‑128 (MQ=255)                                                   ttAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGga              >  1:938420/1‑139 (MQ=255)                                                    tAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAc                          >  2:26716/1‑126 (MQ=255)                                                    tAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAc                          <  1:26716/126‑1 (MQ=255)                                                               cAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGatt  <  1:392537/139‑1 (MQ=255)                                                                      |                                                                                                                                    AAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATT  >  NC_000913/2995245‑2995445

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: ‑