Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,898,986 G→A 11.1% intergenic (+56/‑464) PP_t18 → / → PP_1703 tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,898,9860GA11.1% 121.2 / 5.4 46intergenic (+56/‑464)PP_t18/PP_1703tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base G (21/19);  new base A (2/3);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

GATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTG  >  NC_002947/1898844‑1899127
                                                                                                                                              |                                                                                                                                             
gATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCAAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTGAGGGTCGc                                                                                                                                         >  4:261084/1‑149 (MQ=21)
    tataGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCCTAAGGTCGCttt                                                                                                                                      <  8:165572/148‑1 (MQ=21)
    tataGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGCTCGCCGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTACCAGGGCCCACTGAACACCAGGAAGTAAAAAAGAAAGCGACCTAACGGTCGCtttt                                                                                                                                     >  8:107226/1‑149 (MQ=11)
      taGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCtttttt                                                                                                                                   >  5:276146/1‑149 (MQ=34)
        gcgcAGCGGTTGTTCTGGCGGTTACACTACCGCCCTGTCCCGCCGGGGGTCGCGGGTTCGAGTCCCGGCCGCTCCGCAATATTCCACTTCCAGGGCCCACTCAACACCCGTAATAAAAAAAAAACGCGACTCTAAGGTCGctttcttt                                                                                                                                  <  2:354028/148‑1 (MQ=2)
        gcgcAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCttttttt                                                                                                                                  <  6:213021/148‑1 (MQ=34)
             gCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCt                                                                                                                            <  3:230749/149‑1 (MQ=39)
               ggTATTTCAGTCTTTTAGAATACGGGCCTGTCTCGGGGGGGGTCGCGTGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAATAAAA‑‑‑AAAGCGACCCTAAGGTCGCTTTCTTTTTTACttccgg                                                                                                                       <  4:219085/149‑6 (MQ=2)
                             ttAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCCCTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCggt                                                                                                            >  2:308842/1‑147 (MQ=255)
                               aGAATACTGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGttt                                                                                                           <  1:199201/148‑1 (MQ=255)
                               aGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGttt                                                                                                           <  7:107226/148‑1 (MQ=255)
                                gacccccGGCGTGACAGGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTAc                                                                                                         >  2:274628/6‑149 (MQ=255)
                                gAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACTCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTa                                                                                                          <  7:312749/148‑1 (MQ=255)
                                               aCGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAggg                                                                                          >  3:275321/1‑149 (MQ=255)
                                                    gggggTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAg                                                                                     <  4:339492/149‑1 (MQ=255)
                                                         tCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCGTTTTTTATTTATGCATCATGGAAAAGTTTAAGCCCAGCAAAAAAGGGCTAgattgg                                                                                >  1:179547/1‑144 (MQ=17)
                                                           gcgGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGcc                                                                              >  4:243146/1‑149 (MQ=255)
                                                            cgGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTTTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCa                                                                             <  1:274628/149‑1 (MQ=255)
                                                               gTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCActc                                                                          <  2:300947/149‑1 (MQ=255)
                                                                      tCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTcc                                                                   <  6:276146/149‑1 (MQ=255)
                                                                               cTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAATGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTc                                                          <  1:305633/149‑1 (MQ=255)
                                                                                          tCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCaaa                                               <  3:243146/149‑1 (MQ=255)
                                                                                          tCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCaaa                                               >  3:335000/1‑149 (MQ=255)
                                                                                          tCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCaaa                                               >  7:355824/1‑149 (MQ=255)
                                                                                            cACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAAcc                                             >  5:182568/1‑149 (MQ=255)
                                                                                               ttCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCt                                                                                                                            <  1:211298/67‑1 (MQ=255)
                                                                                               ttCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCt                                                                                                                            >  2:211298/1‑67 (MQ=255)
                                                                                                     ggCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATAAAGCAATCGTTTACGCCCTGCAACAAGGGGCTAGGTGCGCCCCACTTCTTCCCGCAGCGTCAGCCCCTCAAACCCCTCCTCGa                                    >  3:143112/1‑149 (MQ=37)
                                                                                                         cACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCAc                                >  5:273518/1‑149 (MQ=255)
                                                                                                              aaCACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAAccc                                            <  7:299484/132‑1 (MQ=255)
                                                                                                              aaCACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTcc                                        >  5:166625/1‑136 (MQ=255)
                                                                                                              aaCACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTcc                                        <  6:166625/136‑1 (MQ=255)
                                                                                                                 acCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTAAGCCCTGCAACAACGGGCTAGCTGCGCCACACTTATTCCCGCAGCGTCAGCACCTCAAACACGTCCTCGATCACCGGCAtat                        >  8:39305/1‑147 (MQ=255)
                                                                                                                   ccGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTcc                      >  7:344898/1‑149 (MQ=255)
                                                                                                                    cGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCt                     <  4:335000/149‑1 (MQ=255)
                                                                                                                     ggAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCg                              <  1:273062/139‑1 (MQ=255)
                                                                                                                     ggAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCg                              >  2:273062/1‑139 (MQ=255)
                                                                                                                        agaaataaaaaaagCGACCCTAAGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGAAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGa                 >  6:175323/7‑149 (MQ=255)
                                                                                                                         gaaataaaaaaagCGACCCTAAGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGa                 >  6:80645/6‑148 (MQ=255)
                                                                                                                            aaaaaGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATg             >  2:93495/1‑149 (MQ=255)
                                                                                                                             aaaaGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATgg            >  1:226642/1‑149 (MQ=255)
                                                                                                                               aagaaagCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGa          <  2:179547/149‑1 (MQ=255)
                                                                                                                                    agCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAggg       >  7:92770/1‑147 (MQ=255)
                                                                                                                                       gACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAggg       >  1:56204/1‑144 (MQ=255)
                                                                                                                                       gACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAggg       <  2:56204/144‑1 (MQ=255)
                                                                                                                                       gACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTg  <  7:39305/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTG  >  NC_002947/1898844‑1899127

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