Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,346,140	T→C	F426L (TTT→CTT)	nrdA →	ribonucleoside‑diphosphate reductase 1, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,346,140	0	T	C	98.6%	253.6 / NA	72	F426L (TTT→CTT)	nrdA	ribonucleoside‑diphosphate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (37/34);  minor base A (1/0);  total (38/34)
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

GAAAGCCGTTGAGCTGTTCTCGCTGATGATGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGTTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTG  >  NC_000913/2346043‑2346236                                                                                                  |                                                                                                 gAAAGCCGTTGAGCTGTTCTCGCTGATGATGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGctt                                                                                                >  2:653240/1‑100 (MQ=255) gAAAGCCGTTGAGCTGTTCTCGCTGATGATGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGctt                                                                                                <  1:557202/100‑1 (MQ=255)         ttGAGCTGTTCTCGCTGATGATGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGc                                                                                        <  2:2191006/100‑1 (MQ=255)         ttGAGCTGTTCTCGCTGATGATGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGc                                                                                        <  2:2190393/100‑1 (MQ=255)         ttGAGCTGTTCTCGCTGATGATGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGc                                                                                        <  2:150645/100‑1 (MQ=255)         ttGAGCTGTTCTCGCTGATGATGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGc                                                                                        <  1:147586/100‑1 (MQ=255)            aGCTGTTCTCGCTGATGATGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCAt                                                                                     >  2:1365823/1‑100 (MQ=255)            aGCTGTTCTCGCTGATGATGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCAt                                                                                     >  2:1202142/1‑100 (MQ=255)              cTGTTCTCGCTGATGATGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATcg                                                                                   <  2:2563141/100‑1 (MQ=255)                  tctcGCTGATGATGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGcc                                                                               <  1:1175544/100‑1 (MQ=255)                  tctcGCTGATGATGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGcc                                                                               >  2:997888/1‑100 (MQ=255)                    tcGCTGATGATGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAg                                                                             >  2:812072/1‑100 (MQ=255)                       cggttgatgCAGGAACGTGCTGCTACCTGTCGGATCCGGATTCAGAACGGTGACGAGTGCAATACCCACGGGCCGCTGGATCCGGCCATCGCGCCAGTGc                                                                          <  2:516496/96‑1 (MQ=255)                         gatgatGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGt                                                                        <  2:263917/100‑1 (MQ=255)                          atgatgCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTc                                                                       <  1:825512/100‑1 (MQ=255)                            gatgCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAg                                                                     <  1:823727/100‑1 (MQ=255)                            gatgCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAg                                                                     >  1:1233706/1‑100 (MQ=255)                            gatgCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAg                                                                     <  2:2522100/100‑1 (MQ=255)                              tgCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTc                                                                   >  1:1931652/1‑100 (MQ=255)                               gCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCt                                                                  >  1:34251/1‑100 (MQ=255)                                cAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTa                                                                 >  1:2569947/1‑100 (MQ=255)                                cAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTa                                                                 >  1:468668/1‑100 (MQ=255)                                cAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTa                                                                 <  1:501697/100‑1 (MQ=255)                                   gAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCGTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAAcc                                                              <  2:1825273/100‑1 (MQ=255)                                     aCGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCtg                                                            >  2:936435/1‑100 (MQ=255)                                      cGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCtgt                                                           >  1:2075002/1‑100 (MQ=255)                                      cGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCtgt                                                           >  1:882425/1‑100 (MQ=255)                                       gTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTGGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCtgtg                                                          >  1:2525531/1‑100 (MQ=255)                                        tGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGc                                                         >  1:431814/1‑100 (MQ=255)                                         gCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGcc                                                        <  1:528567/100‑1 (MQ=255)                                           gTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTg                                                      >  1:1312378/1‑100 (MQ=255)                                              tACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGgag                                                   >  2:1982712/1‑100 (MQ=255)                                              tACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGgag                                                   >  2:1982711/1‑100 (MQ=255)                                                 cGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATa                                                >  2:674242/1‑100 (MQ=255)                                                       tATCTATATTCAGAACGTTGACCACTGCAATACCCATGGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTg                                          <  1:755632/100‑1 (MQ=255)                                                         tccatatTCAGAATTTTGACCACTGCAATACCCACAGCCCGCGTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCGTGGAGATAGCCCGGcc                                        <  2:916223/97‑1 (MQ=255)                                                         tCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGCGATAGCCCTGcc                                        >  1:1482379/1‑100 (MQ=255)                                                         tCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGcc                                        >  2:473029/1‑100 (MQ=255)                                                          cTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCg                                       <  2:351108/100‑1 (MQ=255)                                                           tataTTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGa                                      >  1:2223619/1‑100 (MQ=255)                                                           tataTTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGa                                      >  1:1245624/1‑100 (MQ=255)                                                            atatTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGAc                                     >  2:574137/1‑100 (MQ=255)                                                              atTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCa                                   <  1:1902755/100‑1 (MQ=255)                                                              atTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCa                                   >  1:2476511/1‑100 (MQ=255)                                                              atTCAGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCa                                   >  1:2506533/1‑100 (MQ=255)                                                               ttCAGAACGTTGACCACTGCAATACCCATACCCCGCTTGATCCGGCCACCGCGCCAGTGCGTCAGTCTAACCTGTGTCTGGAGATAGCCCTGCCGACCaa                                  <  1:2518166/100‑1 (MQ=255)                                                                  aGAACGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAAcc                               >  2:2517604/1‑100 (MQ=255)                                                                    aaCGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGc                             <  2:867472/100‑1 (MQ=255)                                                                    aaCGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGc                             <  1:1254591/100‑1 (MQ=255)                                                                     aCGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCt                            <  1:1176469/100‑1 (MQ=255)                                                                     aCGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCt                            <  1:1177023/100‑1 (MQ=255)                                                                      cGTTGTCCACTGCAATACCCATAGCCCGCTTGATCCGGCCAGCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTg                           <  2:1662822/100‑1 (MQ=255)                                                                      cGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTg                           <  2:8512/100‑1 (MQ=255)                                                                      cGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTg                           <  2:465115/100‑1 (MQ=255)                                                                      cGTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTg                           <  1:9721/100‑1 (MQ=255)                                                                       gTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCGCTGCCGACCCAACCGCTGa                          >  2:1867213/1‑100 (MQ=255)                                                                       gTTGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGa                          >  1:50005/1‑100 (MQ=255)                                                                        ttGACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGaa                         <  1:112749/100‑1 (MQ=255)                                                                          gACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAacg                       >  2:2153585/1‑100 (MQ=255)                                                                          gACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAacg                       >  2:532333/1‑100 (MQ=255)                                                                          gACCACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAacg                       <  1:2258062/100‑1 (MQ=255)                                                                            ccACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAacgac                     >  2:2018695/1‑100 (MQ=255)                                                                             cACTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAacgacn                    <  1:2094404/100‑2 (MQ=255)                                                                              aCTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGt                   <  1:1376163/100‑1 (MQ=255)                                                                               cTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTc                  <  1:266609/100‑1 (MQ=255)                                                                               cTGCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTc                  <  2:881392/100‑1 (MQ=255)                                                                                 gCAATACCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCaa                >  1:1634620/1‑100 (MQ=255)                                                                                     tACCCATAGCCCGCTTGATCCGGCCATCGCGCCAATGCGTCAGTCTAACCTGTACCTGGAGATAGCCCTGCCGAACAAACCGCTGAACGACGTCAacgac            >  1:2453252/1‑100 (MQ=255)                                                                                      aCCCATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAacgacg           >  2:1536384/1‑100 (MQ=255)                                                                                       cccATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTCACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAAcgacga          >  2:1508474/1‑100 (MQ=255)                                                                                       cccATAGCCCGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAAcgacga          >  1:2440909/1‑100 (MQ=255)                                                                                              cccGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGt   >  1:777201/1‑100 (MQ=255)                                                                                               ccGCTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTg  <  2:521534/100‑1 (MQ=255)                                                                                               ccGATTGATCCGGCCATCGCGCCAGTGCGTCCGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGCCGTCAACGACGAGAACGGTg  >  2:2383297/1‑100 (MQ=255)                                                                                                  |                                                                                                 GAAAGCCGTTGAGCTGTTCTCGCTGATGATGCAGGAACGTGCGTCTACCGGTCGTATCTATATTCAGAACGTTGACCACTGCAATACCCATAGCCCGTTTGATCCGGCCATCGCGCCAGTGCGTCAGTCTAACCTGTGCCTGGAGATAGCCCTGCCGACCAAACCGCTGAACGACGTCAACGACGAGAACGGTG  >  NC_000913/2346043‑2346236

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