Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	194,448	G→A	7.5%	D3103N (GAC→AAC)	CJ019_00185 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	194,448	0	G	A	7.5%	167.4 / 2.6	53	D3103N (GAC→AAC)	CJ019_00185	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/24);  new base A (2/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGG  >  ACB143/194314‑194583                                                                                                                                       |                                                                                                                                        ccGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGa                                                                                                                                        <  1:161578/136‑1 (MQ=255)     cACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCt                                                                                                                                     >  6:163543/1‑135 (MQ=255)         caCGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCTACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCg                                                                                                                                >  6:210055/2‑136 (MQ=255)           cgagcgagGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGAGCTCGGTAGCCTAGAACGTCAAACCGACCTTGACTGTGCACGTCAATCAGCCGCTTGCCCACGAACTGATCGTg                                                                                                                              <  4:267184/136‑1 (MQ=14)           cGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTg                                                                                                                              <  5:10297/136‑1 (MQ=255)           cGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTg                                                                                                                              <  4:89199/136‑1 (MQ=255)                cAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCt                                                                                                                         <  5:10339/136‑1 (MQ=255)                          aGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCaacaac                                                                                                               <  3:206366/136‑1 (MQ=255)                               aCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCt                                                                                                            >  2:204771/1‑134 (MQ=255)                               aCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCa                                                                                                          <  3:206910/136‑1 (MQ=255)                                 gAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAgg                                                                                                        <  5:242655/136‑1 (MQ=255)                                  aaGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGt                                                                                                       <  7:208969/136‑1 (MQ=255)                                   aGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACAACCTGa                                                                                                                                   <  4:289303/107‑1 (MQ=255)                                   aGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACAACCTGa                                                                                                                                   <  4:284854/107‑1 (MQ=255)                                   aGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACAACCTGa                                                                                                                                   >  3:289303/1‑107 (MQ=255)                                   aGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACAACCTGa                                                                                                                                   >  3:284854/1‑107 (MQ=255)                                    ggCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCa                                                                                                     >  4:56638/1‑136 (MQ=255)                                          aTCGTCAAGGTCACCATCACGGCCGACCAGTCCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCa                                                                                               <  3:250474/136‑1 (MQ=255)                                          aTCGTCAAGGGCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCa                                                                                               <  8:204294/136‑1 (MQ=255)                                           tCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACAGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCaa                                                                                              <  3:164482/136‑1 (MQ=255)                                                      aCCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGGCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGaaa                                                                                   <  6:186565/136‑1 (MQ=255)                                                         aTCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCa                                                                                <  7:195698/136‑1 (MQ=255)                                                                    ccAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTa                                                                      >  7:87664/1‑135 (MQ=255)                                                                      aGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTacac                                                                   >  2:197222/1‑136 (MQ=255)                                                                           tCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCac                                                               >  1:174943/1‑135 (MQ=255)                                                                           tCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCac                                                               >  2:255674/1‑135 (MQ=255)                                                                             ggTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCacgac                                                            <  1:165978/136‑1 (MQ=255)                                                                                aGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGAcgcg                                                         >  3:69483/1‑136 (MQ=40)                                                                                aGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGAcgcg                                                         >  1:154726/1‑136 (MQ=40)                                                                                aGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGAcgcg                                                         >  3:11462/1‑136 (MQ=40)                                                                                     agaACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAgg                                                    >  6:223518/1‑136 (MQ=40)                                                                                     agaACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACAAGCGCGCCGTACACCCACGACGCGCAAgg                                                    >  3:193592/1‑136 (MQ=255)                                                                                     agaACGTCAAACCGACCTTCAATGTGAAAGTCCAACAGCCGCTGGCACACGCCCTGATCGTGAGCCTGCGCAGCAGCGGGCAGGTCACCATCAAGGCCGGGGAAGCCAGCGCGCCGGGCGCCGACGACgcgctcgg                                                    >  4:240199/1‑132 (MQ=2)                                                                                      gaACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGc                                                   >  7:178094/1‑136 (MQ=40)                                                                                       aaCGTCAAACCGACCTTCACTGTGCACGTCAATCCGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCg                                                  <  2:94806/136‑1 (MQ=255)                                                                                       aaCGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGc                                                   >  8:186498/1‑135 (MQ=40)                                                                                       aaCGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCg                                                  <  1:179483/136‑1 (MQ=40)                                                                                          gTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGAt                                                >  7:224339/1‑135 (MQ=40)                                                                                               aCCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTc                                          >  2:223648/1‑136 (MQ=40)                                                                                                 cGACCTTCAGTGTGCACGTCAATCAGCCGCGGGCCCACGACCTGTTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTTAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTa                                        <  2:154726/136‑1 (MQ=255)                                                                                                        cACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAgg                                   <  7:18913/134‑1 (MQ=255)                                                                                                         aCTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACg                                <  6:257453/136‑1 (MQ=255)                                                                                                                aCGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGc                           >  3:196566/1‑134 (MQ=255)                                                                                                                aCGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGc                           >  7:75471/1‑134 (MQ=255)                                                                                                                     aaTCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCTCGCCGTACACCCACTACGCGCAAGGCGATCACGTCAATAAGGACGATGGCCatat                      >  3:27517/1‑134 (MQ=14)                                                                                                                     aaTCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGAt                      >  4:271184/1‑134 (MQ=255)                                                                                                                     aaTCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCa                    <  4:193592/136‑1 (MQ=255)                                                                                                                        cAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCAGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCt                 >  3:138450/1‑136 (MQ=255)                                                                                                                          gCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTCCACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCtgg               >  5:182225/1‑134 (MQ=255)                                                                                                                          gCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTg               >  1:163025/1‑136 (MQ=40)                                                                                                                          gCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTg               <  7:44454/136‑1 (MQ=40)                                                                                                                                 gCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCa        <  6:219821/136‑1 (MQ=40)                                                                                                                                 gCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCa        <  7:52255/136‑1 (MQ=40)                                                                                                                                       gACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTcgg  <  3:291119/136‑1 (MQ=40)                                                                                                                                       |                                                                                                                                        CCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAACTCGG  >  ACB143/194314‑194583

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