Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	586,073	C→T	pseudogene (17/45 nt)	pauD →	tRNA‑OTHER

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	586,073	0	C	T	100.0%	178.4 / NA	46	pseudogene (17/45 nt)	pauD	tRNA‑OTHER
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (22/24);  total (22/24)

AAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCCGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCT  >  NC_000913/585926‑586211                                                                                                                                                    |                                                                                                                                           aaaaGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGc                                                                                                                                          >  1:76123/1‑150 (MQ=255)     gagGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGtt                                                                                                                                     >  1:198569/1‑151 (MQ=255)     gagGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGtt                                                                                                                                     >  1:231174/1‑151 (MQ=255)      agGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTc                                                                                                                                    <  2:67171/151‑1 (MQ=255)            aaaaCTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATcc                                                                                                                              >  1:423848/1‑151 (MQ=255)              aaCTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTg                                                                                                                            >  2:377897/1‑151 (MQ=255)                cTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGc                                                                                                                           >  2:478658/1‑150 (MQ=255)                cTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGc                                                                                                                           >  1:495193/1‑150 (MQ=255)                    gaCACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGgc                                                                                                                      >  2:154235/1‑151 (MQ=255)                    gaCACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGgc                                                                                                                      <  2:304318/151‑1 (MQ=255)                          tAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGcc                                                                                                                 >  1:566230/1‑150 (MQ=255)                            aGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCAtt                                                                                                              <  1:368316/151‑1 (MQ=255)                             ggCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCAttt                                                                                                             <  2:479213/151‑1 (MQ=255)                                cTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATttctt                                                                                                          >  2:515055/1‑151 (MQ=255)                                cTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATttctt                                                                                                          >  1:49096/1‑151 (MQ=255)                                   gCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCt                                                                                                       >  2:550815/1‑151 (MQ=255)                                           ttGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATg                                                                                               >  1:613810/1‑151 (MQ=255)                                            tGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTAt                                                                                                <  2:570535/149‑1 (MQ=255)                                             gCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGcc                                                                                             <  1:262155/151‑1 (MQ=255)                                             gCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGc                                                                                              <  2:38396/150‑1 (MQ=255)                                             gCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGc                                                                                              <  2:317162/150‑1 (MQ=255)                                                aGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGt                                                                                          >  2:579071/1‑151 (MQ=255)                                                       aCATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATcc                                                                                   >  1:38011/1‑151 (MQ=255)                                                         aTGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGt                                                                                 <  1:325512/151‑1 (MQ=255)                                                                gAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCt                                                                          <  1:13083/151‑1 (MQ=255)                                                                 aGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCtt                                                                         >  1:310510/1‑151 (MQ=255)                                                                        gggACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCtt                                                                  >  1:182960/1‑151 (MQ=255)                                                                          gACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTca                                                                <  2:42665/151‑1 (MQ=255)                                                                                  aGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCa                                                        <  1:484855/151‑1 (MQ=255)                                                                                     gcAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACtt                                                     <  1:599349/151‑1 (MQ=255)                                                                                      cAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACttt                                                    <  1:558677/151‑1 (MQ=255)                                                                                         agaTGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTg                                                 <  2:627675/151‑1 (MQ=255)                                                                                                    tAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTgg                                      >  1:371680/1‑151 (MQ=255)                                                                                                        ttgttgACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCg                                    >  1:416790/1‑149 (MQ=255)                                                                                                                ttAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTc                           >  2:622850/1‑150 (MQ=255)                                                                                                                     aTGTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGgt                     >  1:81018/1‑151 (MQ=255)                                                                                                                       gTCGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGgtgt                   <  1:310060/151‑1 (MQ=255)                                                                                                                         cGTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTc                  <  2:190904/150‑1 (MQ=255)                                                                                                                          gTTCTAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAg                <  1:629214/151‑1 (MQ=255)                                                                                                                              tAGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCGGTCTTCACCGTTTCCGCTTTGCCCTTCACCACACCACTTTTTGCCGCGGTTTGGCGAGATAATTCAACGGTGCCAGGc              >  2:315902/1‑149 (MQ=38)                                                                                                                               aGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGcc             <  1:550826/149‑1 (MQ=255)                                                                                                                               aGGAACTTCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGc              <  1:545182/148‑1 (MQ=255)                                                                                                                                      tCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAAc      <  2:263325/149‑1 (MQ=255)                                                                                                                                      tCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACg     <  2:521049/150‑1 (MQ=255)                                                                                                                                      tCTAAGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACg     <  1:252785/150‑1 (MQ=255)                                                                                                                                         aaGTCGTGGGCTGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCt  <  2:611403/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                           AAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGGAACTTCTAAGTCGTGGGCCGCAGGTTCGAATCCTGCAGGGCGCGCCATTTCTTCCTCATTTATGCCCGTCTTATCCGTTTCCGCTTTGCCCTTCACCACATCACTTTTTGTCGCTGTTTGGCGAGATAATTCAACGGTGTCAGGCCGTAAAACGCCT  >  NC_000913/585926‑586211

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