Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,339,445	A→G	Y95C (TAC→TGC)	fis →	DNA‑binding protein Fis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,339,445	0	A	G	98.6%	281.4 / ‑3.9	69	Y95C (TAC→TGC)	fis	DNA‑binding protein Fis
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (27/41);  total (27/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

GTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTT  >  REL606/3339373‑3339517                                                                         |                                                                         gTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCgg                                                                       >  1:211041/1‑76 (MQ=255)  tGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGc                                                                      <  1:744758/76‑1 (MQ=255)  tGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGc                                                                      <  1:1611137/76‑1 (MQ=255)  tGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGc                                                                      >  2:1268992/1‑76 (MQ=255)  tGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGc                                                                      >  2:1170278/1‑76 (MQ=255)    ntAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCAt                                                                    >  1:2188108/2‑76 (MQ=255)    gTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCAt                                                                    >  2:1499417/1‑76 (MQ=255)     tAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATg                                                                   <  2:761552/76‑1 (MQ=255)     tAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATg                                                                   <  2:586591/76‑1 (MQ=255)     tAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATg                                                                   >  2:2089707/1‑76 (MQ=255)      aaCCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGa                                                                  <  2:1103078/76‑1 (MQ=255)      aaCCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGa                                                                  <  2:1103069/76‑1 (MQ=255)      aaCCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGa                                                                  <  2:765029/76‑1 (MQ=255)       aCCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGCTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGaa                                                                 <  1:255865/76‑1 (MQ=255)         cAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACt                                                               >  1:2024854/1‑76 (MQ=255)            aCCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAAt                                                            <  2:1026415/76‑1 (MQ=255)             cccGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAAtt                                                           <  1:1577787/76‑1 (MQ=255)             cccGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAAtt                                                           <  2:1828051/76‑1 (MQ=255)               cGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCa                                                         >  2:1396722/1‑76 (MQ=255)                gTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAg                                                        >  2:495919/1‑76 (MQ=255)                 tgctgcGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAgg                                                       >  2:666826/1‑76 (MQ=255)                  gctgcGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGt                                                      >  1:701483/1‑76 (MQ=255)                   ctgcGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGtt                                                     <  2:769040/76‑1 (MQ=255)                   ctgcGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGtt                                                     <  1:2096284/76‑1 (MQ=255)                     gcgcTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAg                                                   <  1:1807243/76‑1 (MQ=255)                     gcgcTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAg                                                   >  1:923855/1‑76 (MQ=255)                       gcTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTGATTCAGGTTAGCt                                                 <  2:1254568/76‑1 (MQ=255)                        cTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTa                                                <  1:735695/76‑1 (MQ=255)                        cTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTa                                                <  1:1847417/76‑1 (MQ=255)                        cTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTa                                                <  2:267558/76‑1 (MQ=255)                         tgatgaTGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTaa                                               >  1:187683/1‑76 (MQ=255)                          gatgatGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTaaa                                              >  1:253221/1‑76 (MQ=255)                          gatgatGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTaaa                                              >  1:982676/1‑76 (MQ=255)                          gatgatGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTaaa                                              >  1:1838698/1‑76 (MQ=255)                           atgatgGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAAt                                             <  1:905831/76‑1 (MQ=255)                            tgatgGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATg                                            >  2:87871/1‑76 (MQ=255)                               tgGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCtt                                         <  2:2075931/76‑1 (MQ=255)                               tgGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCtt                                         >  2:1934055/1‑76 (MQ=255)                                gggCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTg                                        <  1:1884272/76‑1 (MQ=255)                                 ggCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGa                                       <  1:793568/76‑1 (MQ=255)                                    aTCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTa                                    <  2:1872227/76‑1 (MQ=255)                                     tCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTaa                                   <  2:1440495/76‑1 (MQ=255)                                       aaCCGTGGTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTaaaa                                 <  1:1140041/76‑1 (MQ=255)                                             ggTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGgcgc                           <  1:714147/76‑1 (MQ=255)                                              gTACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCt                          <  1:1039785/76‑1 (MQ=255)                                               tACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTa                         <  2:950454/76‑1 (MQ=255)                                               tACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTa                         >  2:167725/1‑76 (MQ=255)                                               tACGCTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTa                         >  2:167739/1‑76 (MQ=255)                                                   cTGCGTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCg                     >  2:533844/1‑76 (MQ=255)                                                       gTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCAt                 >  2:1796838/1‑76 (MQ=255)                                                       gTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCAt                 <  1:1376295/76‑1 (MQ=255)                                                       gTAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCAt                 >  1:2028163/1‑76 (MQ=255)                                                        tAAAAAATTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATg                <  2:1131631/76‑1 (MQ=255)                                                           aaaaTTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATgggg             >  2:1600960/1‑76 (MQ=255)                                                            aaaTTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGa            <  2:484584/76‑1 (MQ=255)                                                            aaaTTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGa            <  2:1393949/76‑1 (MQ=255)                                                             aaTTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGaa           >  2:584752/1‑76 (MQ=255)                                                             aaTTGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGaa           <  2:740492/76‑1 (MQ=255)                                                               ttGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAgc         <  2:583562/76‑1 (MQ=255)                                                                tGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAgcg        >  1:223304/1‑76 (MQ=255)                                                                tGAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAgcg        <  1:1523928/76‑1 (MQ=255)                                                                 gAAAAAATGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAgcgc       >  1:2122510/1‑76 (MQ=255)                                                                    aaaaTGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCtt    <  2:315940/76‑1 (MQ=255)                                                                    aaaaTGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCtt    <  2:1918963/76‑1 (MQ=255)                                                                    aaaaTGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCtt    <  2:767682/76‑1 (MQ=255)                                                                     aaaTGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCttt   <  2:132940/76‑1 (MQ=255)                                                                     aaaTGCGGCATGAACTAATTAAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCttt   <  1:1158737/76‑1 (MQ=255)                                                                      aaTGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCtttt  <  2:341735/76‑1 (MQ=255)                                                                      aaTGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCtttt  <  2:1889398/76‑1 (MQ=255)                                                                      aaTGCGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCtttt  <  2:1763085/76‑1 (MQ=255)                                                                         |                                                                         GTGGTAACCAGACCCGTGCTGCGCTGATGATGGGCATCAACCGTGGTACGCTGCGTAAAAAATTGAAAAAATACGGCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTT  >  REL606/3339373‑3339517

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