Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,854,144 +G coding (891/1011 nt) hypE → carbamoyl dehydratase, hydrogenases 1,2,3 maturation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,854,1401.G87.1% 216.7 / 28.2 71coding (887/1011 nt)hypEcarbamoyl dehydratase, hydrogenases 1,2,3 maturation protein
Reads supporting (aligned to +/- strand):  ref base . (5/4);  new base G (32/29);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GACTGGACGCCCTGAACTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCAC  >  NC_000913/2854041‑2854235
                                                                                                    |                                                                                                
gACTGGACGCCCTGAACTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACt                                                                                                   <  1:1324628/100‑1 (MQ=255)
gACTGGACGCCCTGAACTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTACTGGCAGCGTTACATTCCCATCCACt                                                                                                   >  2:33129/1‑100 (MQ=255)
     gACGCTTTGAACTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTgggg                                                                                             <  2:452285/99‑1 (MQ=255)
     gACGCCCTGAACTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTggggg                                                                                            <  1:1264696/100‑2 (MQ=255)
        tccgatctcTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGgaaa                                                                                         <  2:1220135/92‑2 (MQ=255)
         cccTGAACTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGgaaag                                                                                        <  2:1643690/100‑3 (MQ=255)
         cccTGAACTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGgaaag                                                                                        >  1:437509/1‑98 (MQ=255)
         cccTGAACTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGgaaag                                                                                        >  1:437503/1‑98 (MQ=255)
         cacTGAACTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGgaaa                                                                                         <  2:679176/97‑2 (MQ=255)
          ccTGAACTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTAATGGCAGCGTTACATACCCATCCACTGGGGGAAAGa                                                                                        >  2:540989/1‑100 (MQ=255)
            tGAACTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGAcg                                                                                      <  1:1406022/100‑1 (MQ=255)
              aaacttGCCCACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTCGGGGAAAGAcgc                                                                                     <  1:1351424/95‑1 (MQ=255)
               aCTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACgcgg                                                                                   >  2:1490180/1‑100 (MQ=255)
                cTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACgcggc                                                                                  >  1:717549/1‑100 (MQ=255)
                  ttGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGgcgc                                                                                >  2:1205029/1‑100 (MQ=255)
                  tggcgaaCCAAGGCAATCTAGTAAGAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACgcggcg                                                                                 <  1:1434640/94‑1 (MQ=255)
                   tGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCt                                                                               <  2:1366607/100‑1 (MQ=255)
                   tGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCt                                                                               >  2:924372/1‑100 (MQ=255)
                   tGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGCAGCGGCGCTg                                                                              >  1:330058/1‑100 (MQ=255)
                    gCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTg                                                                              <  2:1458722/100‑1 (MQ=255)
                    gCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTg                                                                              >  2:1313768/1‑100 (MQ=255)
                     ccAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGACGCTGa                                                                             >  2:170944/1‑100 (MQ=255)
                      cAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGAt                                                                            >  1:104269/1‑100 (MQ=255)
                        aCGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGAAacggcgataa                                                                             >  2:287967/1‑93 (MQ=255)
                         cGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTgg                                                                         <  2:1683001/100‑1 (MQ=255)
                         cGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTgg                                                                         >  1:965484/1‑100 (MQ=255)
                           aaGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTg                                                                       <  1:1652597/100‑1 (MQ=255)
                            aGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGgaga                                                                      >  2:780982/1‑100 (MQ=255)
                            aGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGa                                                                      >  1:993975/1‑100 (MQ=255)
                                  aCTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAgtggt                                                                <  1:1743005/100‑1 (MQ=255)
                                  aCTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAgtggt                                                                >  2:1335436/1‑100 (MQ=255)
                                    tAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGtggtgg                                                              <  2:888184/100‑1 (MQ=255)
                                        aaTAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACg                                                          <  2:1426846/100‑1 (MQ=255)
                                        aaTAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACg                                                          <  1:387305/100‑1 (MQ=255)
                                        aaTAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGTGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACg                                                          <  1:1124817/100‑1 (MQ=255)
                                         aTAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGt                                                         <  2:864672/100‑1 (MQ=255)
                                         aTAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGt                                                         <  2:717553/100‑1 (MQ=255)
                                           aGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTaa                                                       >  1:232415/1‑100 (MQ=255)
                                            gCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTaaa                                                      >  2:1507569/1‑100 (MQ=255)
                                            gCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTaaa                                                      >  2:911170/1‑100 (MQ=255)
                                              tGTTGAATGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAACACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAgg                                                    <  1:975250/100‑1 (MQ=255)
                                              tGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAgg                                                    <  1:1490169/100‑1 (MQ=255)
                                               gTTGAACGCAACGCGGCAGAGCCAGTGCTGGGAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGgt                                                   <  1:1227492/100‑1 (MQ=255)
                                                   aacgcaacgcGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTc                                               >  1:1066280/1‑100 (MQ=255)
                                                    acgcaacgcGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCg                                              <  2:366068/100‑1 (MQ=255)
                                                       caacgcGGTAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCt                                           <  2:855732/100‑1 (MQ=255)
                                                       caacgcGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCt                                           >  2:861548/1‑100 (MQ=255)
                                                       caacgcGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCt                                           <  2:1524028/100‑1 (MQ=255)
                                                        aacgcGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCtt                                          <  1:1397167/100‑1 (MQ=255)
                                                        aacgcGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCtt                                          >  1:561624/1‑100 (MQ=255)
                                                            cgGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCg                                      <  1:1555055/100‑1 (MQ=255)
                                                             ggCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGt                                    >  1:1773628/1‑100 (MQ=255)
                                                              gCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCgg                                     <  2:395785/99‑1 (MQ=255)
                                                              gCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGt                                    >  2:1337599/1‑100 (MQ=255)
                                                                agagCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCt                                  >  2:495279/1‑100 (MQ=255)
                                                                agagCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCt                                  >  2:414760/1‑100 (MQ=255)
                                                                  agCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGt                                >  1:1663711/1‑100 (MQ=255)
                                                                      aGTGCTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATgg                            <  1:1574951/100‑1 (MQ=255)
                                                                          cTGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTg                        <  1:1756946/100‑1 (MQ=255)
                                                                           tGGCAGCGTTACATTCCCATCCACTGGGGGAAAGATGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTg                        <  2:595331/99‑1 (MQ=255)
                                                                           tGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGa                       >  1:943639/1‑100 (MQ=255)
                                                                           tGGCAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGa                       >  2:1458134/1‑100 (MQ=255)
                                                                              cAGCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAAc                    <  2:1410540/100‑1 (MQ=255)
                                                                                gCGTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGa                  <  1:1325291/100‑1 (MQ=255)
                                                                                  gTTACATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAAc                >  1:667160/1‑100 (MQ=255)
                                                                                     aCATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCt             >  1:948433/1‑100 (MQ=255)
                                                                                      cATTCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTc            >  1:1725638/1‑100 (MQ=255)
                                                                                         tCCCATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGAt         <  1:924366/100‑1 (MQ=255)
                                                                                           ccATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGAttt       <  2:388411/100‑1 (MQ=255)
                                                                                            cATCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTa      >  1:747768/1‑100 (MQ=255)
                                                                                             aTCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTAc     >  2:721396/1‑100 (MQ=255)
                                                                                             aTCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTAc     >  1:525552/1‑100 (MQ=255)
                                                                                              tCCACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTAcc    >  2:480292/1‑100 (MQ=255)
                                                                                               ccACTGGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTAcc    <  1:1445274/99‑1 (MQ=255)
                                                                                               ccACTGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACcac  >  1:445805/1‑100 (MQ=255)
                                                                                                    |                                                                                                
GACTGGACGCCCTGAACTTTGCCAACGAAGGCAAACTAGTAATAGCTGTTGAACGCAACGCGGCAGAGCAAGTGCTGGCAGCGTTACATTCCCATCCACTGGGGAAAGACGCGGCGCTGATTGGTGAAGTGGTGGAACGTAAAGGTGTTCGTCTTGCCGGTCTGTATGGCGTGAAACGAACCCTCGATTTACCAC  >  NC_000913/2854041‑2854235

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: