Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 3,900,666 C→T 65.7% A508V (GCC→GTC)  CJ019_03690 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,900,6660CT65.7% 33.9 / 32.6 35A508V (GCC→GTC) CJ019_03690hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (7/5);  new base T (15/8);  total (22/13)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATGGGCGATTCGCGGCGCGTGGCCTGGGACAAGACCACCCACGAGTTCCCCGCCATCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGCCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGAACAACTGGACCTGCCCTGGGCACCGCTGAACCCGTTCAA  >  ACB143/3900531‑3900795
                                                                                                                                       |                                                                                                                                 
aaaTGGGCGATTCGCGGCGCGTGGCCTGGGACAAGACCACCCACGAGTTCCCCGCCATCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGc                                                                                                                                   >  1:43504/1‑136 (MQ=255)
     ggCGATTCGCGGCGCGTGGCCTGGGACAAGACCACCCACGAGTTCCCCGCCATCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGtcatca                                                                                                                              >  1:293198/1‑136 (MQ=255)
        gATTCGCGGCGCGTGGCCTGGGACAAGACCACCCACGAGTTCCCCGCCATCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCa                                                                                                                           <  2:266670/136‑1 (MQ=255)
        gATTCGCGGCGCGTGGCCTGGGACAAGACCACCCACGAGTTCCCCGCCATCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGCCATCACCa                                                                                                                           <  4:286912/136‑1 (MQ=255)
             gcggcgCGTGGCCTGGGACAAGACCACCCACGAGTTCCCCGCCATCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCg                                                                                                                      <  2:261241/136‑1 (MQ=255)
                            ggACAAGACCACCCACGAGTTCCCCGCCATCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCt                                                                                                        <  6:115121/135‑1 (MQ=255)
                                          acGAGTTCCCCGCCATCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCCGGCGCCTGGGCCGCGAGTtgct                                                                                         >  5:164323/1‑136 (MQ=255)
                                          acGAGTTCCCCGCCATCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTtgct                                                                                         >  7:54561/1‑136 (MQ=255)
                                          acGAGTTCCCCGCCATCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTtgct                                                                                         >  2:123883/1‑136 (MQ=255)
                                                    cGCCATCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGCCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGtggc                                                                               >  2:249612/1‑136 (MQ=255)
                                                     gCCATCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGtggct                                                                              >  5:228763/1‑136 (MQ=255)
                                                       cATCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGtggctgg                                                                            >  2:183425/1‑136 (MQ=255)
                                                         tCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGCCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTggctggct                                                                          <  2:130007/136‑1 (MQ=255)
                                                                cccGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGt                                                                   <  4:236104/136‑1 (MQ=255)
                                                                 ccGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTg                                                                  <  2:281017/136‑1 (MQ=255)
                                                                 ccGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTg                                                                  <  8:194021/136‑1 (MQ=255)
                                                                        gCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGtt                                                           >  8:116419/1‑136 (MQ=255)
                                                                        gCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGtt                                                           >  1:272910/1‑136 (MQ=255)
                                                                                   ggCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTTGCTGGCTGAACAGTGAGCAGTTTGTGGCCAccc                                                >  4:39500/1‑136 (MQ=255)
                                                                                           ccTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGaa                                        >  8:209075/1‑136 (MQ=255)
                                                                                           ccTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGaa                                        >  5:118147/1‑136 (MQ=255)
                                                                                           ccTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCCTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGACGaa                                        >  6:280285/1‑136 (MQ=255)
                                                                                           ccTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGCCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGaa                                        >  2:3830/1‑136 (MQ=255)
                                                                                           ccTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGCCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGaa                                        >  7:187014/1‑136 (MQ=255)
                                                                                           ccTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGCCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGaa                                        >  7:231114/1‑136 (MQ=255)
                                                                                              ggtggCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGCCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGtacaa                                     >  6:9833/1‑136 (MQ=255)
                                                                                                tggCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGAACAACt                                   >  3:222287/1‑136 (MQ=255)
                                                                                                  gCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGCCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGAACAACTg                                  <  5:222435/135‑1 (MQ=255)
                                                                                                  gCCAAGGGCCTGATCACAGACGAACAGTTGCAGCAGGCCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCTTGGCTGTCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCAGAACAACTg                                  <  7:41412/135‑1 (MQ=255)
                                                                                                      aGGGCCTGATCACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGAACAACTGGACCt                             <  3:21582/136‑1 (MQ=255)
                                                                                                      aGGGCCTGATCACCGACGAACAGTTGCAGCAGGCCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGAACAACTGGACCt                             >  2:3571/1‑136 (MQ=255)
                                                                                                           cTGATCACCGACGAACAGTTGCAGCAGGCCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGAACAACTGGACCTGCCCt                        <  6:124686/136‑1 (MQ=255)
                                                                                                               tCACCGACGAACAGTTGCAGCAGGCCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGAACAACTGGACCTGCCCTggg                     >  1:257058/1‑135 (MQ=255)
                                                                                                                cACCGACGAACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGAACAACTGGACCTGCCCTGGGCa                   >  1:247122/1‑136 (MQ=255)
                                                                                                                      cgaACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGAACAACTGGACCTGCCCTGGGCACCGCt              >  4:92771/1‑135 (MQ=255)
                                                                                                                      cgaACAGTTGCAGCAGGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGAACAACTGGACCTGCCCTGGGCACCGCt              >  8:235335/1‑135 (MQ=255)
                                                                                                                                 agcagGTCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGAACAACTGGACCTGCCCTGGGCACCGCTGAACCCGTTCaa  <  7:209075/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                 
AAATGGGCGATTCGCGGCGCGTGGCCTGGGACAAGACCACCCACGAGTTCCCCGCCATCGCCACCCCGGCTCGCACACCACTGGGCCAGCGCCTGGTGGCCAAGGGCCTGATCACCGACGAACAGTTGCAGCAGGCCATCACCAGCCCGGTTCGCAGGCGCCTGGGCCGCGAGTTGCTGCTGCGTGGCTGGCTGAACAGTGAGCAGTTGGTGGCCACCCTGGCCGAACAACTGGACCTGCCCTGGGCACCGCTGAACCCGTTCAA  >  ACB143/3900531‑3900795

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