Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,903,709 G→A 10.8% intergenic (+234/‑37) PP_1703 → / → nirB assimilatory nitrate reductase/sulfite reductase/nitrite reductase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,903,7090GA10.8% 118.7 / 5.0 46intergenic (+234/‑37)PP_1703/nirBassimilatory nitrate reductase/sulfite reductase/nitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base G (19/22);  new base A (3/2);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.57e-01

CCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATGGTCGGCCACCATTGCGTCGAACAACTGGTCAGTCGCGGCGCCCTGCAGC  >  NC_002947/1903567‑1903854
                                                                                                                                              |                                                                                                                                                  
ccGACACAGTCAACCCAGGCAACCCACGCAACCCGATTGTCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCCAAGGCGCCTGGAGCTTCGGCTCCa                                                                                                                                              <  6:300/149‑1 (MQ=255)
 cGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAg                                                                                                                                             <  8:337973/149‑1 (MQ=255)
 cGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCATAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGTCAAAGGCGCCTGGAGCTACGGCTCCat                                                                                                                                             >  7:362322/1‑148 (MQ=255)
   acacAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGgc                                                                                                                                           <  1:20103/149‑1 (MQ=255)
       tggcaacccaggcaacccaCGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑Ctt                                                                                                                                      >  4:268873/2‑149 (MQ=255)
       aggcaacccaggcaacccACGCTACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGGTCCAGGCG‑Ctt                                                                                                                                      >  7:244081/1‑149 (MQ=255)
       aggcaacccaggcaacccACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCACgcg‑cgt                                                                                                                                      >  1:336420/1‑147 (MQ=255)
            acccaggcaacccaCGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCCGAAGCTCCAGgc                                                                                                                                           <  2:297491/140‑1 (MQ=255)
            acccaggcaacccaCGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGACCCGAAGCTCCAGgc                                                                                                                                           >  1:297491/1‑140 (MQ=255)
             cccaggcaacccaCGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCACCCGGGCCCCCTACCGTGATGAATACCGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑Ctttttttt                                                                                                                                <  8:292777/149‑1 (MQ=255)
              ccccccaaccccCCCCACCCTTTTGCCCCGCACCCTTCTTAACCTCCTAAAACAAAAACCCCACTTATCACCGGTTTTCTTTCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑Cttttttttt                                                                                                                               <  1:316065/144‑1 (MQ=255)
                aggcaacccaCGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACGACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGa                                                                                                                             <  7:281433/149‑1 (MQ=255)
                   caacccaCGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCg                                                                                                                          <  1:32809/149‑1 (MQ=255)
                   caacccaCGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCg                                                                                                                          <  1:94627/149‑1 (MQ=255)
                          cGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCg                                                                                                                   <  5:369203/149‑1 (MQ=255)
                                  gATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGa                                                                                                           <  3:268873/149‑1 (MQ=255)
                                    ttGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAg                                                                                                         <  7:221568/149‑1 (MQ=255)
                                                    ttAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCg                                                                                         >  1:287635/1‑149 (MQ=255)
                                                              aTACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACgaga                                                                               <  2:287635/149‑1 (MQ=255)
                                                                  aacaaAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGa                                                                           >  1:272331/1‑149 (MQ=255)
                                                                     aaaGCCCACTGATCAACGGCTATCAACCGGGCCCCCTACCGTGTTTAATAAAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTg                                                                        <  8:244081/149‑1 (MQ=255)
                                                                     aaaGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTg                                                                        <  3:32827/149‑1 (MQ=255)
                                                                     aaaGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTg                                                                        >  1:342304/1‑149 (MQ=255)
                                                                           cACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCAtc                                                                  <  2:342304/149‑1 (MQ=255)
                                                                           cACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGCATCGGCGTATGAAGGCAAAAGCGAGCAGCGGGCAACGAGAGCGACGGGTCAtc                                                                  >  6:308523/1‑149 (MQ=255)
                                                                                    accgCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGACCTTCGGCTCCACCCG‑CATTTTTTTTTATACAGGATAGGCGTAATAAAGGAACAGCGCGCCGCCGGCAACCACAGCCCCTGGGCATCGTAGGCaaa                                                         >  4:186331/4‑148 (MQ=255)
                                                                                    aCGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTgg                                                                       <  1:102057/135‑1 (MQ=255)
                                                                                    aCGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTgg                                                                       >  2:102057/1‑135 (MQ=255)
                                                                                    aCGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATAGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAAc                                                         >  3:296747/1‑149 (MQ=255)
                                                                                     cGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACg                                                        <  1:136549/149‑1 (MQ=255)
                                                                                     cGGCGATCAACCGGGCCCCCTACCGTGATGAATAAAGGCAAAGGCGCCTGGAGCTTCGTCTCCAGGCG‑CTTTTTTTTTGATAGAGGATAGGCTTATGAAAGCAACAGCCAGCAGCGGGCCACGAGAGCACCTCTTCATCGTGGGCaaag                                                        >  2:350205/1‑147 (MQ=255)
                                                                                             aaCCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATAGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATGGTc                                                >  2:355452/1‑149 (MQ=255)
                                                                                               ccGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATGGTCgg                                              >  6:73066/1‑149 (MQ=255)
                                                                                                       ccTTCCGTGATGAATACAGGCAAAGGGGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATGGTCGGCCACCAtt                                      <  3:186331/149‑1 (MQ=255)
                                                                                                                     tACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATGGTCGGCCACCATTGCGTCGAACAACTg                        <  2:272331/149‑1 (MQ=255)
                                                                                                                        aGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATg                                                   >  1:333828/1‑119 (MQ=255)
                                                                                                                        aGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATg                                                   <  2:333828/119‑1 (MQ=255)
                                                                                                                        aGGCAAAGGCGCCTGGAGCTTAGGCTCCATGCG‑CTTTTTTTTTGATCCAGGATCGGCTTAAGAAGGAAACCGCGAGCAGCGGGCAACGACCGCCACTGGTCCTCGGCCGCAAAGGCAAGGTCGGCCAACCTTTCGTCGcacaacag                        >  8:152039/1‑144 (MQ=255)
                                                                                                                          gCAAAGGCGCCTGGATCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTgg                                                                       <  7:88278/97‑1 (MQ=255)
                                                                                                                               ggCGCGTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATGGTCGGCCACCATTGCGTCGAACAACTGGTCAGTcgcg              <  1:367563/149‑1 (MQ=255)
                                                                                                                                gcgcCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGCCTGGTCATCGTCGGCAACGGCATGGTCGGCCCCCATTGCGTCGAACAACTGGTCAGTCgcgg             >  7:230689/1‑149 (MQ=255)
                                                                                                                                 cgcCTGGAGCCGAAGCTCCAGGCG‑CTTTTTTTTTGATAGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCAACCTCGTCAACGCCATGGTCGGCCAACATTGAGACGAACAACTGGTCAGTAgcgga            >  8:249452/1‑148 (MQ=255)
                                                                                                                                  gcCTGGAGCCGAAGCTCCAGGCGCCTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCAt                                                    >  6:350216/1‑109 (MQ=38)
                                                                                                                                  gcCTGGAGCCGAAGCTCCAGGCGCCTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCAt                                                    <  5:350216/109‑1 (MQ=38)
                                                                                                                                          cTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATGGTCGGCCACCATTGCGTCGAACAACTGGTCAGTCGCGGCGCCCTGCAg   >  5:371362/1‑149 (MQ=255)
                                                                                                                                           ttCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATGGTCGGCCACCATTGCGTCGAACAACTGGTCAGTCGCGGCGCCCTGCAgc  >  2:208926/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATGGTCGGCCACCATTGCGTCGAACAACTGGTCAGTCGCGGCGCCCTGCAGC  >  NC_002947/1903567‑1903854

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