Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,381,789	G→A	pseudogene (13/126 nt)	ycjV →	pseudogene, putative ATP‑binding component of a transport system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,381,789	0	G	A	100.0%	126.4 / NA	40	pseudogene (13/126 nt)	ycjV	pseudogene, putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/21);  total (19/21)

ATCGTGATGGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCC  >  NC_000913/1381655‑1381924                                                                                                                                       |                                                                                                                                        aTCGTGATGGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                                     <  2:353835/139‑1 (MQ=255)  tCGTGATGGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTc                                                                                                                                    >  2:93746/1‑139 (MQ=255)                                 tCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATc                                                                                                           >  1:127104/1‑133 (MQ=255)                                 tCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATc                                                                                                           <  2:127104/133‑1 (MQ=255)                                   cGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGga                                                                                                   <  1:1038166/139‑1 (MQ=255)                                     gACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGgaga                                                                                                 >  2:530414/1‑139 (MQ=255)                                       cgcgCAAGA‑GAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTCCGTTAAATGATTATCATGCAGGAGaaaa                                                                                              <  1:688786/139‑1 (MQ=255)                                          gcAAGAGGAAAATAACATTTCCGCCAAAGTTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                                     >  1:689511/1‑98 (MQ=255)                                          gcAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                                     <  2:689511/98‑1 (MQ=255)                                          gcAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAata                                                                                            >  1:363373/1‑139 (MQ=255)                                             agagGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCa                                                                                         <  2:812068/139‑1 (MQ=255)                                                        aCATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCAtt                                                                                >  1:261491/1‑137 (MQ=255)                                                        aCATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCAtt                                                                                <  2:261491/137‑1 (MQ=255)                                                            ttCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAatat                                                                                           <  1:758912/122‑1 (MQ=255)                                                            ttCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGata                                                                          >  1:717398/1‑139 (MQ=255)                                                            ttCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGACACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAatat                                                                                           >  2:758912/1‑122 (MQ=255)                                                                           cGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCAtt                                                           >  1:992297/1‑139 (MQ=255)                                                                                        aCCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGa                                              >  1:717397/1‑139 (MQ=255)                                                                                         ccGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGaa                                             <  2:577670/139‑1 (MQ=255)                                                                                          cGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGaaa                                            <  2:1039371/139‑1 (MQ=255)                                                                                             tGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACgg                                         >  2:621798/1‑139 (MQ=255)                                                                                                       aTGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAAt                               <  2:992297/139‑1 (MQ=255)                                                                                                             tACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTa                         <  2:717397/139‑1 (MQ=255)                                                                                                              acacCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAatat                                                                                           <  1:29827/72‑1 (MQ=255)                                                                                                              acacCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAatat                                                                                           >  2:29827/1‑72 (MQ=255)                                                                                                                     ggTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTAt                                                                                     <  1:151632/71‑1 (MQ=255)                                                                                                                     ggTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTAt                                                                                     >  2:151632/1‑71 (MQ=255)                                                                                                                     ggTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAAt                               <  1:438435/125‑1 (MQ=255)                                                                                                                     ggTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAAt                               >  2:438435/1‑125 (MQ=255)                                                                                                                     ggTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAgg                 >  2:39744/1‑139 (MQ=255)                                                                                                                       ttGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAgggg               >  1:95457/1‑139 (MQ=255)                                                                                                                       ttGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAgggg               >  1:679660/1‑139 (MQ=255)                                                                                                                               aCGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGc                                                                                                       >  2:825167/1‑43 (MQ=255)                                                                                                                               aCGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGc                                                                                                       <  1:825167/43‑1 (MQ=255)                                                                                                                               aCGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCAt       >  1:1056804/1‑139 (MQ=255)                                                                                                                               aCGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCAt       <  2:717398/139‑1 (MQ=255)                                                                                                                               aCGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCAt       <  1:516700/139‑1 (MQ=255)                                                                                                                               aCGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCAt       <  1:39744/139‑1 (MQ=255)                                                                                                                                    ttAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTcccc  <  2:1056804/139‑1 (MQ=255)                                                                                                                                    ttAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTcccc  <  2:851215/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                        ATCGTGATGGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCC  >  NC_000913/1381655‑1381924

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