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.1% 433.8 / NA 112intergenic (+30/‑65)yggI/endAconserved protein YggI/DNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (63/48);  total (64/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 112 total aligned reads displayed.

CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGT  >  NC_000913/3090235‑3090328
                                               |                                              
cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              <  1:819682/50‑1 (MQ=255)
cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCGGATTTg                                              >  2:3159435/1‑50 (MQ=38)
 aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                             <  1:1780345/50‑1 (MQ=255)
 aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                             <  1:1355040/50‑1 (MQ=255)
  gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                            >  2:612042/1‑50 (MQ=255)
  gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                            >  1:1926768/1‑50 (MQ=255)
  gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                            <  1:729086/50‑1 (MQ=255)
   cTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           >  2:244385/1‑50 (MQ=255)
   cTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           <  1:2669951/50‑1 (MQ=255)
      gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGa                                         <  1:3752973/49‑1 (MQ=255)
      gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                        <  1:135438/50‑1 (MQ=255)
      gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                        <  1:2464124/50‑1 (MQ=255)
      gTTGCGAAATAACCATCTGAACCATCAGGAACTTTCCTGATTTGGCTGAt                                        >  2:1677998/1‑50 (MQ=255)
       ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       <  2:1005792/50‑1 (MQ=255)
       ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       <  2:4192130/50‑1 (MQ=255)
         gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                      >  1:1496454/1‑49 (MQ=255)
         gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                      <  1:1906952/49‑1 (MQ=255)
          cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                     <  1:2235568/49‑1 (MQ=255)
          cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    >  1:4463564/1‑50 (MQ=255)
           gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    >  2:3750520/1‑49 (MQ=255)
           gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    >  1:1836753/1‑49 (MQ=255)
           gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCAt                                   <  1:2344420/50‑1 (MQ=255)
           gAAATAACCATCTGAACTATCAGGAACTTTCCAGATCTGGCTGATTGCa                                    >  2:4410549/1‑49 (MQ=255)
            aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  <  1:2544827/50‑1 (MQ=255)
            aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  >  2:2875639/1‑50 (MQ=255)
             aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 >  1:4526942/1‑50 (MQ=255)
             aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 >  2:1986763/1‑50 (MQ=255)
             aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 >  2:2019465/1‑50 (MQ=255)
              aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAcc                                <  1:988950/50‑1 (MQ=255)
              aTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 >  2:3182446/1‑49 (MQ=255)
               tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                               <  2:3334764/50‑1 (MQ=255)
                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaa                             <  2:1421008/50‑1 (MQ=255)
                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                              <  2:4118931/49‑1 (MQ=255)
                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                              <  2:2435881/49‑1 (MQ=255)
                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                              <  2:3951472/49‑1 (MQ=255)
                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                              <  1:2773222/49‑1 (MQ=255)
                  ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGTTTGCATACCaaa                             >  2:1989014/1‑49 (MQ=255)
                  ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            <  2:1340753/50‑1 (MQ=255)
                  ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            <  2:1813930/50‑1 (MQ=255)
                  ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            >  2:1414541/1‑50 (MQ=255)
                   cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            >  1:1528692/1‑49 (MQ=255)
                   cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                           <  2:1522886/50‑1 (MQ=255)
                    aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACa                          <  1:3264775/50‑1 (MQ=255)
                     tCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                         >  2:4430366/1‑50 (MQ=255)
                     tCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                         >  1:153950/1‑50 (MQ=255)
                      cTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                        >  2:480104/1‑50 (MQ=255)
                      cTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                        >  2:570133/1‑50 (MQ=255)
                      cTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                        >  1:2174365/1‑50 (MQ=255)
                       tGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                       >  1:2081959/1‑50 (MQ=255)
                        gAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                      >  2:3674411/1‑50 (MQ=255)
                        gAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                      >  1:849803/1‑50 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     <  2:1024827/50‑1 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     >  2:1092275/1‑50 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     >  1:3465631/1‑50 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     <  2:3681431/50‑1 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     >  1:1159867/1‑50 (MQ=255)
                          aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    >  2:2463744/1‑50 (MQ=255)
                          aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    >  2:2830425/1‑50 (MQ=255)
                          aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    <  2:4358224/50‑1 (MQ=255)
                          aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    <  1:4560110/50‑1 (MQ=255)
                          aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    <  1:385121/50‑1 (MQ=255)
                          aCTATAAGGATCTATCCTGATTTGGCTAATTGCATACCAAAACAGCttgc                    <  2:3991484/50‑3 (MQ=255)
                           cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                   >  1:3676616/1‑50 (MQ=255)
                            tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAATCAGCTTTCGc                  >  2:3402158/1‑50 (MQ=255)
                            tATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  1:383685/1‑50 (MQ=255)
                             aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  2:3300866/1‑49 (MQ=255)
                             aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  2:2193414/1‑49 (MQ=255)
                             aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  1:4213882/1‑49 (MQ=255)
                             aTCAGGAACTTCCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  2:3363275/1‑49 (MQ=39)
                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  <  2:2010490/48‑1 (MQ=255)
                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCt                 <  2:1510414/49‑1 (MQ=255)
                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                <  1:2428957/50‑1 (MQ=255)
                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                >  2:3370035/1‑50 (MQ=255)
                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTATCGc                  <  2:1470149/48‑1 (MQ=39)
                               cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               >  2:1328388/1‑50 (MQ=255)
                               cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               >  1:1623382/1‑50 (MQ=255)
                               cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               >  1:997073/1‑50 (MQ=255)
                                aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg              >  1:391418/1‑50 (MQ=255)
                                 ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt             <  1:1927826/50‑1 (MQ=255)
                                  gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            >  1:1056871/1‑50 (MQ=255)
                                  gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            <  1:2321774/50‑1 (MQ=255)
                                  gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            >  2:1805416/1‑50 (MQ=255)
                                  gAACTTTCCTGATTTGGCTGATTGCATACCAAAACACCTTTCGCTACGtt            >  1:1028682/1‑50 (MQ=255)
                                   aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg           >  1:4309511/1‑50 (MQ=255)
                                   aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg           <  2:220902/50‑1 (MQ=255)
                                   aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg           <  2:223622/50‑1 (MQ=255)
                                    aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc          <  2:575977/50‑1 (MQ=255)
                                     cTTTCCTGATTTGGCTGATTGCCTACCAAAACAGCTTCCGCTACGTTGCt         >  2:2590533/1‑50 (MQ=255)
                                     cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt         <  1:3833377/50‑1 (MQ=255)
                                     cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt         <  2:3072644/50‑1 (MQ=255)
                                      tttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg        >  1:1739245/1‑50 (MQ=255)
                                      tttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTg        >  1:1876683/1‑50 (MQ=255)
                                       ttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg       >  2:1372987/1‑50 (MQ=255)
                                        tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc      >  2:1818020/1‑50 (MQ=255)
                                         ccTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt     >  2:2953345/1‑50 (MQ=255)
                                          cagATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc      >  1:586740/3‑48 (MQ=39)
                                          cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCt     <  2:2242699/49‑1 (MQ=255)
                                          cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc    >  1:2786107/1‑50 (MQ=255)
                                          cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc    <  1:4184875/50‑1 (MQ=255)
                                           tGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCg   >  2:1473320/1‑50 (MQ=255)
                                               |                                              
CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGT  >  NC_000913/3090235‑3090328

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