Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,905,536 C→A 9.7% P583P (CCC→CCA)  nirB → nitrite reductase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,905,5360CA 9.7% 116.9 / 2.2 41P583P (CCC→CCA) nirBnitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base C (22/15);  new base A (2/2);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

AACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGGTTTCGAAACC  >  NC_002947/1905390‑1905683
                                                                                                                                                  |                                                                                                                                                   
aaCATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGAc                                                                                                                                                  >  1:311474/1‑150 (MQ=255)
    ttCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCTCTGCAGGACACCAACGACACCTTCCTGGCCAACATCCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGc                                                                                                                                              <  2:680002/150‑1 (MQ=255)
       tcGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCtca                                                                                                                                           >  1:696083/1‑150 (MQ=255)
                gCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCtcatcgtcatcg                                                                                                                                  <  2:1142513/150‑1 (MQ=255)
                        cAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTg                                                                                                                          >  1:544511/1‑150 (MQ=255)
                             gATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCaa                                                                                                                     >  1:195389/1‑150 (MQ=255)
                                 aTGGACCCGGCGTTGGATCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGaaa                                                                                                                 >  1:369011/1‑150 (MQ=255)
                                    gACCCGGCGTTGGTGCCGCTGCAGGAAACCAACAACACCTACATGGCCAACAAGCAGAAAATCGACACCTATCCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATAc                                                                                                              <  1:841418/150‑1 (MQ=255)
                                       ccGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAGGCGCATCGTCATCGGCCAGGTGGCCAAGAAATACGAc                                                                                                           >  2:1066130/1‑150 (MQ=255)
                                       ccGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGAc                                                                                                           <  2:696083/150‑1 (MQ=255)
                                        cGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGAcc                                                                                                          >  1:106687/1‑150 (MQ=255)
                                                 tGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCa                                                                                                 >  1:591634/1‑150 (MQ=255)
                                                        gCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACAACCTGTACACCAAGATCAc                                                                                          <  2:1269394/150‑1 (MQ=255)
                                                        gCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACAACCTGTACACCAAGATCAc                                                                                          >  1:1269394/1‑150 (MQ=255)
                                                         cAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCAcc                                                                                         <  2:1182467/150‑1 (MQ=255)
                                                           ggACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGCAGCGGAAATAACCCCCGACAAGCTCATCGTCATCTACCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCgg                                                                                       >  1:288086/1‑150 (MQ=255)
                                                             acacCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTg                                                                                     <  2:544511/150‑1 (MQ=255)
                                                               acCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGc                                                                                   >  2:1229745/1‑150 (MQ=255)
                                                               acCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCAGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGa                                                                                              <  1:532720/139‑1 (MQ=255)
                                                               acCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCAGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGa                                                                                              >  2:532720/1‑139 (MQ=255)
                                                                          cTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCAGACAAGCtcatcgtc                                                                                                                                      <  2:828551/88‑1 (MQ=255)
                                                                          cTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCAGACAAGCtcatcgtc                                                                                                                                      >  1:828551/1‑88 (MQ=255)
                                                                                  ccAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCg                                                                >  1:637996/1‑150 (MQ=255)
                                                                                    aaCATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGgc                                                              >  2:178143/1‑150 (MQ=255)
                                                                                        tGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGccc                                                          >  2:1092247/1‑150 (MQ=255)
                                                                                                        cTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGcc                                          >  2:1189704/1‑150 (MQ=255)
                                                                                                         tATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCg                                         >  1:251519/1‑150 (MQ=255)
                                                                                                           ttCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCGCATCGTCATCGGCAAGGTGGGCAAGAAATAAGACATGTACACCAAGATCACAGGAGGCCAGCGCATCGAGCTGTTAGGCGCCCAGTTGAACGCGCTGACGcc                                       >  2:552739/1‑149 (MQ=255)
                                                                                                                   ggCCGCGCATTCCCGGTGGTTAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTggg                               <  1:552739/149‑1 (MQ=255)
                                                                                                                    gCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTgggg                              <  2:106687/150‑1 (MQ=255)
                                                                                                                    gCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTgggg                              <  2:591634/150‑1 (MQ=255)
                                                                                                                        cgcATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGCCCTGTAAATCAATAGCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGaa                          <  2:369011/150‑1 (MQ=255)
                                                                                                                         gcATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAAc                         <  1:1189704/150‑1 (MQ=255)
                                                                                                                            ttCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGa                      <  1:1092247/150‑1 (MQ=255)
                                                                                                                                cGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGa                  >  2:291851/1‑150 (MQ=255)
                                                                                                                                     gtgAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGc                                        >  1:825545/1‑123 (MQ=255)
                                                                                                                                     gtgAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGc                                        <  2:825545/123‑1 (MQ=255)
                                                                                                                                       gAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGcc                                                           >  2:486182/1‑102 (MQ=255)
                                                                                                                                       gAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGcc                                                           <  1:486182/102‑1 (MQ=255)
                                                                                                                                       gAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGGt           >  1:145659/1‑150 (MQ=255)
                                                                                                                                           tAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGGTTTCg       >  2:392913/1‑150 (MQ=255)
                                                                                                                                              cccccGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGGTTTCGaaa    >  2:363147/1‑150 (MQ=255)
                                                                                                                                                cccGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGGTTTCGAAAcc  >  1:400435/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAACTGATCGAGGCCGGTTTCGAAACC  >  NC_002947/1905390‑1905683

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑