Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,090,282 C→T intergenic (+30/‑65) yggI → / → endA conserved protein YggI/DNA‑specific endonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,090,2820CT99.4% 667.5 / NA 171intergenic (+30/‑65)yggI/endAconserved protein YggI/DNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (86/84);  minor base G (0/1);  total (86/85)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 171 total aligned reads displayed.

CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGT  >  NC_000913/3090235‑3090328
                                               |                                              
cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              <  2:940028/50‑1 (MQ=255)
cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              <  2:4991509/50‑1 (MQ=255)
cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              <  2:4718633/50‑1 (MQ=255)
 aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                             >  1:679583/1‑50 (MQ=255)
 aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                             >  2:3379440/1‑50 (MQ=255)
   cTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           >  1:3132974/1‑50 (MQ=255)
   cTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           >  2:3758297/1‑50 (MQ=255)
    tGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           >  2:3476188/1‑49 (MQ=255)
     ggTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGa                                         >  2:664319/1‑50 (MQ=255)
     ggTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGa                                         >  1:138298/1‑50 (MQ=255)
      gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTg                                          <  1:4826927/48‑1 (MQ=255)
      gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGa                                         <  2:1225310/49‑1 (MQ=255)
       ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       <  1:2335949/50‑1 (MQ=255)
       ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       >  2:2646551/1‑50 (MQ=255)
       ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       >  2:3838433/1‑50 (MQ=255)
       ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       <  2:3855896/50‑1 (MQ=255)
       ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       >  1:1496124/1‑50 (MQ=255)
          cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                     <  1:1050440/49‑1 (MQ=255)
          cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                     <  1:5454223/49‑1 (MQ=255)
          cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    <  2:2797265/50‑1 (MQ=255)
          cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    <  2:2715176/50‑1 (MQ=255)
          cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    >  1:2117895/1‑50 (MQ=255)
          cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    >  1:3816122/1‑50 (MQ=255)
           gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    >  1:1396758/1‑49 (MQ=255)
           gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCAt                                   <  2:1255489/50‑1 (MQ=255)
           gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCAt                                   <  2:1792795/50‑1 (MQ=255)
           gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCAt                                   <  2:1703998/50‑1 (MQ=255)
            aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  >  1:5589906/1‑50 (MQ=255)
            aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  <  1:4100586/50‑1 (MQ=255)
            aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  >  2:3689383/1‑50 (MQ=255)
               tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                               <  1:4513555/50‑1 (MQ=255)
                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                              >  1:4163539/1‑49 (MQ=255)
                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                              >  1:2698733/1‑49 (MQ=255)
                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                              >  2:3257388/1‑49 (MQ=255)
                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                              <  1:4708459/49‑1 (MQ=255)
                  ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            <  1:3708146/50‑1 (MQ=255)
                  ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            <  1:4438747/50‑1 (MQ=255)
                   cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            >  2:362358/1‑49 (MQ=255)
                   cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            >  1:2184732/1‑49 (MQ=255)
                   cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                           <  2:4602580/50‑1 (MQ=255)
                   cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                           <  1:1555480/50‑1 (MQ=255)
                   cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                           <  1:2846158/50‑1 (MQ=255)
                   cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                           <  1:3323403/50‑1 (MQ=255)
                    aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            >  2:3962028/1‑48 (MQ=255)
                      cTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                        >  1:2019688/1‑50 (MQ=255)
                      cTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                        >  1:1077390/1‑50 (MQ=255)
                       tGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                       >  1:3767598/1‑50 (MQ=255)
                       tGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                       >  1:2511186/1‑50 (MQ=255)
                        gAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                      >  1:1928723/1‑50 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                         <  2:2550215/46‑1 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                        <  2:2379944/47‑1 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     <  1:1763423/50‑1 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                      <  1:3004445/49‑1 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                       <  2:2192783/48‑1 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                       <  2:5352866/48‑1 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                       <  2:5467394/48‑1 (MQ=255)
                          aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    <  2:1993054/50‑1 (MQ=255)
                          aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    <  2:2445278/50‑1 (MQ=255)
                            tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  2:5560809/1‑50 (MQ=255)
                            tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  1:4605237/1‑50 (MQ=255)
                             aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  2:1917960/1‑49 (MQ=255)
                             aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  1:2646714/1‑49 (MQ=255)
                             aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  2:2696036/1‑49 (MQ=255)
                             aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  2:2818010/1‑49 (MQ=255)
                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCt                 <  2:1970549/49‑1 (MQ=255)
                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                <  1:1227228/50‑1 (MQ=255)
                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                <  2:57888/50‑1 (MQ=255)
                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                <  2:5676369/50‑1 (MQ=255)
                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                >  1:245931/1‑50 (MQ=255)
                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                >  2:2230740/1‑50 (MQ=255)
                               cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               >  2:4980249/1‑50 (MQ=255)
                                aGGTACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg              >  2:748628/1‑50 (MQ=255)
                                aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg              >  1:5512393/1‑50 (MQ=255)
                                 ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt             >  2:2822640/1‑50 (MQ=255)
                                 ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt             >  2:2494863/1‑50 (MQ=255)
                                  gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            >  2:2370843/1‑50 (MQ=255)
                                  gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            <  1:4594334/50‑1 (MQ=255)
                                   aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg           <  1:4931761/50‑1 (MQ=255)
                                    aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc          <  1:5129360/50‑1 (MQ=255)
                                    aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc          <  1:2075236/50‑1 (MQ=255)
                                    aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc          >  1:4280452/1‑50 (MQ=255)
                                      tttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGGTACGTTGCTg        >  2:2953943/1‑50 (MQ=255)
                                      tttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg        >  2:911064/1‑50 (MQ=255)
                                       ttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg       >  2:162214/1‑50 (MQ=255)
                                        tCCTGATTTGGCTGATTGCATCCCAAAACAGCTTTCGCTACGTTGCTg        <  1:3645094/48‑1 (MQ=38)
                                        tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg       <  2:5086886/49‑1 (MQ=255)
                                        tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg        <  1:633224/48‑1 (MQ=255)
                                        tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg        <  1:837789/48‑1 (MQ=255)
                                        tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc      >  1:965977/1‑50 (MQ=255)
                                        tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc      >  1:5310617/1‑50 (MQ=255)
                                        tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc      >  2:1229112/1‑50 (MQ=255)
                                        tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc      >  1:3608182/1‑50 (MQ=255)
                                        tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCCACGTTGCTg        <  1:2391119/48‑1 (MQ=38)
                                         ccTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt     >  2:4166684/1‑50 (MQ=255)
                                         ccTGATGTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc      <  2:3386443/49‑1 (MQ=255)
                                          cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt     <  2:4343736/49‑1 (MQ=255)
                                          cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc    >  2:760818/1‑50 (MQ=255)
                                           tGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCg   >  2:4496237/1‑50 (MQ=255)
                                           tGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCg   >  1:3663321/1‑50 (MQ=255)
                                            gATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt  >  2:717183/1‑50 (MQ=255)
                                               |                                              
CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGT  >  NC_000913/3090235‑3090328

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: