Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,891,999 A→C 34.7% intergenic (+13/‑507) PP_RS08755 → / → napA tRNA‑Asp/Periplasmic nitrate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,891,9990AC34.7% 13.9 / 2.1 35intergenic (+13/‑507)PP_RS08755/napAtRNA‑Asp/Periplasmic nitrate reductase
Reads supporting (aligned to +/- strand):  ref base A (18/4);  new base C (5/8);  total (23/12)
Fisher's exact test for biased strand distribution p-value = 2.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGT  >  KT2440/1891896‑1892126
                                                                                                       |                                                                                                                               
acTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGCGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCCTATTACACATACACGGCCCACTGAAAAACCCGAAGTaaaaaat                                                                                                   >  5:177407/1‑133 (MQ=11)
acTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACAGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCAGCGCCATATTACACATACACGGGCCACAGAACACCCGGAAGTAAAaaag                                                                                                   >  5:177411/1‑134 (MQ=11)
     aTGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGGACACCCCGAAGTAAAAAAGAAAGCGa                                                                                            >  8:293258/1‑136 (MQ=17)
        tataGCGCAGCGGTAGTTCAGTCGGTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCGGCGCCCTATTACACTTGCAGGGCCCCCTGACCCCCCCGAACTAAAAAAGAAAGCCAccg                                                                                         >  7:234449/1‑134 (MQ=11)
            gcgcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCCGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGCCCTTAg                                                                                      >  4:188238/1‑135 (MQ=17)
            gcgcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCGGCGCCCTATTCCACTTCCAGGGCCCCCTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAgg                                                                                     >  7:84300/1‑136 (MQ=16)
            gcgcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTACCCTGCGCCCTATTCCACTTCCCGGGCCCACTGTCCACCCGGCAGTAAAAAAGAAAGCGCCCtt                                                                                        >  4:323521/1‑133 (MQ=11)
              gcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGGCATATTCCACTTCCCGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAggg                                                                                    >  3:227098/1‑135 (MQ=18)
              gcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGGGCCCGGCGCCATATTACACCTCCAGGGCCCCCGGACACCCCCGAAGTAAAAAAGAAAGCGCCCGTAggg                                                                                    >  1:265236/1‑135 (MQ=11)
               cAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCCCTGAACCCCCGGAAGTAAAAAAGAAAGCGACCTTAGGGt                                                                                   >  4:355007/1‑135 (MQ=21)
               cAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCGGCGCCATATTCCACTTCCAGGGCCCCCTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGt                                                                                   >  5:148800/1‑135 (MQ=21)
               cAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCAGGGGCATATTGCACGTCCAGGGCCCACTGGACACCCGGAAGTAAAAAAAAAAGCGCCCATAgggg                                                                                   >  1:15691/1‑134 (MQ=11)
                   ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTGGAGTCCCGTGCCGAGCGCCAAATTACCCTTCCCGGGCCCCCGGAACACCCCGTCGTAAAAAACAAAGCCCCCGTAgggggggg                                                                               >  6:75516/1‑130 (MQ=2)
                   ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGCGTCCCGTCCGCTGCGGCGTATTCCACTTCCAGTGGCCCCTGAACCCCCGGCACTAAAAAACAAAGCCACCATAGGGTCGcg                                                                               >  6:175289/1‑134 (MQ=11)
                   ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTGCGCGGCGCCCAATTCCACGTCCCGGGGCCCCGGCACACCCCGGAGTAAAAAAGAAATCCACCGTATGGtgggg                                                                               >  5:398572/1‑131 (MQ=2)
                   ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGGTGCGCCATATTCCACGTCCAGGGCCCCCTGCACACCCGGAAGTTAAAAAGAAAGCGCCCTGAGGGTCGcg                                                                               >  1:89985/1‑134 (MQ=11)
                   ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCCCTGCACACCCGGAAGTAAAAAAGAAAGCGCCCTTAGGGTCGcg                                                                               >  6:336808/1‑134 (MQ=255)
                   ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCCCTGAACCCCCGGAAGTAAAAAAGAAAGCGTCCGTAGGGGCGCtg                                                                              >  8:28749/1‑135 (MQ=11)
                   ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTACACTTCCAGGGCCCCCTGAACCCCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTAGcg                                                                               >  5:195065/1‑134 (MQ=255)
                   ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATAATCCACTTCCAGGGCCCCCGGAACACCCGGAAGTAAAAAAGAAAGCGACCGGAGGGTCGcgg                                                                              >  7:443321/1‑134 (MQ=11)
                   ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCGGCGCCATATTCCACTTCCAGGGCCCACTGCACACCCCGAAGTAAAAAAGAAATCGACCTTAGGGTCGcg                                                                               >  3:213679/1‑134 (MQ=255)
                   ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGGTCGAGTCCCGTCCGCTGCGGCATATTACACTTCCAGGGCCCACTGCGCACCCGGGAGTAAAAAAGAAATCGCCCTTAGGGTCGcg                                                                               >  6:155968/1‑134 (MQ=11)
                     tAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGGTCGAGTCCCGTGCCCTCCGCCCGATTACACCTACAGGGCCCACTGCACACCCGGAAGTAAAAAATAAATACACCTGTGGGTCGcggg                                                                             >  2:368502/1‑132 (MQ=11)
                           aGTGGGTTCGAATACCGGCCGGCCACGCCGGGGGGCGCGGGGGCGAGTCCCGGCCGCGGCGCCAAATTCCACTTCCCGGCCCCACTGAACACCCGGAAGTAAAAACGAAAGCGACCTTAGGGTCGCTTTTTTTa                                                                        <  2:89985/134‑1 (MQ=11)
                              gggTTACAATCCGGGCCGGGCACGCCGGGGGGGGGGGGTTCGAGTCCCGGCCGCTGCGCAATATTCCACTTCCCGGGCCCACTCAACACCCGCAATTAAAAAAGAACGCGACCTTAGGGTCGCTTTTTTTATTTCt                                                                   <  8:84300/135‑1 (MQ=255)
                                       tcccGGCGTGTCACGCGGGGGGTCGCGGGTGCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTTAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGcatcat                                                            <  6:160279/132‑1 (MQ=14)
                                                 tCAGGCGGGGGGTCGCGGGTTCGAGGCCCGGCCGCTGCGCCATATTCCACTTCCCGGGCCCACTGAACACCCGGAAGTAAAACAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTAc                                                <  7:293258/136‑1 (MQ=18)
                                                 tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCAATATTCCACTTCCAGGGCCCACTGAACACCCGCAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTAc                                                <  7:70754/136‑1 (MQ=255)
                                                                           cccgcccgCCGCCCCAAATCCCCCACCCCGGCCCCCCTCAACACCCGCACTAAAAAAACACCGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTgcg                       <  1:183434/135‑1 (MQ=255)
                                                                                       gcCATATCCCCCATCCCGGGCCCACTGAACACCCGTAAGTAAAAAAGACCGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTcc          <  8:244548/136‑1 (MQ=37)
                                                                                       gcCATAATCCACTGCCCGGCCCCACTGACCACCCGGAAGTAAAAAACAACGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTcc          <  7:106287/136‑1 (MQ=255)
                                                                                                cccTTCCAGGCCCCACAGACCACCCGCAACTAAAAACCACCGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGt  <  5:347271/133‑1 (MQ=255)
                                                                                                cACTTCCCGGGCCCACTCACCACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGt  <  4:395138/135‑1 (MQ=255)
                                                                                                cACTTCCAGGGCCCACTGAACACCCGCAAGTAAAAAACAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGt  <  6:148800/135‑1 (MQ=255)
                                                                                                cACTGCCCGGACCCACTCAACACCCGGAAGAAAAAACCAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGt  <  4:442553/135‑1 (MQ=255)
                                                                                                       |                                                                                                                               
ACTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGT  >  KT2440/1891896‑1892126

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: