Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	5,393,134	G→C	26.2%	A434P (GCC→CCC)	CJ019_05075 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	5,393,134	0	G	C	26.2%	102.9 / 7.0	51	A434P (GCC→CCC)	CJ019_05075	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/8);  new base C (2/11);  total (32/19)
Fisher's exact test for biased strand distribution p-value = 1.04e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGACATACAGGCATTGAAGCACCAAGTGCTAAACCCCGAAGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCTTTGAGCCATTAAAGCTAGTGTTTGAGTTCATCCAGT  >  ACB122/5393000‑5393260                                                                                                                                       |                                                                                                                               ccGACATACAGGCATTGAAGCACCAAGTGCTAAACCCCGAAGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAg                                                                                                                                >  1:258028/1‑135 (MQ=255) ccGACATACAGGCATTGAAGCACCAAGTGCTAAACCCCGAAGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAg                                                                                                                                >  7:250924/1‑135 (MQ=255)       tACAGGCATTGAAGCACCAAGTGCTAAACCCCGAAGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGc                                                                                                                          >  5:65134/1‑135 (MQ=255)              aTTGAAGCACCAAGTGCTAAACCCCGAAGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCATGGCCAGTGGTTCCACAACAAAGCCACGCCCgctgg                                                                                                                   >  1:80844/1‑135 (MQ=255)                  aaGCACCAAGTGCTAAACCCCGAAGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCggtggtggt                                                                                                               >  7:147980/1‑135 (MQ=255)                  aaGCACCAAGTGCTAAACCCCGAAGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTCGCCAGTGGTTCCAGAAGAAACCCACGCCCggtggtggt                                                                                                               >  8:97627/1‑135 (MQ=255)                          aGTGCTAAACCCCGAAGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTGCAGAAGAAAGACACGCCCGGTGGTGGTGAGATa                                                                                                         >  7:143008/1‑133 (MQ=255)                           gTGCTAAACCCCGAAGGGAACCCGTGGGGAAAAGGCGCGCAATCTCACCTGACTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGACCAGTGGTTACAGAAGAAAGCCACGCCCGGTGGTGGTGAGATa                                                                                                         >  3:120610/1‑132 (MQ=255)                             gCTAAACCCCGAAGGGGAACCGTGGGGTAAAGGCGCGCAATCTCCCCTGCCTATTGAGTTGCTTTTTTACAAGGAAGAATTTCTTGGCCAGTGGTTACAGAAGAAAGCCCCGCCCGGGGGTGGTGAGATAATCGa                                                                                                    >  6:244521/1‑135 (MQ=255)                                  aCCCCGAAGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTa                                                                                               >  7:222399/1‑135 (MQ=255)                                   ccccGAAGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTAt                                                                                              >  6:313905/1‑135 (MQ=255)                                    cccGAAGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATc                                                                                             >  5:365683/1‑135 (MQ=255)                                         aGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCg                                                                                         >  1:102231/1‑134 (MQ=255)                                         aGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCg                                                                                         >  1:323177/1‑134 (MQ=255)                                         aGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGGGAGATAATCGAGTTTATCGGCg                                                                                         >  5:140154/1‑134 (MQ=255)                                                 cGTGGTGTAAAGGCGCCAAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAACAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTAt                                                                                <  7:134397/135‑1 (MQ=255)                                                        tAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTg                                                                          >  4:68981/1‑134 (MQ=255)                                                          acgccgcgcACTCTCACCTGCTTATTTAGTTGTTTTTATAATAGGAATAATTTGTTGGCCTGGGTTACCACAATAATCCCACGCCCGGGGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGcc                                                                        <  8:263316/130‑1 (MQ=255)                                                          accgcgcccAACCCCACCTCCCTATTGTGTTGCTTTTTAACACGGACCAATTTCTGGCCCTGGGGTCCCACACCACCCCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGcc                                                                        <  1:150342/131‑1 (MQ=255)                                                              cgcgcAATCTCACCTGCCTATTGAGTTGCTTTTTTTCAAGGACGAATTTCTTTGCCAGTGGTTACAGAAGAAAACCCCGCCCGGTGGTGGTTAGCTAATATAGTTTATCGGCGAGAAGGTATCATTTGCCaaaaa                                                                   >  1:54174/1‑135 (MQ=255)                                                              cgcgcAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCaaaaa                                                                   >  1:246879/1‑135 (MQ=255)                                                              cgcgcAAGCTCACCGGCCTATTGTGTTGCTTTTTAAAAAGGACTAATTTCTTGGCCAGTGGTGCCACAACAACCCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCaaaaa                                                                   <  3:125163/135‑1 (MQ=255)                                                                   aaTCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACttt                                                               >  1:299803/1‑134 (MQ=255)                                                                   aaTCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTg                                                              >  3:256434/1‑135 (MQ=255)                                                                   aaTCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTg                                                              >  4:75024/1‑135 (MQ=255)                                                                   aaTCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTg                                                              >  2:59012/1‑135 (MQ=255)                                                                   aaTCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGGGGTGAGATAATCGAGTTTATCGGCGACAAGGTTTCTTTTGCCAAAAAAttt                                                               >  7:4083/1‑134 (MQ=255)                                                                    aTCTCCCTTGCCTTTTGTGTTTATTTTTTACAAGCACTATTTTTTTGGCCTGTGTATCCAGAACAACCCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGt                                                             <  5:104706/135‑1 (MQ=255)                                                                    aTCTCCCCTGCCTTTTGTGTTGCTTTTTTACCCGGACGATTTTCTGGGCCAGTGGGACCAGAACACAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGt                                                             <  2:299803/135‑1 (MQ=255)                                                                    aTCTCACCTGCTTATTGTGTGGGTTTTTTACACGGCCTTTTTTCTGGGCCAGTGGTTCCACAAGAACCCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGt                                                             <  2:246879/135‑1 (MQ=255)                                                                          ccTGCCTATTGCGTTGCTTTTTTACAAGGACGAATTGCTTGCCCTGGGGTCCCAGAACAACCCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTAc                                                           <  3:333629/131‑1 (MQ=255)                                                                                tATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTa                                                   >  5:72487/1‑133 (MQ=255)                                                                                tATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACATTGTACCCACACTAg                                                  >  5:209329/1‑134 (MQ=255)                                                                                aaTTGAGTAGCTATTTAACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAACCCCCCGCCCGGTCGTGGTGAGATAAACGAGTTTATCGGCGACAACGTATATTTTGCCAAAAACTTTCTACACACACTa                                                   >  3:161238/2‑133 (MQ=255)                                                                                    gAGTTGCTTTTTTCCAAGCTCTTTTTTGTTGCCCAGTGGTACCAGACCAAAGCCACGCCCGGTGCTGGTGATATAATCGAGTTTATCGGCGACAATGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGc                                              <  6:184648/134‑1 (MQ=255)                                                                                    gAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATTATCGAGTTTATCGGCGAGAAGGTATCTTTTGCCCAAACCTTTGTGCCCACCCTAGAAGcc                                             >  3:165872/1‑135 (MQ=255)                                                                                     cgTTGATTTTTAACACGCACTAATTGCTTGGCCCGGGGTCCCAGACCAAACCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGcc                                             <  6:65134/133‑1 (MQ=255)                                                                                      gTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCcac                                           >  3:291976/1‑135 (MQ=255)                                                                                      gTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCca                                            >  1:90246/1‑134 (MQ=255)                                                                                           ttttttACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCt                                      >  7:101808/1‑135 (MQ=255)                                                                                              tttACATGTACGTATTTCTTGCCCAGTGGGTCACGAAGAAACCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCttt                                    <  4:77983/134‑1 (MQ=255)                                                                                              tttACAAGTACTATTTTCTTGCCCAGGGGTTCCAGAACAAACCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCttt                                    <  6:72487/134‑1 (MQ=255)                                                                                              tttACAAGGACTATTTTCTTGGCCAGTGGGTCCATAACAAACCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCttt                                    <  2:102231/134‑1 (MQ=255)                                                                                              ttaacacGGACGAATTTCTTGGCCTGGGGTTCCAGACCAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCttt                                    <  4:34209/131‑1 (MQ=255)                                                                                                  cAAGCACGAAATTCTTGGCCTGGGGTTCCAGAACAACGCCACGCCCGGGGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCTTTGAGc                                <  8:4083/134‑1 (MQ=255)                                                                                                       aCGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCTTTGAGCCAtt                            >  7:39036/1‑133 (MQ=255)                                                                                                         gAATTTCTTGGCCAGTGGTGCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCTTTGAGCCATTAAAg                        <  6:139396/135‑1 (MQ=255)                                                                                                               cTTGGCCAGGGGTTCCAGAAGAACGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCTTTGAGCCATTAAAGCTAgt                   <  4:108715/134‑1 (MQ=255)                                                                                                               cTTGCCCAGGGGTTCCACAACAACCCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCTTTGAGCCATTAAAGCTAgt                   <  4:111139/134‑1 (MQ=255)                                                                                                                       gTGGTTCCAGAACAAAGCCACGCCCGGTGGGGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCTTTGAGCCATTAAAGCTAGTGTTTGAGtt          <  8:250924/135‑1 (MQ=255)                                                                                                                           ttCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTAATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCTTTGAGCCATTAAAGCTAGTGTTTGAGTTCAt       >  3:112176/1‑134 (MQ=255)                                                                                                                               agaagaAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCTTTGAGCCATTAAAGCTAGTGTTTGAGTTCATCCAg   >  5:303094/1‑134 (MQ=255)                                                                                                                                gaagaaAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCTTTGAGCCATTAAAGCTAGTGTTTGAGTTCATCCAGt  >  5:163269/1‑134 (MQ=255)                                                                                                                                       |                                                                                                                               CCGACATACAGGCATTGAAGCACCAAGTGCTAAACCCCGAAGGGAAACCGTGGGGTAAAGGCGCGCAATCTCACCTGCCTATTGAGTTGCTTTTTTACAAGGACGAATTTCTTGGCCAGTGGTTCCAGAAGAAAGCCACGCCCGGTGGTGGTGAGATAATCGAGTTTATCGGCGACAAGGTATCTTTTGCCAAAAACTTTGTACCCACACTAGAAGCCACATGCTTTGAGCCATTAAAGCTAGTGTTTGAGTTCATCCAGT  >  ACB122/5393000‑5393260

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: ‑