Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,495,452 T→C V481A (GTT→GCT)  nirB → nitrite reductase catalytic subunit NirB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,495,4520TC96.4% 186.0 / ‑2.3 56V481A (GTT→GCT) nirBnitrite reductase catalytic subunit NirB
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (27/27);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGTTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACC  >  NC_000913/3495321‑3495587
                                                                                                                                   |                                                                                                                                       
tgcCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAATCAGGGCATCGAAGCTaacaaca                                                                                                                                 >  1:776139/1‑140 (MQ=255)
tgcCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTaacaaca                                                                                                                                 >  1:553665/1‑140 (MQ=255)
  ccATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTaacaacaac                                                                                                                               >  2:424356/1‑140 (MQ=255)
         caaAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCg                                                                                                                        <  1:234397/140‑1 (MQ=255)
         caaAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCg                                                                                                                        <  1:275315/140‑1 (MQ=255)
           aaGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGaa                                                                                                                      <  2:553665/140‑1 (MQ=255)
             ggCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAaca                                                                                                                    <  1:855340/140‑1 (MQ=255)
             ggCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAaca                                                                                                                    <  1:417182/140‑1 (MQ=255)
                tGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACtt                                                                                                                 >  1:1013012/1‑140 (MQ=255)
                 gCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAac                                                                                                                     >  1:700214/1‑135 (MQ=255)
                 gCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAac                                                                                                                     <  2:700214/135‑1 (MQ=255)
                    acacaGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCt                                                                                                             >  2:113602/1‑140 (MQ=255)
                      acaGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTa                                                                                                           <  1:446086/140‑1 (MQ=255)
                        aGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTAtt                                                                                                         >  1:791874/1‑140 (MQ=255)
                         gTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTc                                                                                                        <  2:88398/140‑1 (MQ=255)
                            gcggcgCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTcgcg                                                                                                     >  1:21978/1‑140 (MQ=255)
                              ggcgCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAATCAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTc                                                                                                   <  2:776139/140‑1 (MQ=255)
                              ggcgCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTc                                                                                                   <  2:448519/140‑1 (MQ=255)
                                  ctgaaagctgaaaCCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGa                                                                                               <  2:472514/140‑1 (MQ=255)
                                   tgaaagctgaaaCCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGaa                                                                                              <  1:884715/140‑1 (MQ=255)
                                       agctgaaaCCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGt                                                                                          <  1:113602/140‑1 (MQ=255)
                                                 aaaGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGAt                                                                                <  2:708416/140‑1 (MQ=255)
                                                   aGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGTTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATcc                                                                              >  1:339392/1‑140 (MQ=255)
                                                   aGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGaa                                                                                              >  2:1196812/1‑124 (MQ=255)
                                                   aGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGaa                                                                                              <  1:1196812/124‑1 (MQ=255)
                                                         tACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTg                                                                        >  1:973332/1‑140 (MQ=255)
                                                          aCTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGa                                                                       >  1:1117953/1‑140 (MQ=255)
                                                            tGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAg                                                                     <  1:684924/140‑1 (MQ=255)
                                                            tGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAg                                                                     <  2:973332/140‑1 (MQ=255)
                                                              gCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGc                                                                   >  1:698872/1‑140 (MQ=255)
                                                                tGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCAt                                                                 >  1:198757/1‑140 (MQ=255)
                                                                  cGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTaa                                                              >  1:340079/1‑140 (MQ=255)
                                                                     tGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCAtt                                                                >  1:591036/1‑136 (MQ=255)
                                                                     tGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCAtt                                                                <  2:591036/136‑1 (MQ=255)
                                                                                 gCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAAc                                                >  1:1190129/1‑140 (MQ=255)
                                                                                 gCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAAc                                                <  1:455838/140‑1 (MQ=255)
                                                                                     gTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAActgct                                            >  2:936634/1‑140 (MQ=255)
                                                                                             ggTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAac                                    <  2:791874/140‑1 (MQ=255)
                                                                                                   gAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCa                              <  1:676015/140‑1 (MQ=255)
                                                                                                   gAACGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCa                              <  2:198757/140‑1 (MQ=255)
                                                                                                    aaCGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCaa                             >  1:301351/1‑140 (MQ=255)
                                                                                                    aaCGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCaa                             >  1:980230/1‑140 (MQ=255)
                                                                                                    aaCGCGGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCaa                             >  2:194628/1‑140 (MQ=255)
                                                                                                       gcgGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAgg                          >  1:471225/1‑140 (MQ=255)
                                                                                                       gcgGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAgg                          <  1:1181315/140‑1 (MQ=255)
                                                                                                        cgGAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGc                         >  2:880936/1‑140 (MQ=255)
                                                                                                         ggAACTGGCGAAACAGGGCATCGAAGTTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCt                        >  2:182972/1‑140 (MQ=255)
                                                                                                         ggAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCa                              >  2:197437/1‑134 (MQ=255)
                                                                                                         ggAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCa                              <  1:197437/134‑1 (MQ=255)
                                                                                                         ggAACTGGCGAAACAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCt                        >  2:879183/1‑140 (MQ=255)
                                                                                                                   aaaCAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGa              <  2:698872/140‑1 (MQ=255)
                                                                                                                     aCAGGGCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAg            >  2:139958/1‑140 (MQ=255)
                                                                                                                      cAGGGCATCGAAGCTACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGtt          <  2:340079/140‑1 (MQ=255)
                                                                                                                          gCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGt       >  1:194965/1‑140 (MQ=255)
                                                                                                                          gCATCGAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCGACGGTTGTGAAGTTTGt       >  2:1083063/1‑140 (MQ=255)
                                                                                                                               gAAGCTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAaacc  <  1:936634/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                       
TGCCATCAACAAAGGCTGCCACACAGTTGCGGCGCTGAAAGCTGAAACCAAAGCGGGTACTGGCTGCGGTGGCTGTATCCCGCTGGTCACTCAGGTACTGAACGCGGAACTGGCGAAACAGGGCATCGAAGTTAACAACAACCTGTGCGAACACTTTGCTTATTCGCGTCAGGAACTGTTCCATTTGATCCGCGTTGAAGGCATTAAAACCTTCGAAGAACTGCTGGCGAAACACGGCAAAGGCTACGGTTGTGAAGTTTGTAAACC  >  NC_000913/3495321‑3495587

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: