Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	1,836,360	C→A	39.9%	S497I (AGT→ATT)	barA ←	Signal transduction histidine‑protein kinaseBarA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,836,360	0	C	A	39.9%	26.3 / 49.8	50	S497I (AGT→ATT)	barA	Signal transduction histidine‑protein kinaseBarA
Reads supporting (aligned to +/- strand):  ref base C (16/14);  new base A (10/9);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

CGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTG  >  ACB122/1836228‑1836485                                                                                                                                     |                                                                                                                              cgcgcCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTAATg                                                                                                                             >  8:238347/1‑135 (MQ=255) cgcgcCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGATGCGGTGGCAGCGCTTCGGCATCGTCATGAGCCGTGGGCAGGTTGAGGCTGATCCAGATCTGCGAACCTTCGCCCGGGGTACTg                                                                                                                             >  8:329026/1‑135 (MQ=255)    gcCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTAATGTCg                                                                                                                          >  6:174638/1‑135 (MQ=255)             tGGCCATCGACAAGGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGAtc                                                                                                                 <  8:192513/135‑1 (MQ=255)                            gcggcACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTAATGTCGACACCGATCTCGCCGCCCATCTGt                                                                                                  >  2:273420/1‑135 (MQ=255)                             cggcACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGtt                                                                                                 >  7:300384/1‑135 (MQ=255)                                  cggcggCCCAGCAATGGCTGCGGTGGCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTAATGTCGACACCGATCTCGCCGCCCATCTGTTCGa                                                                                              >  6:80050/1‑133 (MQ=255)                                                tGGCTGCGGTGGCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAAt                                                                              <  5:321594/135‑1 (MQ=255)                                                tGGCTGCGGTGGCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTAATGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAAt                                                                              <  3:330974/135‑1 (MQ=255)                                                          ggCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGcc                                                                      >  6:350054/1‑133 (MQ=255)                                                                gcTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTAATGTCGCCACCGATCTCGCCGCCCATCTGTTCGATCAGGCGGTTGGAAATCACCAGCCCCAGACCg                                                              >  2:78758/1‑135 (MQ=255)                                                                   tCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTAATGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTg                                                           >  2:188500/1‑135 (MQ=255)                                                                     ggCATCTTCATGAGCCTTGGGCAGGTGGAGGCTGATCCATAACTGCGAACCTTCGCCGGGGTTATTGGCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTgcc                                                         <  6:265118/135‑1 (MQ=255)                                                                      gCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTgccg                                                        >  5:346695/1‑135 (MQ=255)                                                                           gccATGAGCCTGGGGCAGGTTGCGGCTGACCCAGAACTGCGACCCTTCGCCCGGGGTACTGTCGACACCGATCGCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCCCCAGCCCCAGACCGGTGCCGCCggg                                                   <  7:255781/133‑1 (MQ=255)                                                                               tGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTgg                                               >  5:196583/1‑135 (MQ=255)                                                                               tGAGCCTTGGGCAGGTGGAGCCTGAGCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCCGACCGGGGCCGCCGGGCTg                                                <  6:346695/134‑1 (MQ=255)                                                                                      tGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTAATGTCGACACCGATCTCGCCGCCCATCTGTTCGAGCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCa                                        <  7:238347/135‑1 (MQ=255)                                                                                         gCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAg                                       >  5:265507/1‑133 (MQ=255)                                                                                            ggTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTAATGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGt                                  <  5:80050/135‑1 (MQ=255)                                                                                               gtaGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTAATGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGt                               <  3:21727/133‑1 (MQ=255)                                                                                                gAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTAATGTCGACACCGATCTCCCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCGGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTg                                >  2:361535/1‑133 (MQ=255)                                                                                                       aTCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTgg                       <  2:58236/135‑1 (MQ=255)                                                                                                        tCCAGACCTGCGAACCTGCGCCGGGGGTACTGTCGCCACCGATCTCGCCGCCCTGCTGTTCGATCAGGCGGTTGGAAATCACCAGCCCCCGCCCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGc                      <  5:176170/135‑1 (MQ=255)                                                                                                        tCCAGAACTGCGACCCTTCGCCCGGGGTAATGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGc                      <  4:280190/135‑1 (MQ=255)                                                                                                        tCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGc                      <  5:135508/135‑1 (MQ=255)                                                                                                        tCCAGAACAGCGAACCTTCGCCCGGGGTACTGGCGACACCGATCTCGCCGCCCATCTGTTCGATCCGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGc                      <  1:258946/135‑1 (MQ=255)                                                                                                        gccAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGCCCGGGGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGc                      <  7:363072/134‑1 (MQ=255)                                                                                                        ccaaGAACTGCGACCCTTCGCGCGGTGTAATGTCGACACCGATCTCGCCGCCCATCTGTGCGATCAGGCGCTTGGAAATCACCCGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGc                      <  1:273420/132‑1 (MQ=255)                                                                                                          cAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGa                                                                                                                   <  5:300355/40‑1 (MQ=255)                                                                                                          cAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGa                                                                                                                   >  6:300355/1‑40 (MQ=255)                                                                                                          cAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCgg                                            >  5:229223/1‑111 (MQ=255)                                                                                                          cAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCgg                                            <  6:229223/111‑1 (MQ=255)                                                                                                          cAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCt                     >  2:331889/1‑134 (MQ=255)                                                                                                          cAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCt                     >  5:345260/1‑134 (MQ=255)                                                                                                          cAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCt                     >  2:60626/1‑134 (MQ=255)                                                                                                          cAGAACTGCGAACCTTCGCCCGGGGTAATGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGc                      >  1:68254/1‑133 (MQ=255)                                                                                                          cAGAACTGCGAACCTTCGCCCGGGGTAATGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGc                      <  2:68254/133‑1 (MQ=255)                                                                                                              ctggCGACCCTCCGCCGGGGGTACTGCCGACCCCGATCTCCCCGCCCTTTTGTGCGAGCCGGCGTTGGTAAAGCCCACCCCCCAGACGGGTGCGCCGGGGGGGGCTGGCCAGTGAGTTGTCAGCCTGGCTGaa                  <  6:255594/130‑1 (MQ=255)                                                                                                                tGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGcc              >  8:166676/1‑135 (MQ=255)                                                                                                                   gAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTgg           <  6:196583/135‑1 (MQ=255)                                                                                                                     aacTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTTTTAGATCAGGCTCTTGGAAATCACCAGCCCCAGCCCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGTCCCGGa          >  6:24397/3‑134 (MQ=255)                                                                                                                     aCCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGCCCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGaa         >  3:290893/1‑135 (MQ=255)                                                                                                                     aCCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGa          >  3:565/1‑134 (MQ=255)                                                                                                                     aCCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGa          >  2:377985/1‑134 (MQ=255)                                                                                                                     aCCTTCGCCCGGGGTAATGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCTCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGa          >  7:127601/1‑134 (MQ=255)                                                                                                                     aCCTTCGCCCGGGGTAATGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGa          >  1:31839/1‑134 (MQ=255)                                                                                                                            cccGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGtgtg  <  6:345260/135‑1 (MQ=255)                                                                                                                            cccGGGGTAATGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGtgtg  <  1:202243/135‑1 (MQ=255)                                                                                                                            cccGGGGTAATGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGtgtg  <  7:37396/135‑1 (MQ=255)                                                                                                                                     |                                                                                                                              CGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCGGCATCGTCATGAGCCTTGGGCAGGTTGAGGCTGATCCAGAACTGCGAACCTTCGCCCGGGGTACTGTCGACACCGATCTCGCCGCCCATCTGTTCGATCAGGCGCTTGGAAATCACCAGCCCCAGACCGGTGCCGCCGGGCTGGCGGGCCAGTGAGTTGTCAGCCTGGCTGAAGGCCTGGAACAGTGTG  >  ACB122/1836228‑1836485

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