Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 933,782 G→A 8.7% A1630A (GCC→GCT)  CJ019_00881 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT2440933,7820GA 8.7% 136.8 / 3.1 46A1630A (GCC→GCT) CJ019_00881hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/19);  new base A (2/2);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGGCAGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAGCGTCGGTGAGGGTGACTTGTAAGGTGCTGCCATCGTTCACG  >  KT2440/933650‑933905
                                                                                                                                    |                                                                                                                           
cTGGGCAGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCg                                                                                                                             <  4:12985/133‑1 (MQ=255)
cTGGGCAGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCg                                                                                                                             >  3:12985/1‑133 (MQ=255)
      aGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGc                                                                                                                     >  5:201859/1‑135 (MQ=255)
      aGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGCTGACACCATCGGCTAg                                                                                                                        <  4:225917/132‑1 (MQ=255)
      aGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGCTGACACCATCGGCTAg                                                                                                                        >  3:225917/1‑132 (MQ=255)
           aCAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGCCACCATCGGCTAGGGCCAGGt                                                                                                                >  3:56589/1‑135 (MQ=255)
           aCAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGt                                                                                                                >  8:86036/1‑135 (MQ=255)
           aCAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGt                                                                                                                >  8:188878/1‑135 (MQ=255)
            cAGGGGGAGGCTGAACAGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTc                                                                                                               <  6:204978/135‑1 (MQ=255)
                gggAGGCTGAACCGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGt                                                                                                            <  5:131748/132‑1 (MQ=255)
                            aGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTgcgggc                                                                                               >  6:66755/1‑135 (MQ=255)
                            aGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTgcgggc                                                                                               >  5:84918/1‑135 (MQ=255)
                            aGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTgcgggc                                                                                               >  4:74276/1‑135 (MQ=255)
                            aGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGGGGTgcgggc                                                                                               >  5:259219/1‑135 (MQ=255)
                                 gTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGCCACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGgca                                                                                           >  6:16474/1‑134 (MQ=255)
                                                acCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTCAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGt                                                                           <  2:10903/135‑1 (MQ=255)
                                                   aGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGAc                                                                        >  2:122558/1‑135 (MQ=255)
                                                        cGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTc                                                                                    <  6:256165/118‑1 (MQ=255)
                                                        cGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTc                                                                                    >  5:256165/1‑118 (MQ=255)
                                                        cGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGc                                                                   >  1:208952/1‑135 (MQ=255)
                                                             aCAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGCGGTGACCGAGCCTGGt                                                              >  1:261869/1‑135 (MQ=255)
                                                             aCAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGtt                                                             >  7:248668/1‑136 (MQ=255)
                                                             aCAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCAGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGt                                                              >  5:197032/1‑135 (MQ=255)
                                                             aCAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCAGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGt                                                              >  5:197035/1‑135 (MQ=255)
                                                               aGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCTGCTAGGGCAAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCTGCCGAGGTGACCGAGCCTGGTTg                                                            >  3:261918/1‑135 (MQ=255)
                                                               aGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTgg                                                           <  7:86036/136‑1 (MQ=255)
                                                                  aCCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTggg                                                          >  4:77121/1‑134 (MQ=255)
                                                                       cgGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGCCCGCGCCGGGTTGGGATACGt                                                    >  5:263330/1‑135 (MQ=255)
                                                                        ggCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTg                                                  >  8:154968/1‑136 (MQ=255)
                                                                          cATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGc                                                 <  3:77121/135‑1 (MQ=255)
                                                                          cATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGc                                                 <  6:84918/135‑1 (MQ=255)
                                                                                gTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAg                                           <  4:138949/135‑1 (MQ=255)
                                                                                    cTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAGCGt                                        <  5:66755/134‑1 (MQ=255)
                                                                                         cGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAGCGTCGGTGa                                  >  2:51107/1‑135 (MQ=255)
                                                                                         cGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAGCGTCGGTGAg                                 <  7:94732/136‑1 (MQ=255)
                                                                                          ggTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGTAGCGTCGGTGAg                                 <  6:244761/135‑1 (MQ=255)
                                                                                          ggTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAGCGTCGGTGAg                                 <  2:156523/135‑1 (MQ=255)
                                                                                          ggTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAGCGTCGGTGAg                                 <  1:51107/135‑1 (MQ=255)
                                                                                               cgcCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGGTGGGATACGTTGCCGGCAGCGTCGGTGAGGGTGa                            >  2:33954/1‑135 (MQ=255)
                                                                                                  cGGGTTCGCCACGACCGGTGAAGGTGACACCATCAGCTAGGGCCCGGTCAGTCGGGTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAGCGTCGGTTAGGGTGACtt                         <  6:197035/135‑1 (MQ=255)
                                                                                                  cGGGTTCGCCACGACCGGTGAAGGTGACACCATCAGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAGCGTCGGTGAGGGTGACtt                         <  6:197032/135‑1 (MQ=255)
                                                                                                  cGGGGTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAGCGTCGGTGAGGGTGACtt                         <  6:201859/135‑1 (MQ=255)
                                                                                                         gCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGGCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAGCGTCGGTGAGGGTGACTTGTAAgg                   <  5:95360/134‑1 (MQ=255)
                                                                                                         gCCACGACCGGTGAAGGTGACACCATCGGCTAGGGACAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAGCGTCGGTGAGGGTGACTTGTAAGGtg                 >  7:31724/1‑136 (MQ=255)
                                                                                                                 cGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAGCGTCGGTGAGGGTGACTTGTAAGGTGCTGCCAt          >  6:72507/1‑135 (MQ=255)
                                                                                                                         gggACACCCTGGGGTAGGGCCAGGTCAGTGGGTGGGGGGGGCGGCAGCAAGGGGGCGGCGGTGACCGAGCCGGGGTGGGATAGGTTGCGGGCAGCGTGGGGGAGGGTGACTTTTAAGGTGCTGCCATCGTTCAcg  <  6:16442/133‑1 (MQ=255)
                                                                                                                         gTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAGCGTCGGTGAGGGTGACTTGTAAGGTGCTGCCATCGTTCAcg  <  2:94101/135‑1 (MQ=255)
                                                                                                                          tGACACCAGCGGCTAGGGCCAGGTCAGGGGGTTGGGCGGGCGGCAGCAAGGCGGCCGAGGTGACCGAGCCTGGTTGGGATAGGTTGCCGGCAGCGTCGGTGAGGGTGACTTGTAAGGTGCTGCCATCGTTCAcg  <  3:74276/134‑1 (MQ=17)
                                                                                                                                    |                                                                                                                           
CTGGGCAGGCGACAGGGTGAGGCTGAACAGGCCGTCGGCACCGACGACACCAGTACCGATGACAGCACCGCCGGCATCGCGTACCTGCACGGTAGCGCCGGGTTCGCCACGACCGGTGAAGGTGACACCATCGGCTAGGGCCAGGTCAGTCGGTTGTGCGGGCGGCAGCAAGTCGGCCGAGGTGACCGAGCCTGGTTGGGATACGTTGCCGGCAGCGTCGGTGAGGGTGACTTGTAAGGTGCTGCCATCGTTCACG  >  KT2440/933650‑933905

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