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.C99.2% 356.0 / ‑3.6 133intergenic (‑51/‑248)ysgA/udphypothetical protein/uridine phosphorylase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (59/73);  total (60/73)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 147 total aligned reads displayed.

AACGATGGTCGAAGCGAGAGGCATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCT  >  NC_000913/4013696‑4013838
                                                                            |                                                                    
aaCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTa                                                                       <  1:4330995/76‑1 (MQ=255)
aaCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTa                                                                       >  2:3708723/1‑76 (MQ=255)
 aCGATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAc                                                                     <  1:2391114/76‑1 (MQ=255)
  tgATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAcc                                                                    >  2:4110009/2‑76 (MQ=255)
   gATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccc                                                                   >  2:1031699/1‑75 (MQ=255)
   gATGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccc                                                                   <  1:2034820/76‑2 (MQ=255)
     tGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccctt                                                                 >  2:4705074/1‑73 (MQ=255)
     tGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccctt                                                                 >  1:2878326/1‑73 (MQ=255)
     tGGTCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTAccctt                                                                 <  1:2194445/76‑4 (MQ=255)
        tCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTtata                                                               >  1:3135495/1‑76 (MQ=255)
        tCGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTtata                                                               <  1:22832/76‑1 (MQ=255)
         cGAAGGGAGAGGCGAGGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaa                                                              <  2:261203/76‑1 (MQ=255)
         cGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaa                                                              <  1:3820872/76‑1 (MQ=255)
         cGAAGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaa                                                              >  1:2893933/1‑76 (MQ=255)
           aaGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaa                                                            >  1:3947962/1‑76 (MQ=255)
            aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaaa                                                           <  2:3178773/76‑1 (MQ=255)
            aGCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATaaaaa                                                           <  2:4046991/76‑1 (MQ=255)
             gCGAGAGGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAAt                                                          <  1:1348947/76‑1 (MQ=255)
                agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTa                                                                       <  1:3961743/60‑1 (MQ=255)
                agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTa                                                                       >  2:3961743/1‑60 (MQ=255)
                agagGCGATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGtt                                                       <  2:2092265/76‑1 (MQ=255)
                    gcgaTGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTgcgc                                                   <  2:3268658/73‑1 (MQ=255)
                     ggaTGCGGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCa                                                  <  2:4862972/74‑1 (MQ=255)
                     cgaTGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCa                                                  >  1:425616/3‑76 (MQ=255)
                      gaggcagcaGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCaa                                                 <  2:962310/73‑1 (MQ=255)
                        tGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATg                                               <  1:384148/76‑1 (MQ=255)
                        tGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATg                                               >  2:2501993/1‑76 (MQ=255)
                         gcagcaGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGt                                              <  2:425616/76‑1 (MQ=255)
                         gcagcaGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGt                                              >  1:4726262/1‑76 (MQ=255)
                         gcagcaGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGt                                              >  2:1002673/1‑76 (MQ=255)
                         gcagcaGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGt                                              <  1:1712867/76‑1 (MQ=255)
                          gagcagGTGCAAATCCGGATTGTGGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGtt                                             <  2:4843600/75‑1 (MQ=255)
                          cagcagGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGtt                                             >  1:1835113/1‑76 (MQ=255)
                           agcagGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTa                                            >  1:511463/1‑76 (MQ=255)
                           agcagGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTa                                            <  2:1875899/76‑1 (MQ=255)
                           agcagGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTa                                            <  2:628294/76‑1 (MQ=255)
                           agcagGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTa                                            <  2:4364430/76‑1 (MQ=255)
                            gcagGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTaa                                           <  2:3876858/76‑1 (MQ=255)
                            gcagGTGCAAATCCAGATTGTTGGGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTaa                                           <  2:246043/76‑1 (MQ=255)
                             cagGGGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAAc                                          <  1:2564259/76‑1 (MQ=255)
                               ggTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACta                                        <  2:1192090/76‑1 (MQ=255)
                               ggTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACta                                        >  2:1603746/1‑76 (MQ=255)
                               ggTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACta                                        >  1:4085659/1‑76 (MQ=255)
                                gTGCAAATCCAGATTGTTGTTTTTTTGCCATGGTATTCTCCTTACCCTTATAAAAATGTTGCTCAATGTTAACtat                                       <  1:2783451/76‑1 (MQ=255)
                                gTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACtat                                       <  1:4543906/76‑1 (MQ=255)
                                 tgaaaaTCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACtata                                      >  2:3563307/4‑76 (MQ=255)
                                 tGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACtata                                      <  2:2992905/76‑1 (MQ=255)
                                 tGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACtata                                      <  2:3790750/76‑1 (MQ=255)
                                   cAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGt                                    <  1:446457/76‑1 (MQ=255)
                                   cAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGt                                    >  2:4045379/1‑76 (MQ=255)
                                    aaaTCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTc                                   <  2:2220161/76‑1 (MQ=255)
                                    aaaTCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTc                                   >  1:3223067/1‑76 (MQ=255)
                                     aaTCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCa                                  <  1:872879/76‑1 (MQ=255)
                                      aTCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAg                                 <  2:1172903/76‑1 (MQ=255)
                                       tCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGc                                >  1:3083792/1‑76 (MQ=255)
                                        ccATATTGTTGTGTTGTTGCGATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATTTTAAGTATAGTCAGCa                               <  2:2545619/76‑1 (MQ=255)
                                        ccAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCa                               >  1:4070493/1‑76 (MQ=255)
                                         cAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCAt                              >  1:2940994/1‑76 (MQ=255)
                                         cAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCAt                              >  1:918769/1‑76 (MQ=255)
                                         cAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCAt                              <  2:2430886/76‑1 (MQ=255)
                                         cAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCAt                              <  1:1766740/76‑1 (MQ=255)
                                         cAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCAt                              <  1:1200546/76‑1 (MQ=255)
                                         cAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCAt                              <  1:1134717/76‑1 (MQ=255)
                                          aGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATg                             >  2:3323022/1‑76 (MQ=255)
                                          aGATGGGTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATg                             <  2:4476673/76‑1 (MQ=255)
                                           gATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGc                            <  2:1593888/76‑1 (MQ=255)
                                           gATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGc                            >  1:3299029/1‑76 (MQ=255)
                                           gATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGc                            >  1:4631232/1‑76 (MQ=255)
                                            aTTGTTGTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGca                           <  2:4035244/76‑1 (MQ=255)
                                             ttgttgTGTTGTTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGcaa                          >  2:2720839/1‑76 (MQ=255)
                                               gttgtgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGcaaca                        <  2:2940994/76‑1 (MQ=255)
                                                ttgtgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGcaacaa                       <  1:4110009/76‑1 (MQ=255)
                                                 tgtgttgtTGCCATGGTAGTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACaaa                      <  2:511921/76‑1 (MQ=255)
                                                  gtgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAAt                     <  2:2261739/76‑1 (MQ=255)
                                                  gtgttgtTGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAAt                     <  2:3549121/76‑1 (MQ=255)
                                                    gttgttGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATca                   >  1:1124302/1‑76 (MQ=255)
                                                    gttgttGCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATca                   >  2:4518256/1‑76 (MQ=255)
                                                     ttgttgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATcac                  <  1:3888032/76‑1 (MQ=255)
                                                     ttgttgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATcac                  >  1:3665731/1‑76 (MQ=255)
                                                       gttgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACAt                >  2:1884045/1‑76 (MQ=255)
                                                       gttgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACAt                >  2:4457009/1‑76 (MQ=255)
                                                        ttgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACAtt               <  2:2280842/76‑1 (MQ=255)
                                                        ttgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACAtt               >  2:4101429/1‑76 (MQ=255)
                                                         tgCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTg              <  2:1516338/76‑1 (MQ=255)
                                                          gCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGc             <  1:1327386/76‑1 (MQ=255)
                                                          gCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGc             >  1:4037404/1‑76 (MQ=255)
                                                          gCCATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGc             >  1:3942443/1‑76 (MQ=255)
                                                           ccATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGcc            >  2:3461061/1‑76 (MQ=255)
                                                           ccATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGcc            >  2:4834473/1‑76 (MQ=255)
                                                            cATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCACCAAATCACATTGCCt           <  1:1002674/76‑1 (MQ=255)
                                                            cATGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCt           <  2:3312867/76‑1 (MQ=255)
                                                             aTGGTATTCTCCGTACCCTTATAAAAATGTTTCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTg          >  2:1962042/1‑76 (MQ=255)
                                                             aTGGTATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTg          >  2:3190249/1‑76 (MQ=255)
                                                                 tATTCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATc      <  1:3496029/76‑1 (MQ=255)
                                                                  atgctcCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCg     <  2:1966841/73‑1 (MQ=255)
                                                                   ttCTCCGTACCCTTATAAAAATGTTGCGCAATGTTATCTATAGTCAGCATGCAACAAATCACATTGCCTGAATCgg    <  2:2901994/76‑1 (MQ=255)
                                                                   ttCTCCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCgg    <  2:2893933/76‑1 (MQ=255)
                                                                    tctcCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGc   >  1:509484/1‑76 (MQ=255)
                                                                    tctcCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGc   >  1:2007861/1‑76 (MQ=255)
                                                                     ctcCGTACCCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCGCATTGCCTGAATCGGCt  >  2:3474318/1‑76 (MQ=255)
                                                                            |                                                                    
AACGATGGTCGAAGCGAGAGGCATGCAGCAGGTGCAAATCCAGATTGTTGTGTTGTTGCCATGGTATTCTCCGTACCTTATAAAAATGTTGCGCAATGTTAACTATAGTCAGCATGCAACAAATCACATTGCCTGAATCGGCT  >  NC_000913/4013696‑4013838

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