Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,968,911	T→C	M23V (ATG→GTG)	rppH ←	RNA pyrophosphohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,968,911	0	T	C	98.4%	210.2 / ‑4.3	61	M23V (ATG→GTG)	rppH	RNA pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (26/34);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

TCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACATTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGATTGTTTCATACTCCAGGAAGACGGTAAACCACT  >  NC_000913/2968780‑2969046                                                                                                                                    |                                                                                                                                        tCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGcc                                                                                                                                 >  2:893814/1‑140 (MQ=255)  cGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGccc                                                                                                                                <  2:268802/140‑1 (MQ=255)        tCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGTAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCg                                                                                                                          >  2:11332/1‑140 (MQ=255)         cTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGa                                                                                                                         >  1:786183/1‑140 (MQ=255)           ttGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGAtt                                                                                                                       <  1:380702/140‑1 (MQ=255)            tGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTa                                                                                                                      <  2:228782/140‑1 (MQ=255)                   tAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCa                                                                                                               >  1:859765/1‑140 (MQ=255)                   tAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCa                                                                                                               <  1:581310/140‑1 (MQ=255)                     aTCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACg                                                                                                             <  1:397316/140‑1 (MQ=255)                     aTCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACg                                                                                                             <  1:1006347/140‑1 (MQ=255)                      tCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGa                                                                                                            <  2:208100/140‑1 (MQ=255)                           cttcttCAAACAATTCACGGTACATCGCCTGCTCTGTGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAAt                                                                                                                 <  2:492279/130‑1 (MQ=255)                           cttcttCAAACAATTCACGGTACATCGCCTGCTCTGTGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAAt                                                                                                                 >  1:492279/1‑130 (MQ=255)                           cttcttCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCt                                                                                                       >  2:5238/1‑140 (MQ=255)                              cttcAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACATTACCTGCCCCTGGCGATTACAAATCACGATACCTACg                                                                                                    >  2:1093109/1‑140 (MQ=255)                              cttcAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACg                                                                                                    >  1:500982/1‑140 (MQ=255)                                tcAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACa                                                                                                                    >  2:854352/1‑122 (MQ=255)                                tcAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACa                                                                                                                    <  1:854352/122‑1 (MQ=255)                                 cAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGttt                                                                                                 >  2:324464/1‑140 (MQ=255)                                   aaCAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATTCCTACGTTTgg                                                                                               <  2:113744/140‑1 (MQ=255)                                     cAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTg                                                                                                >  2:403089/1‑137 (MQ=255)                                     cAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTg                                                                                                <  1:403089/137‑1 (MQ=255)                                     cAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCCACGTTTGGGc                                                                                             <  1:159649/140‑1 (MQ=255)                                        ttCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCg                                                                                            >  2:634090/1‑138 (MQ=255)                                        ttCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCg                                                                                            <  1:634090/138‑1 (MQ=255)                                        ttCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGt                                                                                          >  2:107204/1‑140 (MQ=255)                                             ggTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCa                                                                                     <  2:500982/140‑1 (MQ=255)                                              gTACATCGCCTGCTCTGCGGATTCTCCGGGGTTAATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCAt                                                                                    <  2:191475/140‑1 (MQ=255)                                                aCATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCg                                                                                  <  1:43402/140‑1 (MQ=255)                                                  aTCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCAt                                                                                    <  1:46166/136‑1 (MQ=255)                                                  aTCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCAt                                                                                    >  2:46166/1‑136 (MQ=255)                                                    cGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGtcat                                                                              >  2:132649/1‑140 (MQ=255)                                                     gCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGtcatc                                                                             <  1:324422/140‑1 (MQ=255)                                                     gCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGtcatc                                                                             <  2:859765/140‑1 (MQ=255)                                                       cTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCatcaa                                                                           >  1:1063723/1‑140 (MQ=255)                                                         gCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCatcaatc                                                                         <  2:84311/140‑1 (MQ=255)                                                         gCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCatcaatc                                                                         >  2:5105/1‑140 (MQ=255)                                                         gCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCatcaatc                                                                         >  2:457127/1‑140 (MQ=255)                                                          ctctGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCatcaatca                                                                        >  2:71409/1‑140 (MQ=255)                                                             tGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCg                                                                     <  1:71409/140‑1 (MQ=255)                                                                      tcCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTc                                                            >  2:754159/1‑140 (MQ=255)                                                                       ccGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCa                                                           >  2:278704/1‑140 (MQ=255)                                                                            gTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATa                                                      <  1:893814/140‑1 (MQ=255)                                                                             ttGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATaa                                                     <  2:1063723/140‑1 (MQ=255)                                                                                aTTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGc                                                  >  1:343803/1‑140 (MQ=255)                                                                                aTTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGc                                                  >  1:860008/1‑140 (MQ=255)                                                                                aTTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGc                                                  >  2:33771/1‑140 (MQ=255)                                                                                 ttCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCt                                                 <  1:754159/140‑1 (MQ=255)                                                                                 ttCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCt                                                 <  2:653180/140‑1 (MQ=255)                                                                                               aaTTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGa                                   <  2:343803/140‑1 (MQ=255)                                                                                               aaTTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGa                                   <  2:860008/140‑1 (MQ=255)                                                                                                         aGTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGATTGTTTCATa                         >  1:885525/1‑140 (MQ=255)                                                                                                          gTGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGATTGTTTCATAc                        <  2:885525/140‑1 (MQ=255)                                                                                                           tGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGATTGTTTCa                           <  2:547300/136‑1 (MQ=255)                                                                                                           tGCTGACCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGATTGTTTCa                           >  1:547300/1‑136 (MQ=255)                                                                                                                aCCAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGATTGTTTCATACTCCAgg                  >  1:871255/1‑140 (MQ=255)                                                                                                                  cAAATCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGATTGTTTCATACTCCAGGaa                <  1:5238/140‑1 (MQ=255)                                                                                                                     aTCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGATTGTTTCATACTCCAGGAAGAc             >  1:709406/1‑140 (MQ=255)                                                                                                                      tCGCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGATTGTTTCATACTCCAGGAAGACg            <  1:5105/140‑1 (MQ=255)                                                                                                                        gCCGGGCCCACACTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGATTGTTTCATACTCCAGGAAGACGGt          <  2:786183/140‑1 (MQ=255)                                                                                                                                cacaCTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGATTGTTTCATACTCCAGGAAGACGGTAAACCAct  <  1:107204/140‑1 (MQ=255)                                                                                                                                cacaCTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGATTGTTTCATACTCCAGGAAGACGGTAAACCAct  <  1:33771/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                        TCGAACGTCTTTGCGGCTTAATCCTACTTCTTCAAACAATTCACGGTACATCGCCTGCTCTGCGGATTCTCCGGGGTTGATTCCGCCTTGCGGAAATTGCCAGGAGTGCTGACCAAATCGCCGGGCCCACATTACCTGCCCCTGGCGATTACAAATCACGATACCTACGTTTGGGCGGTAGCCATCGTCATCAATCACCGGACTACCTCAAAATAAAGCTTTATATACGAATGATTGTTTCATACTCCAGGAAGACGGTAAACCACT  >  NC_000913/2968780‑2969046

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