Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	263,172	A→C	pseudogene (41/63 nt)	ykfN ←	pseudogene, YdiA family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	263,172	0	A	C	98.5%	246.3 / ‑4.1	67	pseudogene (41/63 nt)	ykfN	pseudogene, YdiA family
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (29/37);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

CTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAC  >  NC_000913/263077‑263267                                                                                                |                                                                                                cTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTc                                                                                             <  1:576640/100‑1 (MQ=255)   tCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCAt                                                                                           >  2:1037985/1‑100 (MQ=255)       tCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATg                                                                                       <  1:32040/100‑1 (MQ=255)         aCTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGAt                                                                                     <  2:1232419/100‑1 (MQ=255)          cTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATc                                                                                    >  1:899429/1‑100 (MQ=255)             taCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTCAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGt                                                                                 <  1:1334057/100‑1 (MQ=255)             taCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGt                                                                                 <  1:1088306/100‑1 (MQ=255)             taCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGt                                                                                 <  2:85565/100‑1 (MQ=255)              aCATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGtt                                                                                <  1:1457597/100‑1 (MQ=255)                aTCCTTGGCTGCAGGTTTAGTTGTACACCAATCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTa                                                                               <  1:348431/99‑1 (MQ=255)                  ccTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATa                                                                            >  1:332554/1‑100 (MQ=255)                  ccTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATa                                                                            <  2:1159185/100‑1 (MQ=255)                  ccTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATa                                                                            <  1:1568181/100‑1 (MQ=255)                   cTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATAc                                                                           <  1:1114009/100‑1 (MQ=255)                   cTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATAc                                                                           >  2:1192291/1‑100 (MQ=255)                    ttGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACt                                                                          <  1:822281/100‑1 (MQ=255)                       gCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCa                                                                       <  1:1401518/100‑1 (MQ=255)                         tGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCAAACAAATTCGCTCATTATGATCGGTTAATACTTcca                                                                      >  1:270484/1‑97 (MQ=255)                             ggTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCt                                                                 >  2:375259/1‑100 (MQ=255)                              gTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTg                                                                <  1:301231/100‑1 (MQ=255)                              gTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTg                                                                <  1:372276/100‑1 (MQ=255)                                ttAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGt                                                              <  1:375256/100‑1 (MQ=255)                                ttAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGt                                                              >  2:1010177/1‑100 (MQ=255)                                ttAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGt                                                              >  2:680699/1‑100 (MQ=255)                                   gTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGc                                                           <  1:266005/100‑1 (MQ=255)                                      gTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATg                                                        >  1:974437/1‑100 (MQ=255)                                      gTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATg                                                        >  1:150110/1‑100 (MQ=255)                                      gTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATg                                                        >  1:1495999/1‑100 (MQ=255)                                      gTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATg                                                        <  2:724396/100‑1 (MQ=255)                                        acacGACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAt                                                      >  1:100391/1‑100 (MQ=255)                                        acacCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAt                                                      <  1:1134389/100‑1 (MQ=255)                                            cacTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtt                                                  >  2:81290/1‑100 (MQ=255)                                            cacTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtt                                                  <  1:395746/100‑1 (MQ=255)                                              cTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTGTGATCTGTTAATACTTCAACTTCTGGTTGCATGAttttttg                                                >  1:121680/1‑100 (MQ=255)                                              cTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGAttgtttg                                                <  1:1476016/100‑1 (MQ=255)                                                ccTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTc                                              <  2:539527/100‑1 (MQ=255)                                                ccTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTc                                              >  1:880928/1‑100 (MQ=255)                                                 cTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTcc                                             >  1:196320/1‑100 (MQ=255)                                                  tAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCg                                            <  1:828932/100‑1 (MQ=255)                                                  tAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCg                                            <  1:1297598/100‑1 (MQ=255)                                                  tAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCg                                            <  2:754083/100‑1 (MQ=255)                                                    aaTTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTa                                          >  2:1368905/1‑100 (MQ=255)                                                    aaTTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTa                                          <  2:992583/100‑1 (MQ=255)                                                     aTTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTaa                                         >  1:471952/1‑100 (MQ=255)                                                      tttAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTaaa                                        <  2:1262716/100‑1 (MQ=255)                                                       ttAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTaaaa                                       <  1:473915/100‑1 (MQ=255)                                                            gtgtTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATa                                  >  2:211474/1‑100 (MQ=255)                                                              gtTGGCAATGTGTTCAATAAAGCTCGCACAAATTCGCTCATTATGATCGGTTAATACTTCCACTTCTGGTTGCATGATTGTTTGTCCGTCAAAAGATAAc                                >  1:499370/1‑100 (MQ=255)                                                              gtTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAAc                                >  2:743574/1‑100 (MQ=255)                                                              gtTGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGATGCATGATTGTTTGTCCGTAAAAACAt                                   >  1:1351927/1‑97 (MQ=255)                                                               ttcgctaTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAAc                                <  1:1461100/93‑1 (MQ=255)                                                                tGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAAcgc                              >  1:451197/1‑100 (MQ=255)                                                                tGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAAcgc                              <  1:1345906/100‑1 (MQ=255)                                                                tGGCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAAcgc                              >  2:138365/1‑100 (MQ=255)                                                                 ggCAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAAcgcg                             <  2:986535/100‑1 (MQ=255)                                                                  ctaaTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACgcgc                            <  1:696797/98‑1 (MQ=255)                                                                   cAATGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGcc                           <  2:899436/100‑1 (MQ=255)                                                                    aaTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACgcgcag                          >  2:696803/1‑98 (MQ=255)                                                                    aaTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCt                          >  2:812361/1‑100 (MQ=255)                                                                    aaTGTGTTCAATAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCt                          >  2:812422/1‑100 (MQ=255)                                                                             aaTAAAGCTCGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGgtag                 <  2:332557/100‑1 (MQ=255)                                                                                      cGAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGgcg        >  1:35893/1‑100 (MQ=255)                                                                                       gAACAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGgcgc       <  2:171674/100‑1 (MQ=255)                                                                                        aaCAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCa      <  2:1373802/100‑1 (MQ=255)                                                                                          cAAATTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCAtt    >  1:776460/1‑100 (MQ=255)                                                                                            aaTTCGTTCATTATGATCGGTCAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAc  <  1:937447/100‑1 (MQ=255)                                                                                            aaTTCGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAc  <  1:316489/100‑1 (MQ=255)                                                                                                |                                                                                                CTTCGCTCACTATACATCCTTGGCTGCAGGTTTAGTTGTACACCACTCCTAAATTTAATGTGTTGGCAATGTGTTCAATAAAGCTCGAACAAATTAGCTCATTATGATCGGTTAATACTTCAACTTCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTAC  >  NC_000913/263077‑263267

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