Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,496,373 A→C E788A (GAA→GCA)  nirB → nitrite reductase, large subunit, NAD(P)H‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,496,3730AC98.2% 200.4 / ‑4.1 57E788A (GAA→GCA) nirBnitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (21/35);  total (22/35)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CCGTACTGCCGACAAACTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAGCGCGAACAGCACCGTCCGGCAACGCCG  >  NC_000913/3496233‑3496512
                                                                                                                                            |                                                                                                                                           
ccGTACTGCCGACAAACTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTAc                                                                                                                                     <  2:620037/149‑1 (MQ=255)
ccGTACTGCCGACAAACTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACtg                                                                                                                                   <  1:793105/151‑1 (MQ=255)
ccGTACTGCCGACAAACTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACt                                                                                                                                    <  1:418077/150‑1 (MQ=255)
        ccGACAAACTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTg                                                                                                                           <  2:616098/151‑1 (MQ=255)
            cAAACTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACt                                                                                                                       <  2:232155/151‑1 (MQ=255)
             aaaCTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTg                                                                                                                      <  2:283786/151‑1 (MQ=255)
               aCTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGa                                                                                                                     >  2:455491/1‑150 (MQ=255)
                  gACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGaaa                                                                                                                   <  2:734070/149‑1 (MQ=255)
                    cgcgTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGt                                                                                                               <  1:778343/151‑1 (MQ=255)
                     gcgTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTc                                                                                                              <  1:498343/151‑1 (MQ=255)
                         aCCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATa                                                                                                          <  1:540812/151‑1 (MQ=255)
                           cGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACg                                                                                                        >  1:43419/1‑151 (MQ=255)
                               ccGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGt                                                                                                    <  1:374339/151‑1 (MQ=255)
                                      tAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTgcgc                                                                                              <  2:68695/150‑1 (MQ=255)
                                           aaCCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTc                                                                                        <  1:265826/151‑1 (MQ=255)
                                            aCCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCg                                                                                       <  2:683466/151‑1 (MQ=255)
                                            aCCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCg                                                                                       <  1:588619/151‑1 (MQ=255)
                                             ccTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGc                                                                                      >  1:142752/1‑151 (MQ=255)
                                             ccTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGc                                                                                      <  2:593030/151‑1 (MQ=255)
                                              cTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCt                                                                                     >  1:855060/1‑151 (MQ=255)
                                               tCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCtt                                                                                    >  2:293645/1‑151 (MQ=255)
                                                cGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCt                                                                                     <  2:420135/149‑1 (MQ=255)
                                                cGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTc                                                                                   >  2:712915/1‑151 (MQ=255)
                                                 gAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTc                                                                                   >  2:82871/1‑150 (MQ=255)
                                                          aTCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTtc                                                                          >  1:836506/1‑150 (MQ=255)
                                                           tCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTtca                                                                         >  2:360932/1‑150 (MQ=255)
                                                            cGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTtcat                                                                        <  1:213956/150‑1 (MQ=255)
                                                             gATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTtcatca                                                                      <  1:384858/151‑1 (MQ=255)
                                                             gATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTtcatca                                                                      <  2:768848/151‑1 (MQ=255)
                                                                tACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAAc                                                                    >  2:915951/1‑150 (MQ=255)
                                                                tACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAAc                                                                    >  1:142792/1‑150 (MQ=255)
                                                                tACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAAc                                                                    >  2:369644/1‑150 (MQ=255)
                                                                         gCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGc                                                          <  1:800973/151‑1 (MQ=255)
                                                                         gCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGc                                                          >  2:862123/1‑151 (MQ=255)
                                                                          cAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCg                                                         <  1:649246/151‑1 (MQ=255)
                                                                          cAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCg                                                         <  1:196469/151‑1 (MQ=255)
                                                                               aTCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGAcc                                                    >  2:922568/1‑151 (MQ=255)
                                                                                      acgacAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGt                                             <  2:842277/151‑1 (MQ=255)
                                                                                      acgacAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGt                                             >  1:910290/1‑151 (MQ=255)
                                                                                            aaCTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGc                                           >  2:679263/1‑147 (MQ=255)
                                                                                                ggggCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATgg                                     <  2:725867/149‑1 (MQ=255)
                                                                                                 gggCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGc                                  <  1:315648/151‑1 (MQ=255)
                                                                                                  ggCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGcc                                 <  2:43419/151‑1 (MQ=255)
                                                                                                   gCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCa                                <  1:313149/151‑1 (MQ=255)
                                                                                                    cTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGcc                                 >  1:370395/1‑149 (MQ=255)
                                                                                                       aaCGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAgc                            <  1:270927/151‑1 (MQ=255)
                                                                                                         cGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAgcgc                          <  1:424166/151‑1 (MQ=255)
                                                                                                          gCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAGcgcg                         <  2:112522/151‑1 (MQ=255)
                                                                                                           cacaTCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAGCGCGa                        <  2:414052/151‑1 (MQ=255)
                                                                                                             caTCTGGAAGAAGAGATGGCGCGCCTGCGTGCAGCGGAACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAgcgc                          >  2:761413/1‑147 (MQ=255)
                                                                                                                    aagaagAGATGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAGCGCGAACAGCACCg               >  2:750132/1‑151 (MQ=255)
                                                                                                                          agaTGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAGCGCGAACAGCACCGTCCGGc         >  2:298822/1‑151 (MQ=255)
                                                                                                                            aTGGCGCGCCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAGCGCGAACAGCACCGTCCGGCaa       <  1:854224/151‑1 (MQ=255)
                                                                                                                               gcgcgcCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAGCGCGAACAGCACCGTCCGGCAACGc    >  1:263769/1‑151 (MQ=255)
                                                                                                                                cgcgcCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGCACGTGCAGATGGTGCCAGAGCGCGAACAGCACCGTCCGGCAACGcc   <  1:574210/151‑1 (MQ=255)
                                                                                                                                cgcgcCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAGCGCGAACAGCACCGTCCGGCAACGcc   <  2:140949/151‑1 (MQ=255)
                                                                                                                                 gcgcCTGCGTGCAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAGCGCGAACAGCACCGTCCGGCAACGCCg  >  2:360654/1‑151 (MQ=255)
                                                                                                                                            |                                                                                                                                           
CCGTACTGCCGACAAACTGACGCGTACCGCACCGTGGTTAGAAAACCTCGAAGGCGGCATCGATTACCTGAAAGCAGTGATCATTGACGACAAACTGGGGCTGAACGCACATCTGGAAGAAGAGATGGCGCGCCTGCGTGAAGCGGTACTGTGTGAGTGGACTGAAACGGTCAATACGCCGTCTGCGCAGACTCGCTTCAAACACTTCATCAACAGCGACAAGCGTGACCCGAACGTGCAGATGGTGCCAGAGCGCGAACAGCACCGTCCGGCAACGCCG  >  NC_000913/3496233‑3496512

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑