Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,005,399 +C coding (1044/1092 nt) ygeW → putative carbamoyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,005,3981.C99.2% 332.5 / ‑3.6 129coding (1043/1092 nt)ygeWputative carbamoyl transferase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (45/83);  total (46/83)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 153 total aligned reads displayed.

AAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGA  >  NC_000913/3005323‑3005468
                                                                             |                                                                      
aaGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCgatgat                                                                         >  1:1340066/1‑76 (MQ=255)
 aGGCTTATACGAAAAATACCGTATCTCTACCTACAAAGAAGCCATCTGGAAGCCTTTTATCATCGCCGCGATGATc                                                                       <  2:4327361/76‑1 (MQ=255)
 aGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATc                                                                       <  2:1802603/76‑1 (MQ=255)
   gCGTATTCGAAAAATACCGTATCGCTACCTACTAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATcct                                                                     <  1:1041066/76‑3 (MQ=255)
   gCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATcct                                                                     >  1:3848111/1‑74 (MQ=255)
   gCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATcct                                                                     >  2:1376419/1‑74 (MQ=255)
   gCGTATTCGAAAAATACCGGATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATcct                                                                     >  1:4534491/1‑74 (MQ=255)
    cGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATcctg                                                                    <  2:2467424/76‑4 (MQ=255)
    cGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATcctg                                                                    <  1:823201/76‑4 (MQ=255)
    cGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATcctg                                                                    >  1:2027561/1‑73 (MQ=255)
       aTTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTcc                                                                  <  1:2619904/76‑1 (MQ=255)
       aTTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTcc                                                                  <  1:4384392/76‑1 (MQ=255)
       aTTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTcc                                                                  >  2:751009/1‑76 (MQ=255)
         tCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCg                                                                >  1:2734550/1‑76 (MQ=255)
          cGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGt                                                               >  2:4039466/1‑76 (MQ=255)
          cGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGt                                                               >  1:290979/1‑76 (MQ=255)
           gAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTa                                                              >  1:1507374/1‑76 (MQ=255)
           gAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTa                                                              >  1:946437/1‑76 (MQ=255)
            aaaaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTaa                                                             >  1:4010662/1‑76 (MQ=255)
            aaaaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTaa                                                             <  1:266748/76‑1 (MQ=255)
            aaaaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTaa                                                             <  2:3153185/76‑1 (MQ=255)
             aaaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTaaa                                                            <  1:172280/76‑1 (MQ=255)
             aaaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTaaa                                                            <  2:4069914/76‑1 (MQ=255)
             aaaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTaaa                                                            <  1:656253/76‑1 (MQ=255)
              aaaTCCCGTATCGCTGCCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAAt                                                           <  2:4691903/76‑1 (MQ=255)
              aaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAAt                                                           <  2:3143730/76‑1 (MQ=255)
              aaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAAt                                                           >  1:1384799/1‑76 (MQ=255)
                aTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATAc                                                         <  2:2304119/76‑1 (MQ=255)
                 tACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACg                                                        <  1:2021042/76‑1 (MQ=255)
                 acccGTACCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACg                                                        <  2:696828/74‑1 (MQ=255)
                  aCCGTATCGCTACGTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGc                                                       <  2:1229888/76‑1 (MQ=255)
                   ccGGATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGcc                                                      <  2:644455/76‑1 (MQ=255)
                    cGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCa                                                     >  1:4621313/1‑76 (MQ=255)
                      tATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCaaa                                                   <  2:1993490/76‑1 (MQ=255)
                      tATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCaaa                                                   <  2:886310/76‑1 (MQ=255)
                       aTCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAc                                                  <  2:1588167/76‑1 (MQ=255)
                        tCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAcc                                                 <  2:3611888/76‑1 (MQ=255)
                        tCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAcc                                                 >  1:656050/1‑76 (MQ=255)
                        tCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAcc                                                 <  1:3001792/76‑1 (MQ=255)
                         cGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCa                                                <  1:3613632/76‑1 (MQ=255)
                         cGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCa                                                <  1:216873/76‑1 (MQ=255)
                          catacctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAg                                               <  2:4590590/74‑1 (MQ=255)
                            tacctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGt                                             >  1:3564938/1‑76 (MQ=255)
                            tacctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGt                                             <  1:4700278/76‑1 (MQ=255)
                             acctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTg                                            <  1:3338982/76‑1 (MQ=255)
                              cctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGc                                           <  2:2469314/76‑1 (MQ=255)
                              cctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGc                                           <  2:2502560/76‑1 (MQ=255)
                              cctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGc                                           <  1:3262013/76‑1 (MQ=255)
                               ctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCa                                          >  1:3245457/1‑76 (MQ=255)
                                 acAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACt                                        <  1:4254841/76‑1 (MQ=255)
                                  cAAAGGAGCCAGCTGGAGGCCTTATAGCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTg                                       <  2:4301009/76‑1 (MQ=255)
                                  cAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACtc                                       >  2:563671/1‑75 (MQ=255)
                                  cAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTg                                       <  1:3039993/76‑1 (MQ=255)
                                  cAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTg                                       <  2:1227015/76‑1 (MQ=255)
                                  aaaaGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTg                                       <  1:139982/75‑1 (MQ=255)
                                   aaaGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGc                                      <  2:1465614/76‑1 (MQ=255)
                                     agaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTc                                    >  1:4191050/1‑76 (MQ=255)
                                       aagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGa                                  <  1:3181067/76‑1 (MQ=255)
                                       aagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGa                                  <  2:1694338/76‑1 (MQ=255)
                                       aagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGa                                  >  1:1278213/1‑76 (MQ=255)
                                            aGCTGGAAGCCTTATATCATCGCCGCGGTGAGACTGGCACGTTAATACGCCAAAACAGGCGCAGTGCTAGaacaac                             >  1:562126/1‑76 (MQ=255)
                                            aGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAAc                             <  2:2229088/76‑1 (MQ=255)
                                             gccggAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAAct                            <  2:333012/73‑1 (MQ=255)
                                              cTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActg                           <  2:1112411/76‑1 (MQ=255)
                                              cTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActg                           >  1:771542/1‑76 (MQ=255)
                                               tGGAAGTCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgc                          >  1:1633217/1‑76 (MQ=255)
                                                ggAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgct                         >  1:3891009/1‑76 (MQ=255)
                                                ggAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgct                         >  2:1810070/1‑76 (MQ=255)
                                                 gAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActgctg                        <  1:1542621/76‑1 (MQ=255)
                                                 gAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGCGCCCCGCTCGAGCAActgctg                        >  2:331876/1‑76 (MQ=255)
                                                   aGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaa                      <  2:1507373/76‑1 (MQ=255)
                                                     ccTTGTGTCAGCGCCGGGGTGAGCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCACCTGCTGAaag                    <  2:251920/76‑1 (MQ=255)
                                                      cTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaaga                   >  2:3219750/1‑76 (MQ=255)
                                                        tataTCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaagaag                 >  1:2495342/1‑76 (MQ=255)
                                                         atatCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAgc                <  1:2246/76‑1 (MQ=255)
                                                         atatCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAgc                >  2:4167268/1‑76 (MQ=255)
                                                          tatCAGCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAgcg               <  1:4821895/76‑1 (MQ=255)
                                                           atcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAgcgc              >  2:1622561/1‑76 (MQ=255)
                                                           atcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAgcgc              <  2:4068532/76‑1 (MQ=255)
                                                            tcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCa             <  2:2096536/76‑1 (MQ=255)
                                                             catcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCaa            >  1:3187712/1‑76 (MQ=255)
                                                             catcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCaa            >  1:1112508/1‑76 (MQ=255)
                                                             catcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCaa            >  2:1717935/1‑76 (MQ=255)
                                                              atcTCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAg           <  2:4834746/76‑1 (MQ=255)
                                                               tcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCACCTGCTGAAAGAAGCGCACGa          <  2:2002511/76‑1 (MQ=255)
                                                                cgccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGaa         <  2:3918232/76‑1 (MQ=255)
                                                                cgccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGaa         <  2:3245457/76‑1 (MQ=255)
                                                                 gccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAc        >  2:101213/1‑76 (MQ=255)
                                                                  ccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAcg       <  2:797724/76‑1 (MQ=255)
                                                                   cgcgATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGGAACTGCTGAAAGAAGCGCAAGAAcgc      <  2:1376900/76‑1 (MQ=255)
                                                                   cgcgATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAcgc      <  1:3788678/76‑1 (MQ=255)
                                                                   cgcgATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAcgc      >  2:3867884/1‑76 (MQ=255)
                                                                    gcgATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAcgcg     <  2:3012431/76‑1 (MQ=255)
                                                                    gcgATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAcgcg     >  2:2672012/1‑76 (MQ=255)
                                                                    gcgATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAcgcg     >  1:4535335/1‑76 (MQ=255)
                                                                    gcgATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAcgcg     <  1:1877120/76‑1 (MQ=255)
                                                                     tgatgatcCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGt    >  2:4098886/2‑76 (MQ=255)
                                                                     cgATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGt    <  1:3607991/76‑1 (MQ=255)
                                                                      gatgatcCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTg   >  1:1602039/1‑76 (MQ=255)
                                                                      gatgatcCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTg   >  2:3692350/1‑76 (MQ=255)
                                                                             |                                                                      
AAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGA  >  NC_000913/3005323‑3005468

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