Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,013,772 +C intergenic (‑53/‑246) ysgA ← / → udp hypothetical protein/uridine phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,013,7701.C98.5% 354.8 / ‑1.6 131intergenic (‑51/‑248)ysgA/udphypothetical protein/uridine phosphorylase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (58/71);  total (59/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

AACGATGGTCGAAGCGAGAGGCATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTC  >  NC_000913/4013696‑4013839
                                                                             |                                                                      
aaCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTa                                                                         <  1:979279/76‑1 (MQ=255)
 ccGAGGGTGGAAGCGAGAGGCGGTGCAGCAGGTGCCCAGCCAGATTGTTGTTTTGTTGCCATGGTATTCTCCGTAc                                                                       <  1:831402/75‑1 (MQ=255)
 aCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCCTGGTATTCTCCGTAc                                                                       >  2:121054/1‑76 (MQ=255)
  cGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAcc                                                                      <  1:2260185/76‑1 (MQ=255)
   gATGGTCGAAGCGGGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccc                                                                     <  1:5698577/76‑2 (MQ=255)
   gATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccc                                                                     <  1:250939/76‑2 (MQ=255)
   gATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccc                                                                     <  2:4209006/76‑2 (MQ=255)
    aTGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccct                                                                    <  1:260926/76‑3 (MQ=255)
     tGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATGGTTGTGTTGTTGCCATGGTATTGTCCGTAccctt                                                                   <  2:791744/76‑4 (MQ=255)
         cGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaa                                                               >  1:1297728/1‑76 (MQ=255)
          gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaa                                                              <  1:1934974/76‑1 (MQ=255)
          gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaa                                                              <  2:1481556/76‑1 (MQ=255)
          gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaa                                                              >  1:5282044/1‑76 (MQ=255)
          gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaa                                                              >  2:3203476/1‑76 (MQ=255)
          gAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaa                                                              >  1:2791574/1‑76 (MQ=255)
           aaGCGAGAGGCGATGCAGCAGGTGCAAATCCCGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaa                                                             <  1:1224137/76‑1 (MQ=255)
           aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaa                                                             <  2:299599/76‑1 (MQ=255)
           aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaa                                                             >  2:747607/1‑76 (MQ=255)
            aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaaa                                                            <  2:4719208/76‑1 (MQ=255)
            aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaaa                                                            >  2:1008813/1‑76 (MQ=255)
            aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaaa                                                            <  1:802022/76‑1 (MQ=255)
             gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt                                                           <  2:3491319/76‑1 (MQ=255)
             gCGAGAGGCGATGCAGCAGGGGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt                                                           >  1:1249086/1‑76 (MQ=255)
              cGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATg                                                          >  2:1597425/1‑76 (MQ=255)
              cGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATg                                                          >  1:346773/1‑76 (MQ=255)
               gagaGGCGATGCAGCATGTTCAAATCCAGATTGTTGTGTTGTTGCTATGGTATTCTCCCTACCCTTATATAAATGt                                                         <  1:1035332/76‑1 (MQ=255)
                    gcgaTGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTgcgc                                                    >  2:2530566/4‑76 (MQ=255)
                    gcgaTGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTgcgc                                                    <  1:5338146/73‑1 (MQ=255)
                     cgaTGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCa                                                   >  1:105825/3‑76 (MQ=255)
                     cgaTGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCa                                                   >  2:693404/3‑76 (MQ=255)
                        tGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATg                                                >  1:1631986/1‑76 (MQ=255)
                        tGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATg                                                >  1:1936276/1‑76 (MQ=255)
                        tGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAcc                                                                      >  1:4792134/1‑54 (MQ=255)
                        tGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAcc                                                                      <  2:4792133/54‑1 (MQ=255)
                         gcagcaGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGt                                               <  1:4195048/76‑1 (MQ=255)
                         gcagcaGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGt                                               <  2:5012260/76‑1 (MQ=255)
                         gcagcaGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGt                                               <  1:3682187/76‑1 (MQ=255)
                          cagcagGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGtt                                              >  1:2714832/1‑76 (MQ=255)
                             cagGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAAc                                           <  1:4883695/76‑1 (MQ=255)
                              agGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACt                                          >  2:3607656/1‑76 (MQ=255)
                              agGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACt                                          <  2:1423052/76‑1 (MQ=255)
                               ggTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACta                                         <  2:1263414/76‑1 (MQ=255)
                                 tGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACgata                                       >  2:5809123/1‑76 (MQ=255)
                                  gCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAg                                      <  1:938048/76‑1 (MQ=255)
                                  gCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAg                                      <  1:3879329/76‑1 (MQ=255)
                                   cAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGt                                     >  2:2814193/1‑76 (MQ=255)
                                   cAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGt                                     >  2:5027580/1‑76 (MQ=255)
                                    aaaTCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTc                                    <  1:3857520/76‑1 (MQ=255)
                                     aagccAGATTGTTGTGTTGTTGCCATGGTATCCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAGCTATAGTCa                                   <  1:2946163/73‑1 (MQ=255)
                                     aaTCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCa                                   <  2:4033520/76‑1 (MQ=255)
                                      aTCCAGATTGTTGTGTTGTTGCGTTGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAg                                  <  2:4521426/76‑1 (MQ=255)
                                      aTCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAg                                  <  2:4981608/76‑1 (MQ=255)
                                      aTCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAg                                  >  2:2012141/1‑76 (MQ=255)
                                       tCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGc                                 <  2:4228635/76‑1 (MQ=255)
                                       tCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGc                                 <  1:3497624/76‑1 (MQ=255)
                                        ccAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCa                                <  1:2530566/76‑1 (MQ=255)
                                         tggATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCAt                               <  2:5430725/74‑1 (MQ=255)
                                         cAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCAt                               >  2:5106437/1‑76 (MQ=255)
                                         cAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCAt                               >  1:403169/1‑76 (MQ=255)
                                          aGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATg                              >  2:3394206/1‑76 (MQ=255)
                                           gATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGc                             <  1:4546677/76‑1 (MQ=255)
                                              tgttgtgttgTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGcaac                          >  2:2940365/1‑76 (MQ=255)
                                              tgttgtgttgTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGcaac                          <  2:1717886/76‑1 (MQ=255)
                                              tgttgtgttgTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAAGGTTAACTATAGTCAGCATGcaac                          <  2:2714832/76‑1 (MQ=255)
                                               gttgtgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGcaaca                         >  1:1666588/1‑76 (MQ=255)
                                                 tgtgttgtTGCCATTGTATTCTCCGTACCCTTATAAAAATGTTGCGTAATGTTAACTATAGTCAGCATGCGACaaa                       <  2:2694133/76‑1 (MQ=255)
                                                 tgtgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGcaacaca                       >  2:23133/1‑74 (MQ=255)
                                                  gtgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAAt                      <  2:829964/76‑1 (MQ=255)
                                                   tgttgtTGCTATGGCATTGTCCGTGCCATTTTGAAAATGTTGCGCAATGTTAACTATAGTCAGCACGCAACAAATc                     >  1:2057404/1‑76 (MQ=255)
                                                   tgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATc                     <  2:3363499/76‑1 (MQ=255)
                                                   tgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATc                     >  2:5831998/1‑76 (MQ=255)
                                                     ttgttgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATcac                   <  2:3500744/76‑1 (MQ=255)
                                                      tgttgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATcaca                  <  1:3170283/76‑1 (MQ=255)
                                                       gttgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACAt                 <  1:790814/76‑1 (MQ=255)
                                                       gttgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACAt                 >  2:3023559/1‑76 (MQ=255)
                                                         tgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTg               >  1:4873066/1‑76 (MQ=255)
                                                         tgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTg               <  1:3066372/76‑1 (MQ=255)
                                                          gCCATGGTATTCTCCGTACCCGTATAAAAATGTTGCGCAATGTGAACTATAGGCAGCATGCAACAAATCACATTGc              >  1:4070570/1‑76 (MQ=255)
                                                            cATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCt            <  2:5017939/76‑1 (MQ=255)
                                                             aTGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATc                     >  2:4350766/1‑66 (MQ=255)
                                                             aTGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTg           <  1:4522241/76‑1 (MQ=255)
                                                             aTGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTg           >  1:5582623/1‑76 (MQ=255)
                                                              tGGTATTCTCCTTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGa          <  2:2945640/76‑1 (MQ=255)
                                                              tGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGa          <  1:4164481/76‑1 (MQ=255)
                                                               ggTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGaa         >  1:5582431/1‑76 (MQ=255)
                                                                gTATTCTCCGTCCCCTTATAAAAAGGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAAt        <  1:2608719/76‑1 (MQ=255)
                                                                 tATTCTCCGTGCCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATc       <  2:3236361/76‑1 (MQ=255)
                                                                 tATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATc       >  2:453149/1‑76 (MQ=255)
                                                                 tATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATc       <  2:1631986/76‑1 (MQ=255)
                                                                 tATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATc       <  2:1762888/76‑1 (MQ=255)
                                                                 tATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATc       <  1:1370510/76‑1 (MQ=255)
                                                                 tATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATc       <  1:1070157/76‑1 (MQ=255)
                                                                  aTTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCg      <  2:1116550/76‑1 (MQ=255)
                                                                   ttCTCCGTACCCTTAAAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCgg     <  2:2716657/76‑1 (MQ=255)
                                                                    tctcCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGc    <  1:3680648/76‑1 (MQ=255)
                                                                    tctcCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGc    <  1:610343/76‑1 (MQ=255)
                                                                     ctcCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCt   >  2:543101/1‑76 (MQ=255)
                                                                      tcCGTACCCTTCTAAAATGGTTGCGCAATGTTACTTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTc  <  2:5212650/76‑1 (MQ=255)
                                                                      tcCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTc  >  2:2043679/1‑76 (MQ=255)
                                                                      tcCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTc  >  2:240431/1‑76 (MQ=255)
                                                                             |                                                                      
AACGATGGTCGAAGCGAGAGGCATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCTC  >  NC_000913/4013696‑4013839

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