Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,130,088	T→G	pseudogene (56/1101 nt)	yghO ←	pseudogene, IS‑interrupted

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,130,088	0	T	G	100.0%	149.0 / NA	42	pseudogene (56/1101 nt)	yghO	pseudogene, IS‑interrupted
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (21/21);  total (21/21)

GGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTC  >  NC_000913/3129943‑3130234                                                                                                                                                  |                                                                                                                                                   ggTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAgcgtcg                                                                                                                                               <  2:418553/151‑1 (MQ=255) ggTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAgcgtcg                                                                                                                                               >  2:347983/1‑151 (MQ=255)          ttGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGt                                                                                                                                      >  2:636229/1‑151 (MQ=255)               cgGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGtgatg                                                                                                                                  >  2:164299/1‑150 (MQ=255)                ggTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGAc                                                                                                                                >  2:407226/1‑151 (MQ=255)                      aCGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGaaa                                                                                                                           <  2:93404/150‑1 (MQ=255)                            aaCTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGAt                                                                                                                     >  1:692399/1‑150 (MQ=255)                                aTCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAAt                                                                                                                 <  1:150333/150‑1 (MQ=255)                                  cTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGcc                                                                                                              >  1:85162/1‑151 (MQ=255)                                    gCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCtt                                                                                                            <  1:621810/151‑1 (MQ=255)                                     ccccGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGcc                                                                                                              >  1:537907/1‑148 (MQ=255)                                     ccccGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCtt                                                                                                            >  1:353950/1‑150 (MQ=255)                                     ccccGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCtt                                                                                                            >  1:424115/1‑150 (MQ=255)                                     ccccGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCtt                                                                                                            <  2:353996/150‑1 (MQ=255)                                          cTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTc                                                                                                      >  2:475119/1‑151 (MQ=255)                                            ttttttGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCa                                                                                                     <  2:269197/150‑1 (MQ=255)                                                 tGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCAttttt                                                                                                <  1:148446/150‑1 (MQ=255)                                                 tGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCAttttt                                                                                                <  2:733504/150‑1 (MQ=255)                                                 tGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCAttttt                                                                                                <  2:665869/150‑1 (MQ=255)                                                 tGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCAttttt                                                                                                <  1:347195/150‑1 (MQ=255)                                                     acccacGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTaa                                                                                           <  1:347949/151‑1 (MQ=255)                                                     acccacGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTaa                                                                                           <  1:407177/151‑1 (MQ=255)                                                      cacacGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTACGGTCATTTttatta                                                                                            >  1:165579/3‑149 (MQ=255)                                                            ccTGCCAGATGATATGGTCCGTCCCCAGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTaagga                                                                                        >  2:207383/1‑147 (MQ=255)                                                             cTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTg                                                                                   >  1:489165/1‑151 (MQ=255)                                                               gCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCg                                                                                 >  1:628914/1‑151 (MQ=255)                                                                   gatgatATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACa                                                                             <  1:186580/151‑1 (MQ=255)                                                                      gatATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGt                                                                          <  1:737442/151‑1 (MQ=255)                                                                         atGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTc                                                                         >  2:185133/1‑149 (MQ=255)                                                                              ccgtccCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATACGGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGa                                                                    >  1:252768/1‑149 (MQ=255)                                                                                                 aGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATc                                               <  1:311379/151‑1 (MQ=255)                                                                                                      aaTGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTTATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTgg                                            <  2:509241/149‑1 (MQ=255)                                                                                                         gCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCa                                       <  2:368228/151‑1 (MQ=255)                                                                                                                 cGTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAg                                    >  2:441858/1‑146 (MQ=255)                                                                                                                  gTTCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAgg                              <  1:573024/151‑1 (MQ=255)                                                                                                                    tCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAggg                             <  1:508951/150‑1 (MQ=255)                                                                                                                    tCAATGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAggg                             <  1:96893/150‑1 (MQ=255)                                                                                                                      aaTGAATAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTtc                          >  1:762398/1‑151 (MQ=255)                                                                                                                           aTAAAGGGGGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTcatcat                     >  2:494573/1‑151 (MQ=255)                                                                                                                            tAAAGGGTGTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATg                    >  1:297406/1‑151 (MQ=255)                                                                                                                                   ggTATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtcatc             >  2:74872/1‑151 (MQ=255)                                                                                                                                     tATCCAGTTTGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGtcatca            >  2:512846/1‑150 (MQ=255)                                                                                                                                             ttGGAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCt    <  1:71776/150‑1 (MQ=255)                                                                                                                                               ggAGCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTc  <  2:81677/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GGTATCTATTTGCGCGGTAATACGCCCAACTATCTGCCCCGCTTTTTTTGCCACCCACGCCTGCCAGATGATATGGTCCGTCCCCGGATTTTTCGCAGACAAATGCTCGTTGCGTTCAATGAATAAAGGGGGTATCCAGTTTGGATCGTCGGGATAAAGTGATGACGGAAAAGCGATAAATGCCTTAAGGTCATTTTTATTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTC  >  NC_000913/3129943‑3130234

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