Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 2,630,169 A→G 8.7% D31G (GAC→GGC)  CJ019_02487 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,630,1690AG 8.7% 137.4 / 2.8 46D31G (GAC→GGC) CJ019_02487hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (14/28);  new base G (2/2);  total (16/30)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTCGTGAATCTCGCGCCGGGCGCTGAACAGGGAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGAGCAGCTGGAGGATGCTGCGCGTTCGCTGGGGACCCTTACCACACCGCT  >  ACB143/2630036‑2630301
                                                                                                                                     |                                                                                                                                    
cATTCGTGAATCTCGCGCCGGGCGCTGAACAGGGAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACa                                                                                                                                    <  5:117049/136‑1 (MQ=255)
cATTCGTGAATCTCGCGCCGGGCGCTGAACAGGGAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACa                                                                                                                                    <  2:256444/136‑1 (MQ=255)
      tGAATCTCGCGCCGGGCGCTGAACAGGGAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAgg                                                                                                                              >  5:129386/1‑136 (MQ=255)
        aaTCTCGCGCCGGGCGCTGAACAGGGAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCATTACCGACAGCGACGcc                                                                                                                            >  2:75038/1‑136 (MQ=255)
        aaTCTCGCGCCGGGCGCTGAACAGGGAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGcc                                                                                                                            >  3:262371/1‑136 (MQ=255)
         aTCTCGCGCCGGGCGCTGAACAGGGAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCt                                                                                                                           <  8:8584/136‑1 (MQ=255)
           ctcGCGCCGGGCGCTGAACAGGGAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTg                                                                                                                         >  1:290771/1‑136 (MQ=255)
             cgcgCCGGGCGCTGAACAGGGAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCa                                                                                                                       >  2:162352/1‑136 (MQ=255)
             cgcgCCGGGCGCTGAACAGGGAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCa                                                                                                                       >  7:211219/1‑136 (MQ=255)
              gcgcCGGGCGCTGAACAGGGAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTg                                                                                                                         <  5:287632/133‑1 (MQ=255)
                tccGGGCTCTGAACAGGGAGCATCGTATGTCTGCCATGCTAGACCATGTAGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTTGCCTGTCGTGAACGCCTGGCTGCCAATACCGGCAGCCAGGCATTGCACGa                                                                                                                    >  8:28314/2‑136 (MQ=255)
                gcCGGGAGCTGAACAGGGAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGGCAGCGAGGCCTTGCACGa                                                                                                                    >  8:32402/1‑136 (MQ=255)
                     gcgcTGAACCGGGAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAAAACCGACAGCGAGGCCTTGCACGACTTGc                                                                                                               <  2:292332/136‑1 (MQ=255)
                                ggAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTTGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTg                                                                                                    <  5:275245/136‑1 (MQ=255)
                                     aTCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTgcggcg                                                                                               <  4:73618/136‑1 (MQ=255)
                                        gTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGcc                                                                                             <  6:262776/135‑1 (MQ=255)
                                           tGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTgc                                                                                          <  6:211566/135‑1 (MQ=255)
                                             tCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTgcgc                                                                                        <  5:132475/135‑1 (MQ=255)
                                             tCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTgcgc                                                                                        <  8:131354/135‑1 (MQ=255)
                                                         gACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAg                                                                             <  6:190431/134‑1 (MQ=255)
                                                                        cAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGc                                                            >  2:295609/1‑136 (MQ=255)
                                                                        cAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGc                                                            >  5:198021/1‑136 (MQ=255)
                                                                                 aCCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTg                                                    <  8:186972/135‑1 (MQ=255)
                                                                                 aCCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTg                                                    <  4:74235/135‑1 (MQ=255)
                                                                                 aCCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTg                                                    <  3:178688/135‑1 (MQ=255)
                                                                                 aCCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTg                                                    <  3:25063/135‑1 (MQ=255)
                                                                                 aCCTTGCAGGTGCGCTTGCTGGCCTGTCGGGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGt                                                     <  6:121044/134‑1 (MQ=255)
                                                                                   cTTGCCGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGag                                                 <  1:73591/136‑1 (MQ=255)
                                                                                   cTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGGCAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGag                                                 <  7:32402/136‑1 (MQ=255)
                                                                                   cTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGag                                                 >  3:100280/1‑136 (MQ=255)
                                                                                   cTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGag                                                 <  1:9475/136‑1 (MQ=255)
                                                                                   cTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGATGCCAATACCGGCAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGAGCCTGCGCAGTCTGCTGAGGCCTCTGCGTGGGCTGCCTGGCGTGGag                                                 <  7:28314/136‑1 (MQ=255)
                                                                                       ccgttGCGCTTGAGGGCCTGTCGTGAACGCCTGTCTGCCAATACCTACAGCGAGGCCTTGCACTACTTGCGTACCACCTTGCGTCGCCTGCGCAGTCTGATGAGTCCGCTGCGTGGGCTGCCTGGCGTGGatcagc                                             <  2:172140/132‑1 (MQ=255)
                                                                                        aGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGAGCAGCt                                            >  1:13542/1‑136 (MQ=255)
                                                                                        aGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAAAACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGAGCAGCt                                            >  5:57925/1‑136 (MQ=255)
                                                                                         ggTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACCGCGAGGCCTTGCACGACATGCGTACCCCCATGCCGCGCCTGCGCAGGCAGCTGAGGCCGCAGCGTGGGCTGCCTGGCGGGGAGCAGCTg                                           >  4:43921/1‑136 (MQ=255)
                                                                                                ttGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGAGCAGCTGGAGGAtg                                    >  3:91241/1‑136 (MQ=255)
                                                                                                             gTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGAGCAGCTGGAGGATGCTGCGCGTTCGc                        <  6:262150/135‑1 (MQ=255)
                                                                                                             gTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGAGCAGCTGGAGGATGCTGCGCGTTCGc                        <  4:22295/135‑1 (MQ=255)
                                                                                                              tGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGAGCAGCTGGAGGATGCTGCGCGTTCGCTg                      <  2:130073/136‑1 (MQ=255)
                                                                                                                aaCGCCTGGCTGCCAATACCGACCGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGGCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGAGCAGCTGGAGGATGCTGCGCGTTCGCTggg                    <  8:106669/136‑1 (MQ=255)
                                                                                                                aaCGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGAGCAGCTGGAGGATGCTGCGCGTTCGCTggg                    <  4:91241/136‑1 (MQ=255)
                                                                                                                aaCGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGAGCAGCTGGAGGATGCTGCGCGTTCGCTggg                    <  2:13542/136‑1 (MQ=255)
                                                                                                                   gCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGAGCAGCTGGAGGATGCTGCGCGTTCGCTGGGGAc                 <  3:48767/136‑1 (MQ=255)
                                                                                                                                  ccGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGAGCAGCTGGAGGATGCTGCGCGTTCGCTGGGGACCCTTACCACACCGCt  <  6:110520/136‑1 (MQ=255)
                                                                                                                                  ccGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGACTGGCGTGGAGCAGCTGGAGGATGCTGCGCGTTCGCTGGGGACCCTTACCACACCGCt  >  4:39815/1‑136 (MQ=255)
                                                                                                                                     |                                                                                                                                    
CATTCGTGAATCTCGCGCCGGGCGCTGAACAGGGAGCATCGTATGTCTGCCATGCTCGACCATGTCGTTGCCCAGATACTGACCTTGCAGGTGCGCTTGCTGGCCTGTCGTGAACGCCTGGCTGCCAATACCGACAGCGAGGCCTTGCACGACTTGCGTACCACCTTGCGGCGCCTGCGCAGTCTGCTGAGGCCGCTGCGTGGGCTGCCTGGCGTGGAGCAGCTGGAGGATGCTGCGCGTTCGCTGGGGACCCTTACCACACCGCT  >  ACB143/2630036‑2630301

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