Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 3,959,315 T→G 42.5% intergenic (+76/‑138) CJ019_03739 → / → rlhA_2 hypothetical protein/23S rRNA 5‑hydroxycytidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24403,959,3150TG42.5% 42.8 / 11.6 49intergenic (+76/‑138)CJ019_03739/rlhA_2hypothetical protein/23S rRNA 5‑hydroxycytidine synthase
Reads supporting (aligned to +/- strand):  ref base T (5/22);  new base G (19/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.07e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCACTGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGAAAATA  >  KT2440/3959182‑3959445
                                                                                                                                      |                                                                                                                                  
ttCACTGGCCTTTGACGGCTTGCGGGGGGAAGTGAGGGCGAGGGCGGCCAAGCGCTGAGCGGGGCTTGGGGGGGGGCGCAGCGCCCCGACAGCGGTGATTGGTCGCAAAATCGCGGGG‑GCCGCGCTGCGGCCAtt                                                                                                                                   <  8:111770/135‑1 (MQ=255)
ttCACGGGCCTTGGACGGCATGCGGGGGGAAGTGCCGGCTCTGGCGGCAAAGGGCTGAGCTGCGGTTGTGGGGGTGCGCAGCGCCCCGACATCGCGGTTTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCtt                                                                                                                                   <  8:104571/135‑1 (MQ=255)
 tCACTGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCttt                                                                                                                                  >  4:348875/1‑135 (MQ=255)
  cACTGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGGCCGTTc                                                                                                                                 >  8:53699/1‑135 (MQ=255)
     tGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCGTTcgcg                                                                                                                              >  6:100260/1‑135 (MQ=255)
       gCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTTCCGCCCGTTAGCGAc                                                                                                                            >  8:273704/1‑135 (MQ=255)
             gAAGGCATGCGTGGGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGGGGCGCTTGGGGGGGGGCGCAGCGCCCCGACATCGCTGATTGGTCGAAAAATCGCCGGG‑GCCGCTCTGCGGCCCTTACGCGACGCTa                                                                                                                        <  8:235552/133‑1 (MQ=255)
               aGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GGCGCGCTGCGGCCCGTTCGGCCCGCGAGGcc                                                                                                                    >  1:89392/1‑135 (MQ=255)
                 gccggCGGGTGCAAGTCCCGGCGATGGCGGCCATGCGCTGAGCTGCGCTTGTGGGGGGGCGCAGCGCCCCGACAGCGCTGATGGGACCCAAAATCGCCGGG‑GCCGCGCTGCGGCCCTTTCGCGACGCTAGGCCGc                                                                                                                  <  3:140368/131‑1 (MQ=255)
                   ttGCGTGGGGAAGGGCCGGCGATGGCGGCCACGCGCTGAGCTGCGCTGGGGGGGGGGCGCAGCGCCCCGACTTCGCTGATTGGCCGAAAAAGCGCCGGG‑GCCGCTCGGCGGCCCTTGCGCGACGCTAGGCCGCTc                                                                                                                <  6:254728/134‑1 (MQ=255)
                    tGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTcc                                                                                                               >  4:174908/1‑135 (MQ=255)
                      cGTGTGGATGTGATGGCTTTGCATGCCTAGCGTGGAGCTGAGCTTGGGTGGGTGGGCAGCGCCCAGACTGGGCTGATT‑GTAACAAAATTCCCGGGTGCCGCTGTGGGGCCCTTTCGCGACGCTAGGGCGCTCCca                                                                                                             <  8:367671/135‑1 (MQ=255)
                      cGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCca                                                                                                             <  8:315292/135‑1 (MQ=255)
                           ggAAGTGTCGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCGTTCGCGCCGCGAGGCCGCTCCCCCACtt                                                                                                        >  7:104172/1‑135 (MQ=38)
                               gTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCTTTAGCGCCGCGAGGGCGCGCCCACAGCTgtgt                                                                                                    >  5:67687/1‑132 (MQ=255)
                                ggccGGCGAGGGCTGCCACGCGCTGAGCGGCGCTTGGGGGGGGGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGcgcg                                                                                                   <  8:117708/132‑1 (MQ=38)
                                  aCGGCGATGGCTGCCAAGCGCGGAGCGGCGCTTGGGGGGGGGCGCAGCGCCCCGACATCGCTGATTGGGCGCAAAAGCGCCGGG‑GCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTa                                                                                                 <  1:33706/135‑1 (MQ=38)
                                     gCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GGCGCACAGCGGCCCGTTAGGGAGCCTAGTCCGCTCCCCCACAAGCGCGGAccc                                                                                              >  1:300727/1‑135 (MQ=255)
                                     gCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCGGTCGCGCCGCAAGGCCGCGCCCCCACTTGcgcgggacc                                                                                              >  3:245553/1‑130 (MQ=38)
                                        ttGGCGGCCACGCGCGGCGGGGCGCTTGGGGGGGGGCGCCGCGCCCCGACATCGCTGATTGGGCGAAAAATCGCGGGG‑GCCGCGCGGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCa                                                                                           <  2:295063/134‑1 (MQ=255)
                                         gggcggcCACGCGCTGCGGGGGGCGTGGGGGGGTGCGCAGCGCCCGGACTTCGCTGTTTGGTCGAAAAACCCCCGGG‑CCCGCGCGGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAg                                                                                          <  4:68968/134‑1 (MQ=255)
                                            cTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GGCGCGCTGCGGCCCGGTCGCGCCGCGAGGGCGCGCCCACCCTTGCGCGGACCCCCAGcc                                                                                        >  1:295063/1‑134 (MQ=255)
                                               ccACGCGCTGCGCTGCGGTTGGGGGGGGGCCCCGCGCCCCGAAAGCGCTGATTGGGCCAAAAACCGCCGGG‑GCCGCTCTGCGGCCCTTGCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATc                                                                                    <  2:16837/135‑1 (MQ=25)
                                                      ctgcgctgcgTTTGTGGGGGGGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCAcc                                                                             <  4:245553/135‑1 (MQ=255)
                                                         aggggcgcTGGGGGGGGGGCTCAGGGCCCCGACAGCGCTTTTGGGTCGCAATAGCGCGGGG‑GCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACctct                                                                          <  2:143707/131‑1 (MQ=25)
                                                            tGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCCCGCTTGGGCCCGTTAGCGCCGCGTGGGCGCGCCCCCCCGTGCGCGGTCCCGCCGCCCATCCCCCCCCgctccg                                                                       >  5:157961/1‑133 (MQ=255)
                                                               gcTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCGGCGGCCCGTTCGCGCCGCGAGGGCGCGCCCACAGTTTCGCGTGCCCGCAGCCAAACCCGCCCCGCGCGGccc                                                                    >  8:80630/1‑133 (MQ=255)
                                                                   gtgGGGGTGCGCAGCGCCCCGACATCGCTGATTGCTCGCAAAATCGCCGGG‑GCCGCGCTGCGGCCCCTTCGCGCCCCGACGCCGCTCCCACACTGGCGCGGACCCCCACCCAATCACGCCCCCCCCCGCAcccc                                                                 >  2:304237/1‑134 (MQ=255)
                                                                     ggggggGCGCAGCGCCCGGACTGCGCTGTTTGGGCCCAAATTCGCCGGG‑GCCGCTGTGCGGCCCGTGCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAat                                                              <  7:330572/135‑1 (MQ=37)
                                                                     gggggTGCGCAGCGCCCCGACATCGCTGATGGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAat                                                              <  3:300830/135‑1 (MQ=255)
                                                                     gggggTGCGAAGCGCCCTGACATCGCTTATTGTTCCAAAAATCGCCGGG‑GCCGCTGTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCaa                                                               <  6:131346/134‑1 (MQ=38)
                                                                      ggggTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCGTTGGCGACGCGAGGCCGGGCCCCCCGTTGCGCGTGCCCGCCGCCAAACCCGCCCCTCGCGGcccccc                                                                 >  6:154273/1‑131 (MQ=255)
                                                                         gTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GGCGCGCAGCGGCCCGTTAGCGACGGCAGGCCGCGCCCACACTTTCGCGGACCCGCCCCCAATACCGCCCCGCTCGGCCCCCCCaggatg                                                          >  4:66144/1‑135 (MQ=255)
                                                                           gcgcAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCGTTCGCGCCGCTAGGCCGCCCCCCCAGTTGCGCGTGCCCGCCGCCAATCCCGCCCCGATCGGCCCCCCAATGCTgcc                                                        >  4:290573/1‑135 (MQ=37)
                                                                           gcgcAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCGCTGCGGCCCTTTCCCGTCGCTAGGCCGCTCCCACAGTTTCGCGGACCCGCCCCCAATCCCCCCCCTCTCCGCACCCCAAGGATgcc                                                        >  8:317473/1‑135 (MQ=255)
                                                                              cAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCGTTTGCGGCGCGAGGGCGCTCCCCCACTTTGGCGGCCCCGCCGCCCAGCCCGCCCCGCTCGGCCCCCCCatgatggggcc                                                     >  4:376456/1‑130 (MQ=255)
                                                                              cAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCGTTAGCGCCGCGAGGCCGCTCCCACAGTTGCGCGGGCCCGCAGCCCATCCCGCCCCGCGCGGCCCCCCCATGCTtccgcc                                                     >  8:76480/1‑135 (MQ=25)
                                                                              cAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCGCGGCGGCCCGTTCGCGCCGCGCGGCCGCGCCCCCCCTTGCGCGTCCCCGCCGCCCCTCCCGCACCGCTCGGCCCCCCCCTGGGgccgcc                                                     >  8:104928/1‑135 (MQ=255)
                                                                                    cccGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCGTTCGCGCCGCGAGGCCGCTGCCCCCGTTGCGCGTGCCCGCAGGCCATCCCGCCCCGCGCCGCCCCCCCATGATGCCGCCCCCCCg                                               >  8:242483/1‑135 (MQ=25)
                                                                                    cccGACAGCGCGGATTGGGCGCAAAATCGCCGGG‑GCCGCTCGGCGGCCCTTGCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCg                                               <  4:264275/135‑1 (MQ=255)
                                                                                       taCATCGGTGTTTGGTCGCAAAATCGGCGGG‑GCCGCTGGGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGttt                                            <  3:41131/134‑1 (MQ=38)
                                                                                                 aTTGGTCGCAAAATCGCCGGG‑GCCGCGCTGCGGGCCGTTAGGGCCGCGAGGGCGCGCCCCCACTTTCGCGGACCCCCCGCCAAACCCGCCCCACACGGCCCCCCCCAGATGGCCCCACCCCGGTTGAccccaac                                   >  4:113909/1‑129 (MQ=255)
                                                                                                   tGGGCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAAt                                <  3:179493/135‑1 (MQ=255)
                                                                                                         aaaacgcGCCGGG‑CCCGCGCGGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGc                          <  5:181569/129‑1 (MQ=255)
                                                                                                          aaaaTCGCCGGG‑GCCGCTCTGCGGCCCGTTCGCGACCCTAGGCCGCGCCCCCACTTGCGCGGTCCCGCCGCCAATCCCGCACCTCTCGGCACCCCCCTGATGCCCCCACCCCGGTTCACCCCTACAAAGccccc                          >  6:84999/1‑132 (MQ=25)
                                                                                                                 cgggg‑gCCGGTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCg                  <  7:273704/133‑1 (MQ=255)
                                                                                                                          gCTGTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGcc          <  3:174908/135‑1 (MQ=255)
                                                                                                                                 ggCCGTTCCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGAAAat   <  5:154273/135‑1 (MQ=255)
                                                                                                                                 ggCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGAAAat   <  4:315793/135‑1 (MQ=255)
                                                                                                                                 ggCCCGTCCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGAAAat   <  7:53699/135‑1 (MQ=255)
                                                                                                                                   cccTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGAAAata  >  6:357739/1‑134 (MQ=255)
                                                                                                                                   cccTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGGTTCACCGCTACAATGCCCGCCTAAACCGGGACAGCCCGAAAata  >  5:273821/1‑134 (MQ=255)
                                                                                                                                      |                                                                                                                                  
TTCACTGGCCTTTGAAGGCATGCGTGTGGAAGTGACGGCGATGGCTGCCAAGCGCTGAGCTGCGCTTGTGGGGGTGCGCAGCGCCCCGACATCGCTGATTGGTCGCAAAATCGCCGGG‑GCCGCTCTGCGGCCCTTTCGCGACGCTAGGCCGCTCCCACAGTTGCGCGTACCCGCAGCCAATCCCGCACCTCTCGGCACCCCAATGATGCCGCCAGCCCGTTTCACCGCTACAATGCCCGCCTAAACCGTGACAGCCCGAAAATA  >  KT2440/3959182‑3959445

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