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%	127.3 / NA	40	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 (21/19);  total (21/19)

AATATAATTAAAAATAAAATTTTTGCGT‑AAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATT  >  NC_000913/2995197‑2995448                                                                                                                       |                                                                                                                                       aaTTTAATTAAAAATAAAATTTTTGCGTAAAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGc                                                                                                                    <  1:345138/139‑1 (MQ=255)                                         acaGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAAc                                                                            >  1:1009772/1‑139 (MQ=255)                                                ttAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAAGTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCa                                                                                             <  2:790451/115‑1 (MQ=255)                                                ttAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCa                                                                                             >  1:790451/1‑115 (MQ=255)                                                         cATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAaat                                                            >  1:1419886/1‑139 (MQ=255)                                                                       aaTATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAat                                                                                                         >  1:1134298/1‑80 (MQ=255)                                                                       aaTATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAat                                                                                                         <  2:1134298/80‑1 (MQ=255)                                                                         tatTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTc                                            >  1:222659/1‑139 (MQ=255)                                                                            tACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCt                                                                                                                              >  2:291786/1‑54 (MQ=255)                                                                            tACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCt                                                                                                                              <  1:291786/54‑1 (MQ=255)                                                                            tACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAggg                                                                                                                         <  2:1369463/59‑1 (MQ=255)                                                                            tACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAggg                                                                                                                         >  1:1369463/1‑59 (MQ=255)                                                                                 ttAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGAttt                                    >  1:639892/1‑139 (MQ=255)                                                                                     aTTATATTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGtt                                >  2:651902/1‑139 (MQ=255)                                                                                         tatTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTgg                                                   >  1:320896/1‑116 (MQ=255)                                                                                         tatTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTgg                                                   <  2:320896/116‑1 (MQ=255)                                                                                         tatTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAcc                            >  1:550553/1‑139 (MQ=255)                                                                                         tatTCAAATCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAAcc                            >  1:1316049/1‑139 (MQ=255)                                                                                                aTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAg                     <  2:550553/139‑1 (MQ=255)                                                                                                aTCATTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAg                     <  2:222659/139‑1 (MQ=255)                                                                                                   aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGaa                                                                                                          <  2:682722/51‑1 (MQ=255)                                                                                                   aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGaa                                                                                                          >  1:682722/1‑51 (MQ=255)                                                                                                   aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTaaa                                                             >  2:941739/1‑96 (MQ=255)                                                                                                   aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTaaa                                                             <  1:941739/96‑1 (MQ=255)                                                                                                   aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTc                       <  2:652214/134‑1 (MQ=255)                                                                                                   aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTc                       >  1:652214/1‑134 (MQ=255)                                                                                                   aTTAAACTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAgggg                  <  2:1316049/139‑1 (MQ=255)                                                                                                         cTTGAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAAc            <  2:1009772/139‑1 (MQ=255)                                                                                                            gAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACa                                                                                                             >  1:1241563/1‑39 (MQ=255)                                                                                                            gAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACa                                                                                                             <  2:1241563/39‑1 (MQ=255)                                                                                                            gAGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAAt         >  1:1590282/1‑139 (MQ=255)                                                                                                             aGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAgggg                  <  2:579795/129‑1 (MQ=255)                                                                                                             aGCCAGGGAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAgggg                  >  1:579795/1‑129 (MQ=255)                                                                                                                  gggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAAt         <  1:835352/133‑1 (MQ=255)                                                                                                                  gggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAAt         >  2:835352/1‑133 (MQ=255)                                                                                                                  gggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGattat   <  2:1419886/139‑1 (MQ=255)                                                                                                                   ggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCata                                                                                                 <  2:1179483/44‑1 (MQ=255)                                                                                                                   ggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGtt                                >  2:472895/1‑109 (MQ=255)                                                                                                                   ggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGtt                                <  1:472895/109‑1 (MQ=255)                                                                                                                   ggAGTTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGattatt  >  2:439048/1‑139 (MQ=255)                                                                                                                       |                                                                                                                                       AATATAATTAAAAATAAAATTTTTGCGT‑AAAAAAATACCACAGGCATTAAAAAATCATGAGATGATTAAAATATTACAATTAGATTATATTCAAATCATTAAACTTGAGCCAGGGAGCTATAAGTCCTCAGGGCTTGCCAAAACAGAATTATTCATATCCAAAGAATGTGTTTCAAACATTTTTTCAACTAAATAATGTCTGGCTTGTTCGCAGATTTCGTTAACCTTGTCAGGGGAGAAACAATGATTATT  >  NC_000913/2995197‑2995448

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