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%	131.9 / NA	42	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 (20/22);  total (20/22)

GGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCC  >  NC_000913/1381663‑1381924                                                                                                                               |                                                                                                                                        gggAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGAGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGATAATGGTCCGTGCTg                                                                                                                             >  1:497643/1‑139 (MQ=255)  ggAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTgg                                                                                                                            <  1:590175/139‑1 (MQ=255)                    atTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCa                                                                                                          <  2:258429/139‑1 (MQ=255)                           cGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGga                                                                                                   <  2:146760/139‑1 (MQ=255)                           cGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGga                                                                                                   <  2:22474/139‑1 (MQ=255)                               cgcgCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGaaa                                                                                               >  2:316224/1‑139 (MQ=255)                                            aaTAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                                     <  2:667951/88‑1 (MQ=255)                                            aaTAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGt                                                                                                                                     >  1:667951/1‑88 (MQ=255)                                                                  gCGTGGCAGAATTAACCGGTGCGGAATGTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAAGGATTATCATGCAGGAGAAAATATCACTATTCATATTGATATGACGAAATGTCAt                                                            <  2:497643/139‑1 (MQ=255)                                                                       gCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGt                                                                                                                           <  1:557778/71‑1 (MQ=255)                                                                       gCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGt                                                                                                                           <  1:729190/71‑1 (MQ=255)                                                                       gCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGt                                                                                                                           >  2:557778/1‑71 (MQ=255)                                                                       gCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGt                                                                                                                           >  2:729190/1‑71 (MQ=255)                                                                       gCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATg                                                                                                        >  2:105980/1‑90 (MQ=255)                                                                       gCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATg                                                                                                        <  1:105980/90‑1 (MQ=255)                                                                                      gCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGa                                        >  2:36008/1‑139 (MQ=255)                                                                                             ttATGCTCTACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCa                                 <  2:135369/139‑1 (MQ=255)                                                                                                    ctACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAAt                                                                 >  1:511715/1‑100 (MQ=255)                                                                                                    ctACACCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAAt                                                                 <  2:511715/100‑1 (MQ=255)                                                                                                      acacCACGGTTGGGGGCACGAGTTAATGTTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTc                             <  2:386141/134‑1 (MQ=255)                                                                                                      acacCACGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGAAAGTc                             >  1:386141/1‑134 (MQ=255)                                                                                                       catcacGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTa                                                                                                             <  1:163169/50‑1 (MQ=255)                                                                                                       catcacGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTa                                                                                                             >  2:163169/4‑53 (MQ=255)                                                                                                       caccacGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAatat                                                                                           >  2:231953/1‑71 (MQ=255)                                                                                                       caccacGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAatat                                                                                           <  1:231953/71‑1 (MQ=255)                                                                                                       caccacGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTaaa                       >  1:575741/1‑139 (MQ=255)                                                                                                            cGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTa                                                                                                             >  1:703274/1‑48 (MQ=255)                                                                                                            cGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTa                                                                                                             <  2:703274/48‑1 (MQ=255)                                                                                                            cGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAAt                               <  1:514804/126‑1 (MQ=255)                                                                                                            cGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAAt                               >  2:514804/1‑126 (MQ=255)                                                                                                            cGGTTGGGGGCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAg                  >  2:610778/1‑139 (MQ=255)                                                                                                                 gggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATGATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGa             >  1:712998/1‑139 (MQ=255)                                                                                                                 gggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGga                                                                                                   >  2:305781/1‑53 (MQ=255)                                                                                                                 gggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGga                                                                                                   <  1:305781/53‑1 (MQ=255)                                                                                                                 gggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGa             >  1:854746/1‑139 (MQ=255)                                                                                                                 gggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGa             >  2:287991/1‑139 (MQ=255)                                                                                                                  ggggCACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGaa            >  2:862033/1‑139 (MQ=255)                                                                                                                      cACGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCa        <  1:287991/139‑1 (MQ=255)                                                                                                                       aCGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCAt       <  1:316224/139‑1 (MQ=255)                                                                                                                       aCGAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCAt       <  2:854746/139‑1 (MQ=255)                                                                                                                         gAGTTAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTc     >  1:453188/1‑139 (MQ=255)                                                                                                                            ttAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATCCTGGGGGAAGGCATTcccc  <  2:672561/139‑1 (MQ=255)                                                                                                                            ttAATGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTcccc  <  2:575741/139‑1 (MQ=255)                                                                                                                               |                                                                                                                                        GGGAATACGACCGGAAGATATTCATCCGGACGCGCAAGAGGAAAATAACATTTCCGCCAAAATTAGCGTGGCAGAATTAACCGGTGCGGAATTTATGCTCTACACCACGGTTGGGGGCACGAGTTAGTGGTCCGTGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCC  >  NC_000913/1381663‑1381924

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