Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,903,708 C→A 12.0% intergenic (+233/‑38) 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,7080CA12.0% 121.1 / 4.8 46intergenic (+233/‑38)PP_1703/nirBassimilatory nitrate reductase/sulfite reductase/nitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base C (18/22);  new base A (4/2);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATGGTCGGCCACCATTGCGTCGAACAACTGGTCAGTCGCGGCGCCCTGCAGC  >  NC_002947/1903560‑1903854
                                                                                                                                                    |                                                                                                                                                   
gAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCATCCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTc                                                                                                                                                     <  1:354987/149‑1 (MQ=255)
gAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTc                                                                                                                                                     <  1:147580/149‑1 (MQ=255)
       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                                                                       >  2:102057/1‑135 (MQ=255)
                                                                                           aCGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTgg                                                                       <  1:102057/135‑1 (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                                                   <  2:333828/119‑1 (MQ=255)
                                                                                                                               aGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATg                                                   >  1:333828/1‑119 (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                                                    <  5:350216/109‑1 (MQ=38)
                                                                                                                                         gcCTGGAGCCGAAGCTCCAGGCGCCTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCAt                                                    >  6:350216/1‑109 (MQ=38)
                                                                                                                                                 cTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATGGTCGGCCACCATTGCGTCGAACAACTGGTCAGTCGCGGCGCCCTGCAg   >  5:371362/1‑149 (MQ=255)
                                                                                                                                                  ttCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATGGTCGGCCACCATTGCGTCGAACAACTGGTCAGTCGCGGCGCCCTGCAgc  >  2:208926/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAAGAGGCCGACACAGGCAACCCAGGCAACCCACGCAACCCGATTGGCACGCTCCCTGCTTAACCTCGGATACAACAAAGCCCACTGATCAACGGCGATCAACCGGGCCCCCTACCGTGATGAATACAGGCAAAGGCGCCTGGAGCTTCGGCTCCAGGCG‑CTTTTTTTTTGATCGAGGATCGGCTTATGAAGGCAACAGCGAGCAGCGGGCAACGAGAGCGACTGGTCATCGTCGGCAACGGCATGGTCGGCCACCATTGCGTCGAACAACTGGTCAGTCGCGGCGCCCTGCAGC  >  NC_002947/1903560‑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: ‑