Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 2,525,483 T→A 9.3% V52E (GTA→GAA)  CJ019_02380 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,525,4830TA 9.3% 134.2 / 3.2 43V52E (GTA→GAA) CJ019_02380hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/15);  new base A (2/2);  total (26/17)
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.95e-01

ATACACATGCGTTGTTCGCCCTGGGCGACTTGCGCGGTAAATTGGCCCAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGACATCAAGTTGAATGTGTATGGCGTAG  >  ACB143/2525348‑2525611
                                                                                                                                       |                                                                                                                                
aTACACATGCGTTGTTCGCCCTGGGCGACTTGCGCGGTAAATTGGCCCAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGt                                                                                                                                  <  1:3457/136‑1 (MQ=255)
     caTGCGTTGTTCGCCCTGGGCGACTTGCGCGGTAAATTGGCCCAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTc                                                                                                                              >  7:77996/1‑135 (MQ=255)
          gttgttCGCCCTGGGCGACTTGCGCGGTAAATTGGCCCAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGa                                                                                                                         >  6:102354/1‑135 (MQ=255)
            tgttCGCCCTGGGCGACTTGCGCGGTAAATTGGCCCAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGCAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGAc                                                                                                                        >  8:27388/1‑134 (MQ=255)
             gttCGCCCTGGGCGACTTGCGCGGTAAATTGGCCCAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTGCATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACaa                                                                                                                      <  5:85171/135‑1 (MQ=255)
             gttCGCCCTGGGCGACTTGCGCGGTAAATTGGCCCAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACaa                                                                                                                      <  3:232250/135‑1 (MQ=255)
              ttCGCCCTGGGCGACTTGCGCGGTAAATTGGCCCAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGAAGTCGATGACAagc                                                                                                                    >  6:338940/1‑136 (MQ=255)
              ttCGCCCTGGGCGACTTGCGCGGTAAATTGGCCCAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGAAGTCGATGACAag                                                                                                                     >  6:344978/1‑135 (MQ=255)
                 gCCCTGGGCGACTTGCGCGGTAAATTGGCCCAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATAGACACCGAGAGTGCACTGGTAGTCGATGACAagcagc                                                                                                                 >  6:155762/1‑136 (MQ=255)
                          gACTTGCGCTGTAAATTGGCCCCACTGTTTCAGGTCCGCTTCATCTATGTCACGGACCAGAGTCCCGAAGGCCTTTACATGGCCGAAATCGACAGCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAAc                                                                                                        <  5:55477/136‑1 (MQ=255)
                               gcgcGGTAAATTGGCCCAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGaaa                                                                                                   >  8:151190/1‑136 (MQ=255)
                               gcgcGGTAAATTGGCCCAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGGGCACTGGTAGTCGATGACAAGCAGCGCCATGAACTGaa                                                                                                    >  1:404188/1‑135 (MQ=255)
                                     tAAATTGGCCCAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGc                                                                                             <  4:357921/136‑1 (MQ=255)
                                               cAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTc                                                                                   <  2:397512/136‑1 (MQ=255)
                                                            gCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTAtt                                                                      <  2:168062/136‑1 (MQ=255)
                                                                        aTGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCcgcg                                                           >  3:400622/1‑135 (MQ=255)
                                                                         tGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCcg                                                             <  4:30419/132‑1 (MQ=255)
                                                                           tCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAgg                                                       <  4:227370/136‑1 (MQ=255)
                                                                                gAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGc                                                   >  7:234299/1‑135 (MQ=255)
                                                                                gAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGAAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCa                                                  <  5:338940/136‑1 (MQ=255)
                                                                                gAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGAAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCa                                                  <  5:344978/136‑1 (MQ=255)
                                                                                  aCAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAg                                                >  2:405380/1‑136 (MQ=255)
                                                                                      agTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTg                                             >  6:220242/1‑135 (MQ=255)
                                                                                              aGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATgcgc                                    >  2:354552/1‑136 (MQ=255)
                                                                                                        aTGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGACa                          >  8:15931/1‑136 (MQ=255)
                                                                                                        aTGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGACa                          >  7:36988/1‑136 (MQ=255)
                                                                                                             cGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGACATCAAg                     >  8:3471/1‑136 (MQ=255)
                                                                                                             cGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGACATCAAg                     >  6:103130/1‑136 (MQ=255)
                                                                                                             cGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGACATCAAg                     >  7:1909/1‑136 (MQ=255)
                                                                                                              gtaaTCGACACCGAGATTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGACATCAAGt                    <  6:273712/134‑1 (MQ=255)
                                                                                                               aaaTCGACACCGAGAGTGCACTGGTAGTCGATGACaa                                                                                                                      <  8:30802/37‑1 (MQ=255)
                                                                                                               aaaTCGACACCGAGAGTGCACTGGTAGTCGATGACaa                                                                                                                      >  7:30802/1‑37 (MQ=255)
                                                                                                               aaaTCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGa                            <  3:103535/127‑1 (MQ=255)
                                                                                                               aaaTCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGa                            >  4:103535/1‑127 (MQ=255)
                                                                                                               aaaTCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGACATCAAGtt                   >  4:340314/1‑136 (MQ=255)
                                                                                                                aaTCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGACATCAAGtt                   <  2:255124/135‑1 (MQ=255)
                                                                                                                  tCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGtg                                                                                <  3:104676/72‑1 (MQ=255)
                                                                                                                  tCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGtg                                                                                >  4:104676/1‑72 (MQ=255)
                                                                                                                  tCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGACATCAAGTTGaa                >  2:259732/1‑136 (MQ=255)
                                                                                                                  tCGACAACGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCGTCGAAGGTGGCAATCTGGAAATGCGATTTCGCGCCCTAATGTTGaa                >  2:390249/1‑136 (MQ=255)
                                                                                                                     acacCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGGCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGGGGCAAGCTGGAAATGCGCTTTCGCGACATCAAGTTGAAtgt             <  6:307176/136‑1 (MQ=255)
                                                                                                                         cGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGACATCAAGTTGAATGTGTAt         >  3:5731/1‑136 (MQ=255)
                                                                                                                             agTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGACAGTCGGCGACCACTTTCGTGCTGCCGTATTGCACAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGACATCAAGTTGAATGTGTATGGCg     <  4:48255/136‑1 (MQ=255)
                                                                                                                                gCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGACATCAAGTTGAATGTGTATGGCGTAg  >  4:9879/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                
ATACACATGCGTTGTTCGCCCTGGGCGACTTGCGCGGTAAATTGGCCCAACTGTTTCAGGGCCGCTTCATCTATGTCACCGAACAGAGTCCCGAAGGCCTGTACATGGCCGAAATCGACACCGAGAGTGCACTGGTAGTCGATGACAAGCAGCGCCTTGAACTGAAAGTCGGCGACCACTTTCGTGCTGCCGTATTGCCCAGCCGCGAAGGTGGCAAGCTGGAAATGCGCTTTCGCGACATCAAGTTGAATGTGTATGGCGTAG  >  ACB143/2525348‑2525611

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