Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,847,724	(C)6→7	coding (285/1044 nt)	selD ←	selenophosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,847,718	1	.	C	97.8%	149.7 / ‑2.4	45	coding (291/1044 nt)	selD	selenophosphate synthase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (23/21);  total (23/22)
Fisher's exact test for biased strand distribution p-value = 4.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGACATGCATAGCGTCCACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑C‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGT  >  NC_000913/1847619‑1847810                                                                                                     |                                                                                              acgacATGCATAGCGTCCACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTg                                                                                                <  1:1025780/100‑1 (MQ=255)   gacATGCATAGCGTCCACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑c‑c                                                                                            >  1:225907/1‑100 (MQ=255)   gacATGCATAGCGTCCACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑c                                                                                              <  2:1238275/99‑1 (MQ=255)     gatctaTAGCGTCCACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑c‑ccc                                                                                          <  1:269181/95‑1 (MQ=255)      aTGCATAGCGTCCACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑c‑cccc                                                                                         <  2:1599574/100‑1 (MQ=255)       tGCATAGCGTCCACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑c‑ccccc                                                                                        >  2:228485/1‑100 (MQ=255)        gCATAGCGTCCACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑c‑cccccc                                                                                       >  2:1806723/1‑99 (MQ=255)          aTAGCGTCCACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑c‑cccagatc                                                                                     >  2:269182/1‑95 (MQ=255)              cGTCCACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATcgcg                                                                                  <  1:1699444/100‑1 (MQ=255)              cGTCCACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATcgcg                                                                                  >  2:979220/1‑100 (MQ=255)              cGTCCACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATcgcg                                                                                  >  2:348762/1‑100 (MQ=255)                tCCACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGaa                                                                                >  2:91415/1‑100 (MQ=255)                 ccACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAg                                                                               <  1:1625969/100‑1 (MQ=255)                 ccACCTTCGGTCACTACGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAg                                                                               <  2:1788262/100‑1 (MQ=255)                       tCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCg                                                                         >  1:1434550/1‑100 (MQ=255)                        cGGTCACTTCGCTGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGc                                                                        >  1:202359/1‑100 (MQ=255)                         ggTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCt                                                                       >  2:1166918/1‑100 (MQ=255)                          gTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCt                                                                       <  2:1705380/99‑1 (MQ=255)                           tCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGa                                                                     <  1:5178/100‑1 (MQ=255)                            cACTTCGTGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGAt                                                                    <  2:1021940/100‑1 (MQ=255)                                 cgcgGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGt                                                               <  2:1199453/100‑1 (MQ=255)                                      gCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTc                                                          >  2:82951/1‑100 (MQ=255)                                        aaTTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTcgc                                                        <  1:203737/100‑1 (MQ=255)                                         aTTTCTGGGGAAAGGTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTcgcc                                                       <  1:1151618/100‑1 (MQ=255)                                           ttCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTcgccgc                                                     <  1:979212/100‑1 (MQ=255)                                             cTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCaa                                                   >  1:995658/1‑100 (MQ=255)                                               ggggAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATg                                                 >  2:1063723/1‑100 (MQ=255)                                               ggggAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCACaa                                                   >  2:180276/1‑98 (MQ=255)                                               ggggAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGAAGATAGTCGCCGCAAt                                                  >  2:824912/1‑99 (MQ=255)                                                     aGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCa                                           <  1:867224/100‑1 (MQ=255)                                                      gTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCATTGCGGcgac                                          <  2:286696/100‑4 (MQ=255)                                                      gTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCaa                                          >  1:1540089/1‑100 (MQ=255)                                                      gTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAAATATCGCTGATGGAGTTAATAGCCACAATACGGCCaa                                          >  2:278451/1‑100 (MQ=255)                                                            tAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGa                                    >  2:1662875/1‑100 (MQ=255)                                                              tctggCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGaaa                                  <  2:69122/97‑1 (MQ=255)                                                              aTCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGaaa                                  <  2:137509/100‑1 (MQ=255)                                                                cGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAgg                                >  1:1143187/1‑100 (MQ=255)                                                                 ggCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGaaaaga                               >  1:69122/1‑97 (MQ=255)                                                                  gCCAGCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATAGAAAGGAt                              >  2:235402/1‑100 (MQ=255)                                                                      gCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATc                          >  2:2596/1‑100 (MQ=255)                                                                      gCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATc                          <  1:841678/100‑1 (MQ=255)                                                                      gCCGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATc                          >  2:616178/1‑100 (MQ=255)                                                                        cGAGGATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCaa                        <  1:1479441/100‑1 (MQ=255)                                                                           ggATCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGa                     <  2:1044598/100‑1 (MQ=255)                                                                             aTCGCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATc                   >  2:1547555/1‑100 (MQ=255)                                                                                gCAATCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGc                >  2:990884/1‑100 (MQ=255)                                                                                   aTCGCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATa             <  2:1434561/100‑1 (MQ=255)                                                                                      gCCATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAaag          >  1:295655/1‑100 (MQ=255)                                                                                       ccATAATCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAaaga         <  2:905301/100‑1 (MQ=255)                                                                                         aTAATCGGTTTG‑CTCTCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAaagaag       <  2:225908/100‑1 (MQ=255)                                                                                           aaTCGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACAATCGGCATAAAGAAGTc     >  2:1676235/1‑100 (MQ=255)                                                                                             tCGGTTTCCC‑CCCCCATCGCGAACATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCg    <  1:1481616/99‑1 (MQ=255)                                                                                              cGGTTTGCC‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCggt  <  1:348760/100‑1 (MQ=255)                                                                                                     |                                                                                              ACGACATGCATAGCGTCCACCTTCGGTCACTTCGCGGGCAATTTCTGGGGAAAGTTTGTTAATCGGCCAGCCGAGGATCGCAATCGCCATAATCGGTTTG‑C‑CCCCCATCGCGAAGATATCGCTGATGGCGTTAGTCGCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGT  >  NC_000913/1847619‑1847810

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