Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,651,303	T→G	pseudogene (13/741 nt)	yhiS →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,651,303	0	T	G	100.0%	196.7 / NA	53	pseudogene (13/741 nt)	yhiS	pseudogene
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (28/25);  total (28/25)

AAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAACCCATTAATTCCACATTCTCTGGAGAA  >  NC_000913/3651156‑3651449                                                                                                                                                    |                                                                                                                                                   aaaGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTa                                                                                                                                                 <  2:856176/151‑1 (MQ=255)   aGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTAcc                                                                                                                                               >  2:652171/1‑151 (MQ=255)       aaaTCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAg                                                                                                                                           <  2:18826/151‑1 (MQ=255)        aaTCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAgg                                                                                                                                          <  1:503823/151‑1 (MQ=255)              aGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTataat                                                                                                                                    >  1:114795/1‑151 (MQ=255)               ggTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTataat                                                                                                                                    >  2:525144/1‑150 (MQ=255)                 tAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTAtaataa                                                                                                                                  >  1:730979/1‑150 (MQ=255)                  aaTTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATaataaat                                                                                                                                >  2:567825/1‑151 (MQ=255)                  aaTTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATaataaa                                                                                                                                 <  2:362240/150‑1 (MQ=255)                      atatCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACat                                                                                                                            >  2:14153/1‑151 (MQ=255)                        atCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACatat                                                                                                                          <  1:162176/151‑1 (MQ=255)                          cAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATc                                                                                                                         <  2:712307/150‑1 (MQ=255)                          cAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATc                                                                                                                         >  1:95233/1‑150 (MQ=255)                             gAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACgg                                                                                                                     <  2:858500/151‑1 (MQ=255)                                gATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCg                                                                                                                   >  2:632607/1‑150 (MQ=255)                                  tataACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACAGCGat                                                                                                                 >  1:392569/1‑150 (MQ=255)                                   ataACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGatat                                                                                                               >  2:702846/1‑151 (MQ=255)                                      aCGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTta                                                                                                            <  1:877634/151‑1 (MQ=255)                                        gCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTtat                                                                                                           >  1:497547/1‑150 (MQ=255)                                             aCATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGc                                                                                                     <  1:137604/151‑1 (MQ=255)                                              cATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGca                                                                                                    <  1:827223/151‑1 (MQ=255)                                                  tATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCacaa                                                                                                 <  1:115420/150‑1 (MQ=255)                                                   aTTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCacaac                                                                                                >  1:383395/1‑150 (MQ=255)                                                      gtgtgtAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACaaa                                                                                             <  2:711632/150‑1 (MQ=255)                                                             gaagaaAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATa                                                                                     <  1:857900/151‑1 (MQ=255)                                                                 aaaaTCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGc                                                                                 <  1:104967/151‑1 (MQ=255)                                                                   aaTCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCt                                                                                >  1:543379/1‑150 (MQ=255)                                                                        aTTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACg                                                                          <  1:743631/151‑1 (MQ=255)                                                                                tCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGacac                                                                  >  1:918665/1‑151 (MQ=255)                                                                                 cTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAg                                                                      <  1:626885/146‑1 (MQ=255)                                                                                  tAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGacacac                                                                >  2:338677/1‑151 (MQ=255)                                                                                      tATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGa                                                             <  1:106779/150‑1 (MQ=255)                                                                                       aTCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGa                                                             >  2:862266/1‑149 (MQ=255)                                                                                            aTTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGcc                                                      >  1:783311/1‑151 (MQ=255)                                                                                              ttAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCtt                                                    >  1:711883/1‑151 (MQ=255)                                                                                              ttAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCtt                                                    >  1:547011/1‑151 (MQ=255)                                                                                               tAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCtt                                                    <  2:442318/150‑1 (MQ=255)                                                                                                aaaCATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGc                                                  >  1:50327/1‑151 (MQ=255)                                                                                                   cATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGa                                                >  2:918533/1‑150 (MQ=255)                                                                                                   cATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGAt                                               <  1:114084/151‑1 (MQ=255)                                                                                                   cATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGAt                                               >  1:228553/1‑151 (MQ=255)                                                                                                      cTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAAt                                            >  1:808258/1‑151 (MQ=255)                                                                                                           cATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGga                                       <  2:876641/151‑1 (MQ=255)                                                                                                               aaaCAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGgaagaa                                   <  1:806215/151‑1 (MQ=255)                                                                                                                aaCAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGaac                                  <  1:404631/151‑1 (MQ=255)                                                                                                                 aCAAATGATGAATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGaaca                                 <  1:131509/151‑1 (MQ=255)                                                                                                                       gatgaATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAAcc                           >  2:507290/1‑151 (MQ=255)                                                                                                                          gaATAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAACCCa                         <  1:401101/150‑1 (MQ=255)                                                                                                                             tAAGGAATTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAACCCATTa                      >  2:911837/1‑150 (MQ=255)                                                                                                                                  aaTTTTATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAACCCATTAATTcc                 >  2:254709/1‑150 (MQ=255)                                                                                                                                       tATGTCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAACCCATTAATTCCACATtc           >  2:729397/1‑151 (MQ=255)                                                                                                                                           tCAATTGACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAACCCATTAATTCCACATTCTCTg       >  2:753969/1‑151 (MQ=255)                                                                                                                                                 gACGTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAACCCATTAATTCCACATTCTCTGGAGaa  <  1:786397/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   AAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTTTATGTCAATTGACTTTACCCCAGGTATAATAAATACATATCACGGCGATATTTATAACTGCACAACAAATACCGATAATGCTAAAACGCCAGACACACCGAAGTGGCCTTGCGATAATTGGGAAGAACAACAACCCATTAATTCCACATTCTCTGGAGAA  >  NC_000913/3651156‑3651449

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