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.C98.7% 409.3 / ‑1.3 154coding (1043/1092 nt)ygeWputative carbamoyl transferase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (68/83);  total (71/83)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 165 total aligned reads displayed.

AAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGA  >  NC_000913/3005323‑3005468
                                                                             |                                                                      
aaGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGACAGCTGGAAGCCTTATATAAGCGCCGCgatgat                                                                         >  1:5064638/1‑76 (MQ=255)
       aTTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTcc                                                                  <  2:4694821/76‑1 (MQ=255)
        ttCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTccc                                                                 <  2:1244360/76‑1 (MQ=255)
        ttCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTccc                                                                 <  2:1149090/76‑1 (MQ=255)
         tCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCg                                                                >  1:5835186/1‑76 (MQ=255)
         tCGAAAAATACCGTATCGCGACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCg                                                                <  1:1651883/76‑1 (MQ=255)
          cGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGt                                                               >  2:473992/1‑76 (MQ=255)
           gAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTa                                                              >  2:2567204/1‑76 (MQ=255)
            aaaaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTaa                                                             >  1:492679/1‑76 (MQ=255)
            aaaaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTaa                                                             <  2:4047153/76‑1 (MQ=255)
              aaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTaaa                                                            >  1:545925/1‑75 (MQ=255)
              aaaTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAAt                                                           <  2:3612148/76‑1 (MQ=255)
                aTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATAc                                                         >  2:2858725/1‑76 (MQ=255)
                aTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATAc                                                         >  1:4308791/1‑76 (MQ=255)
                aTACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATAc                                                         <  1:895438/76‑1 (MQ=255)
                aTACCGTATCGCTACCTACAAAGAAGCCAGCTGCAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATAc                                                         >  1:3621996/1‑76 (MQ=255)
                 tACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTACATACg                                                        >  2:3147689/1‑76 (MQ=255)
                 tACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACg                                                        <  1:47307/76‑1 (MQ=255)
                  aCCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGc                                                       <  2:2317831/76‑1 (MQ=255)
                    cGTATCGCGCCCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCa                                                     <  2:1543222/76‑1 (MQ=255)
                     gTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCaa                                                    <  1:1998329/76‑1 (MQ=255)
                     gTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCaa                                                    <  1:177485/76‑1 (MQ=255)
                     gTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCaa                                                    >  1:4991252/1‑76 (MQ=255)
                       aTCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAc                                                  >  2:86379/1‑76 (MQ=255)
                       aTCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAc                                                  <  1:3488361/76‑1 (MQ=255)
                       aTCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAc                                                  >  1:5568367/1‑76 (MQ=255)
                       aTCGCTACCTACAAAGAAGCCAGCTGGAAGACTTATATCAGCGCCGCGATGATCCTGTTCCGCACATACGACAAAc                                                  >  2:304438/1‑76 (MQ=255)
                       aTCGCTACCTACAAAGAAGCCAACTGGAAGCCTTATATCATCGCCGCGATGATCCTCTCCCGCAAATACGCCAAAc                                                  >  2:691875/1‑76 (MQ=255)
                        tCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAcc                                                 <  1:5231456/76‑1 (MQ=255)
                        tCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAcc                                                 <  2:4429583/76‑1 (MQ=255)
                        tCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAAcc                                                 >  2:3996997/1‑76 (MQ=255)
                        tCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCAGCGATGATCCTGTCCCGTAAATACTCTAAAcc                                                 >  1:162386/1‑76 (MQ=255)
                        tCGCTAAATACAAACAAGCCAGCGGGAAGCCGTATATCATCGCCGCGAGGAGCCTGGCCAGCAAACACGCCAAAcc                                                 >  2:5426690/1‑76 (MQ=255)
                         cGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCACACCa                                                >  1:4402712/1‑76 (MQ=255)
                         cGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCa                                                <  2:4696254/76‑1 (MQ=255)
                         cGCTACCTAATACTAAGCCAGCTCGAAGCCTTATATCATCGGCGCGATGATCCTGTCCCGTAAATACGCCAAACCa                                                <  1:3411833/76‑1 (MQ=255)
                          gCTACCTACAAAGCAGCCAGCTGGACGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAACTACGCCACAcccg                                               >  2:1499847/1‑74 (MQ=255)
                          gCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAACTACGCCAAACCAg                                               >  2:2552178/1‑76 (MQ=255)
                          gCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAg                                               >  2:2512924/1‑76 (MQ=255)
                          gCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAg                                               >  1:1906546/1‑76 (MQ=255)
                            tacctacAAAGAAGCTAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAggg                                             >  2:2395193/1‑75 (MQ=255)
                            tacctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGt                                             <  2:4648614/76‑1 (MQ=255)
                            taccgacAAAGAGGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGt                                             <  1:1801857/76‑1 (MQ=255)
                              cctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGc                                           >  1:3412525/1‑76 (MQ=255)
                               ctacAACGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAACTACGCCAAACCCGGTGCa                                          >  1:4004575/1‑76 (MQ=255)
                               ctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCa                                          >  1:3533314/1‑76 (MQ=255)
                                tacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCAc                                         <  1:5288813/76‑1 (MQ=255)
                                 acAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACt                                        >  2:2657723/1‑76 (MQ=255)
                                 acAAAGAAGCCAGCTGGAAGCCTTATAGCATCGCCGCGATGAACCTGTCCCGTAAATACTCCAAACCAGCTGCAcc                                        >  2:5274577/1‑75 (MQ=255)
                                  cAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTg                                       <  1:2778563/76‑1 (MQ=255)
                                  cAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTg                                       <  1:3304412/76‑1 (MQ=255)
                                  cAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTg                                       <  2:2251395/76‑1 (MQ=255)
                                   aaaGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGc                                      <  1:4126736/76‑1 (MQ=255)
                                    aggaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGTTGCACTGCt                                     <  1:2932653/74‑1 (MQ=255)
                                    aagaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCt                                     <  2:5483849/76‑1 (MQ=255)
                                     agaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTc                                    <  2:1757213/76‑1 (MQ=255)
                                     agaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTc                                    <  1:1436095/76‑1 (MQ=255)
                                     agaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTc                                    <  2:3280145/76‑1 (MQ=255)
                                     agaagCCAGCGGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTc                                    <  2:2182864/76‑1 (MQ=255)
                                      gaagCCAGCTGGAGTCCTTATATCATCGCTGCGATGATCCTGTCCCGTAAATACGCCACACCAGGTGCGCTGCTca                                   <  2:4446279/76‑2 (MQ=255)
                                      gaagCCAGCTGGAAGCCTTCTATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCg                                   >  1:4815050/1‑76 (MQ=255)
                                      gaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCg                                   >  2:2796961/1‑76 (MQ=255)
                                      gaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCg                                   >  1:4766219/1‑76 (MQ=255)
                                       aagCCTGCTGGAGTCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGa                                  <  1:1098711/76‑1 (MQ=255)
                                       aagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGa                                  <  1:3058462/76‑1 (MQ=255)
                                         gCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGc                                >  1:1726649/1‑76 (MQ=255)
                                         gCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGc                                <  1:3157751/76‑1 (MQ=255)
                                          ccAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCa                               >  1:460910/1‑76 (MQ=255)
                                            aGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAAc                             >  1:2568240/1‑76 (MQ=255)
                                             gCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAAct                            >  1:598613/1‑76 (MQ=255)
                                              cTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActg                           <  2:712150/76‑1 (MQ=255)
                                              cTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAActg                           <  1:4347644/76‑1 (MQ=255)
                                               tGGAAGCCTTATATCATCGCCGCGATGATACTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAAccgc                          >  2:5787195/1‑76 (MQ=255)
                                                  aaGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGa                       >  2:1310155/1‑76 (MQ=255)
                                                   aGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaa                      >  1:2870516/1‑76 (MQ=255)
                                                    gCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaaa                     <  1:1052889/76‑1 (MQ=255)
                                                    gCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGaaa                     <  1:86379/76‑1 (MQ=255)
                                                      cgtataTCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaaga                   <  2:3663655/74‑1 (MQ=255)
                                                      cTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAACCCAGGTGCTCTGCTCGAGCAAGTGCTGAaaga                   <  2:5064637/76‑1 (MQ=255)
                                                      cTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaaga                   >  2:240942/1‑76 (MQ=255)
                                                      cTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaaga                   >  1:2958660/1‑76 (MQ=255)
                                                       ttATATCATCGCCTCGACGACGCGGTACCGCAAATACGCCAAACCACGAAAACTGCGCCAGCAACTGCTGgaagaa                  >  2:3225493/1‑76 (MQ=255)
                                                       ttATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAaagaa                  >  2:1871281/1‑76 (MQ=255)
                                                         atatCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAgc                <  1:2200854/76‑1 (MQ=255)
                                                         atatCACCGCCGCGATGATCCTGTCCCGTAAATACGCCAAGCCAGGTGCACTGCTCGAGCAACTTCTGAAAGAAgc                <  1:657347/76‑1 (MQ=255)
                                                          tatCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAgcg               >  2:2002740/1‑76 (MQ=255)
                                                          tatCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAgcg               >  1:5426412/1‑76 (MQ=255)
                                                            tcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCa             >  2:1909752/1‑76 (MQ=255)
                                                            tcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCa             >  2:1966691/1‑76 (MQ=255)
                                                             cagcgccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGATCAGCTGCTGAAAGAAGCGCaa            <  2:1610777/73‑1 (MQ=255)
                                                               tcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGa          >  1:3927803/1‑76 (MQ=255)
                                                                 gccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAc        <  1:5686946/76‑1 (MQ=255)
                                                                 gccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAc        >  1:473892/1‑76 (MQ=255)
                                                                  ccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAcg       <  2:436366/76‑1 (MQ=255)
                                                                   cgcgATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAcgc      <  1:1310155/76‑1 (MQ=255)
                                                                   cgcgATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAAcgc      <  1:4946381/76‑1 (MQ=255)
                                                                     cgATGATCCTGTCCCGTAAATTCGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGt    <  2:5568365/76‑1 (MQ=255)
                                                                      gatgatcCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTg   <  1:2063604/76‑1 (MQ=255)
                                                                       atgatcCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGa  <  2:1580229/76‑1 (MQ=255)
                                                                       atgatcCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGa  >  1:1675688/1‑76 (MQ=255)
                                                                             |                                                                      
AAGGCGTATTCGAAAAATACCGTATCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGA  >  NC_000913/3005323‑3005468

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