Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,905,511 G→A 7.5% R575H (CGC→CAC)  nirB → nitrite reductase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,905,5110GA 7.5% 149.2 / 3.1 53R575H (CGC→CAC) nirBnitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base G (29/20);  new base A (2/2);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

TGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAAC  >  NC_002947/1905363‑1905660
                                                                                                                                                    |                                                                                                                                                     
tgtgACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCgcgc                                                                                                                                                      >  1:1092521/1‑150 (MQ=255)
 gtgACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCa                                                                                                                                                     <  1:151911/150‑1 (MQ=255)
        cTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCgg                                                                                                                                              >  1:903339/1‑150 (MQ=255)
         tGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCggt                                                                                                                                             >  2:1045618/1‑150 (MQ=255)
          gCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCggtg                                                                                                                                            <  2:903339/150‑1 (MQ=255)
            aaGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCggtggt                                                                                                                                          <  2:498100/150‑1 (MQ=255)
             aGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAATCAGCCGAGCACGGGCCGGGCGTTGGTGCCGCTGCTGGACACCCACGACACTTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCAGGGCGGTCCCGttgttg                                                                                                                                         >  1:934942/1‑150 (MQ=255)
                     gTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAAcc                                                                                                                                 <  1:717991/150‑1 (MQ=255)
                          cAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGa                                                                                                                            >  2:257345/1‑150 (MQ=255)
                               ttCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGCGGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGc                                                                                                                       >  1:836441/1‑150 (MQ=255)
                                     ccTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCtcatcg                                                                                                                 >  1:1161834/1‑150 (MQ=255)
                                              ggAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGGGAAAAAACCCCCGACAAGCTCATCGTCATCGGcc                                                                                                        >  1:1094779/1‑150 (MQ=255)
                                               gAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCa                                                                                                       >  1:565859/1‑150 (MQ=255)
                                                aaCCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAg                                                                                                      <  1:257345/150‑1 (MQ=255)
                                                 aCCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAgg                                                                                                     >  1:1294417/1‑150 (MQ=255)
                                                  ccAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCACATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGt                                                                                                    >  1:1521472/1‑150 (MQ=255)
                                                    aGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTgg                                                                                                  <  2:1294417/150‑1 (MQ=255)
                                                        gATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCaa                                                                                              >  2:658845/1‑150 (MQ=255)
                                                             tGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAAt                                                                                         <  2:1161834/150‑1 (MQ=255)
                                                               gACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATAc                                                                                       <  1:1045618/150‑1 (MQ=255)
                                                               gACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCACATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATAc                                                                                       <  2:1521472/150‑1 (MQ=255)
                                                                aCCCGGCGCTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTAGTGAAGCGAAGACCTGGTCGTGAAAGAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAAGACg                                                                                      >  1:1191036/1‑150 (MQ=255)
                                                                  ccGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCggtg                                                                                                                                            >  1:1539111/1‑94 (MQ=255)
                                                                  ccGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCggtg                                                                                                                                            <  2:1539111/94‑1 (MQ=255)
                                                                        ttGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGTAATACTACCTGTac                                                                              <  2:808425/150‑1 (MQ=255)
                                                                             gCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAAcccc                                                                                                                               <  1:450565/96‑1 (MQ=255)
                                                                             gCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAAcccc                                                                                                                               >  2:450565/1‑96 (MQ=255)
                                                                                    cAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAAccc                                                                                                                                <  2:1415681/88‑1 (MQ=255)
                                                                                    cAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAAccc                                                                                                                                >  1:1415681/1‑88 (MQ=255)
                                                                                       gACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCtcatcg                                                                                                                 <  1:59674/100‑1 (MQ=255)
                                                                                       gACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCtcatcg                                                                                                                 >  2:59674/1‑100 (MQ=255)
                                                                                       gACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCtcatc                                                                                                                  >  2:1322390/1‑99 (MQ=255)
                                                                                       gACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCtcatc                                                                                                                  <  1:1322390/99‑1 (MQ=255)
                                                                                       gACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCGCATCGTCATCGGCCAGGTGGCCAAGCAAAACGACCCGCACACCAAGATTCCCGGt                                                               >  2:36170/1‑150 (MQ=255)
                                                                                              aCGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCggtggt                                                                                                                                          >  2:330484/1‑68 (MQ=255)
                                                                                              aCGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCggtggt                                                                                                                                          <  1:330484/68‑1 (MQ=255)
                                                                                                    ccTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCg                                                  >  1:1229563/1‑150 (MQ=255)
                                                                                                    ccTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCg                                                  >  2:585424/1‑150 (MQ=255)
                                                                                                    ccTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCACATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCg                                                  >  1:518383/1‑150 (MQ=255)
                                                                                                     cTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGa                                                 >  1:925860/1‑150 (MQ=255)
                                                                                                      ttCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGAc                                                >  2:1565542/1‑150 (MQ=255)
                                                                                                       tCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGAcc                                               >  1:1431393/1‑150 (MQ=255)
                                                                                                               aaCATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGt                                                                                <  2:946004/109‑1 (MQ=255)
                                                                                                               aaCATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGt                                                                                >  1:946004/1‑109 (MQ=255)
                                                                                                                          aaaCGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCa                            <  1:585424/150‑1 (MQ=255)
                                                                                                                               gCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCg                                         <  2:719255/132‑1 (MQ=255)
                                                                                                                               gCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCg                                         >  1:719255/1‑132 (MQ=255)
                                                                                                                                   cTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGc                                                                                                                       <  2:1382034/50‑1 (MQ=255)
                                                                                                                                   cTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGc                                                                                                                       >  1:1382034/1‑50 (MQ=255)
                                                                                                                                   cTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGcc                   >  1:95879/1‑150 (MQ=255)
                                                                                                                                        cGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGa              <  2:934942/150‑1 (MQ=255)
                                                                                                                                         ggtggtGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGAt             >  2:830278/1‑150 (MQ=255)
                                                                                                                                              tgCCGCACATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTggg        <  2:518383/150‑1 (MQ=255)
                                                                                                                                               gCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTgggg       >  1:181168/1‑150 (MQ=255)
                                                                                                                                                  gcgcATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCt                                                                                  <  1:70826/72‑1 (MQ=255)
                                                                                                                                                  gcgcATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCt                                                                                  >  2:70826/1‑72 (MQ=255)
                                                                                                                                                  gcgcATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGa    >  1:34887/1‑150 (MQ=255)
                                                                                                                                                   cgcATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGaa   <  1:830278/150‑1 (MQ=255)
                                                                                                                                                   cgcATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGaa   <  2:1431393/150‑1 (MQ=255)
                                                                                                                                                    gcATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAAc  <  2:1092521/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGTGACATCTGCAAGCCGGCGGTGGGCAACATTCTCGCCTCGTGCTGGAACCAGCCGATCATGGACCCGGCGTTGGTGCCGCTGCAGGACACCAACGACACCTTCATGGCCAACATGCAGAAAAACGGCACCTATTCGGTGGTGCCGCGCATTCCCGGTGGTGAAATAACCCCCGACAAGCTCATCGTCATCGGCCAGGTGGCCAAGAAATACGACCTGTACACCAAGATCACCGGTGGCCAGCGCATCGACCTGTTCGGCGCCCAGTTGCACGAGCTGCCGCTGATCTGGGGCGAAC  >  NC_002947/1905363‑1905660

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: