Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,210,811	C→T	pseudogene (1712/2001 nt)	yehQ →	pseudogene

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,210,811	0	C	T	98.6%	247.6 / NA	71	pseudogene (1712/2001 nt)	yehQ	pseudogene
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (37/33);  minor base A (1/0);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCT  >  NC_000913/2210714‑2210906                                                                                                  |                                                                                                gTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGc                                                                                               >  2:684874/1‑100 (MQ=255)   cAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTa                                                                                             >  2:217055/1‑100 (MQ=255)     gTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTAcc                                                                                           <  1:289370/100‑1 (MQ=255)     gTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTAcc                                                                                           <  1:732557/100‑1 (MQ=255)      ttttGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCt                                                                                          >  2:697909/1‑100 (MQ=255)         tgttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCt                                                                                       >  1:190364/1‑100 (MQ=255)         tgttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCt                                                                                       >  1:581766/1‑100 (MQ=255)         tgttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCt                                                                                       >  1:1352636/1‑100 (MQ=255)           ttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGc                                                                                     <  2:738147/100‑1 (MQ=255)              aaGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAg                                                                                  >  2:370194/1‑100 (MQ=255)                      ttGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGc                                                                          <  1:431606/100‑1 (MQ=255)                       tGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGcc                                                                         <  1:557214/100‑1 (MQ=255)                       tGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGcc                                                                         <  1:1142046/100‑1 (MQ=255)                        gACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGGTTTGCCt                                                                        >  1:687765/1‑100 (MQ=255)                          cacTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtgt                                                                      >  1:1290835/1‑100 (MQ=255)                           acTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtgtg                                                                     <  2:1491658/100‑1 (MQ=255)                           acTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtgtg                                                                     <  2:1396195/100‑1 (MQ=255)                           acTGGAACCACGGGGGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtgtg                                                                     >  2:232537/1‑100 (MQ=255)                            cTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtgtcc                                                                    <  1:1068422/100‑3 (MQ=255)                            cTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGc                                                                    >  1:446710/1‑100 (MQ=255)                            cTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGc                                                                    >  1:527269/1‑100 (MQ=255)                            cTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCAGTGc                                                                    <  1:1257433/100‑1 (MQ=255)                                aaCCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTc                                                                <  1:1487284/100‑1 (MQ=255)                                   cACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTAc                                                             <  1:378380/100‑1 (MQ=255)                                   cACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTAc                                                             <  1:161036/100‑1 (MQ=255)                                   cACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTAc                                                             <  1:1087351/100‑1 (MQ=255)                                   cACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTAc                                                             >  1:659649/1‑100 (MQ=255)                                      gggTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGc                                                          <  1:1002397/100‑1 (MQ=255)                                        gTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtc                                                        >  2:1350281/1‑100 (MQ=255)                                           atgatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtcatc                                                     >  2:1206869/1‑100 (MQ=255)                                             gatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAg                                                   <  2:92565/100‑1 (MQ=255)                                             gatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAg                                                   >  2:18380/1‑100 (MQ=255)                                              atgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGt                                                  >  1:1494197/1‑100 (MQ=255)                                                gaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTg                                                >  1:1217675/1‑100 (MQ=255)                                                gaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTg                                                >  2:372131/1‑100 (MQ=255)                                                gaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTg                                                >  2:319447/1‑100 (MQ=255)                                                 aCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGc                                               <  1:1526730/100‑1 (MQ=255)                                                  cAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCt                                              >  1:973582/1‑100 (MQ=255)                                                  cAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCt                                              >  2:476304/1‑100 (MQ=255)                                                    aaaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGa                                            >  2:1052321/1‑100 (MQ=255)                                                     aaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGAt                                           <  2:1474758/100‑1 (MQ=255)                                                       aCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATAc                                         <  2:581764/100‑1 (MQ=255)                                                             gCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGcc                                   <  1:476305/100‑1 (MQ=255)                                                              ccTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCa                                  >  1:711572/1‑100 (MQ=255)                                                               cTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAg                                 <  2:712671/100‑1 (MQ=255)                                                                 gggCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGc                               >  2:127956/1‑100 (MQ=255)                                                                 gggCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGc                               >  2:128302/1‑100 (MQ=255)                                                                  ggCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCg                              <  1:501901/100‑1 (MQ=255)                                                                    cgcTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGtt                            >  2:1303078/1‑100 (MQ=255)                                                                     gcTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTa                           <  1:128302/100‑1 (MQ=255)                                                                     gcTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTa                           <  1:127956/100‑1 (MQ=255)                                                                     gcTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTa                           <  2:1561942/100‑1 (MQ=255)                                                                     gcTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTa                           >  1:1467259/1‑100 (MQ=255)                                                                      cTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTAc                          <  2:90471/100‑1 (MQ=255)                                                                           cGCAGAAACTACGCCTGTGGCACAGATACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTAtt                     >  2:485014/1‑100 (MQ=255)                                                                             cAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCa                   >  1:1397944/1‑100 (MQ=255)                                                                             cAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCa                   <  1:24863/100‑1 (MQ=255)                                                                             cAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCAGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCa                   >  2:876809/1‑100 (MQ=255)                                                                                    tACGCCTGTGGCACTGCTACCTTCTGCCATCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTc            <  2:759198/100‑1 (MQ=255)                                                                                      cGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCa          <  2:1292872/100‑1 (MQ=255)                                                                                      cGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCa          >  1:454172/1‑100 (MQ=255)                                                                                      cGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCa          >  1:522917/1‑100 (MQ=255)                                                                                      cGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCa          >  1:898362/1‑100 (MQ=255)                                                                                      cGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCa          >  2:164999/1‑100 (MQ=255)                                                                                        ccTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCGGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATa        <  2:659647/100‑1 (MQ=255)                                                                                        ccTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATa        <  1:370194/100‑1 (MQ=255)                                                                                           gtgGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACg     >  2:543167/1‑100 (MQ=255)                                                                                             ggCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGc   <  2:1352629/100‑1 (MQ=255)                                                                                             ggCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGc   <  2:190364/100‑1 (MQ=255)                                                                                             ggCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGc   >  1:206535/1‑100 (MQ=255)                                                                                              gCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCt  >  1:1232795/1‑100 (MQ=255)                                                                                                  |                                                                                                GTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCT  >  NC_000913/2210714‑2210906

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