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,2820CT95.3% 620.8 / 11.5 170intergenic (+30/‑65)yggI/endAconserved protein YggI/DNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base T (72/90);  total (77/93)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 170 total aligned reads displayed.

CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGT  >  NC_000913/3090235‑3090328
                                               |                                              
cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              <  2:2416221/50‑1 (MQ=255)
cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              <  2:2894018/50‑1 (MQ=255)
cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              <  1:909014/50‑1 (MQ=255)
cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              <  1:5035561/50‑1 (MQ=255)
cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              <  2:4356369/50‑1 (MQ=255)
cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              <  1:2755300/50‑1 (MQ=255)
cAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTg                                              >  1:2332532/1‑50 (MQ=255)
 aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                             >  1:1238449/1‑50 (MQ=255)
 aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                             >  2:190512/1‑50 (MQ=255)
 aGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTgg                                             <  1:2554279/50‑1 (MQ=255)
  gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                            <  1:4266218/50‑1 (MQ=255)
  gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                            <  1:3656055/50‑1 (MQ=255)
  gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGc                                            <  1:2872478/50‑1 (MQ=255)
  gCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGc                                            <  2:2061421/50‑1 (MQ=255)
   cTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           >  2:4236959/1‑50 (MQ=255)
    tGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           >  1:542131/1‑49 (MQ=255)
    tGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCt                                           >  2:1026873/1‑49 (MQ=255)
     ggTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGa                                         >  2:4492629/1‑50 (MQ=255)
      gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTg                                          <  2:4394877/48‑1 (MQ=255)
      gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGa                                         <  1:4591906/49‑1 (MQ=255)
      gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                        >  1:3392199/1‑50 (MQ=255)
      gTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAt                                        <  1:2213367/50‑1 (MQ=255)
      gTTGCGAAATAACCAGCTGAGCTATCAGGAACTTTCCTGATTTGGCTGAt                                        <  2:5169606/50‑1 (MQ=255)
       ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       <  2:876914/50‑1 (MQ=255)
       ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       <  2:282952/50‑1 (MQ=255)
       ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       <  1:2161560/50‑1 (MQ=255)
       ttGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGAtt                                       >  1:2504919/1‑50 (MQ=255)
        tGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                      >  1:4117892/1‑50 (MQ=255)
         gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                      <  2:5067042/49‑1 (MQ=255)
         gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                      <  2:3215621/49‑1 (MQ=255)
         gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                      >  2:2947588/1‑49 (MQ=255)
         gCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTg                                      <  1:2813196/49‑1 (MQ=255)
          cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGc                                     <  2:4582513/49‑1 (MQ=255)
          cGAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    >  1:4388402/1‑50 (MQ=255)
           gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    >  2:2377597/1‑49 (MQ=255)
           gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCa                                    >  1:1272894/1‑49 (MQ=255)
           gAAATAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCAt                                   <  1:4163422/50‑1 (MQ=255)
           gAAATAACCATCTGAACTATCAGGAACTTTCCTGAATTGGCTGATTGCAt                                   <  1:556301/50‑1 (MQ=255)
            aaaTAGCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCGTa                                  <  2:4126689/50‑1 (MQ=255)
            aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  >  2:5026463/1‑50 (MQ=255)
            aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  <  2:1650005/50‑1 (MQ=255)
            aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  >  2:4664179/1‑50 (MQ=255)
            aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  >  2:4648688/1‑50 (MQ=255)
            aaaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATa                                  >  2:3254887/1‑50 (MQ=255)
             aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 <  1:2149513/50‑1 (MQ=255)
             aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 >  2:1557871/1‑50 (MQ=255)
             aaTAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATAc                                 <  1:4737541/50‑1 (MQ=255)
              aTAACCATCTGAACTATCAGGAGCTTTCCTGATTTGGCTGATTGCATAcc                                <  1:4442478/50‑1 (MQ=255)
               tAACCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                               <  1:1879261/50‑1 (MQ=255)
                aaCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                               >  2:476560/1‑49 (MQ=255)
                aaCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCa                               >  2:5207046/1‑49 (MQ=255)
                aaCCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCa                               >  1:4503466/1‑49 (MQ=255)
                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                              >  1:3861937/1‑49 (MQ=255)
                 aCCATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaa                              >  1:2830135/1‑49 (MQ=255)
                  ccATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            <  2:2799917/50‑1 (MQ=255)
                   cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            >  2:564178/1‑49 (MQ=255)
                   cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                           <  2:316818/50‑1 (MQ=255)
                   cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                           <  1:1746920/50‑1 (MQ=255)
                   cATCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAAc                           <  1:4054547/50‑1 (MQ=255)
                    aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            >  2:946728/1‑48 (MQ=255)
                    aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCaaaa                            >  1:3208135/1‑48 (MQ=255)
                    aTCTGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACa                          <  2:2327847/50‑1 (MQ=255)
                       tGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                       >  2:2166910/1‑50 (MQ=255)
                       tGAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCt                       >  2:3475879/1‑50 (MQ=255)
                        gAACTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                      >  2:5074300/1‑50 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAg                         <  1:1974081/46‑1 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGc                        <  1:4184668/47‑1 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     <  2:2250194/50‑1 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCttt                     >  1:4028848/1‑50 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                      <  1:3675153/49‑1 (MQ=255)
                         aaCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCtt                      <  2:4051058/49‑1 (MQ=255)
                          aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    <  2:3014081/50‑1 (MQ=255)
                          aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    >  2:269588/1‑50 (MQ=255)
                          aCTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTc                    <  2:716990/50‑1 (MQ=255)
                           cTATCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCg                   <  1:2335091/50‑1 (MQ=255)
                             aTCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  >  2:4486629/1‑49 (MQ=255)
                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGc                  <  1:1624721/48‑1 (MQ=255)
                              tCAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTa                <  1:3322660/50‑1 (MQ=255)
                               cAGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               >  1:4367973/1‑50 (MQ=255)
                                aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               <  1:1847342/49‑1 (MQ=255)
                                aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTAc               <  1:852710/49‑1 (MQ=255)
                                aGGAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACg              >  1:2041069/1‑50 (MQ=255)
                                 ggAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGt             <  2:523527/50‑1 (MQ=255)
                                  gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            <  2:2047895/50‑1 (MQ=255)
                                  gAACTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGtt            <  1:2704801/50‑1 (MQ=255)
                                   aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg           <  2:4016613/50‑1 (MQ=255)
                                   aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg           <  1:3912406/50‑1 (MQ=255)
                                   aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg           <  2:1434419/50‑1 (MQ=255)
                                   aaCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTg           >  2:2350717/1‑50 (MQ=255)
                                    aCTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGc          >  2:102400/1‑50 (MQ=255)
                                     cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt         >  2:670550/1‑50 (MQ=255)
                                     cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt         <  1:1212680/50‑1 (MQ=255)
                                     cTTTCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCt         <  2:1203583/50‑1 (MQ=255)
                                       ttCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTgg       >  1:1344570/1‑50 (MQ=255)
                                        tCCTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGc      >  2:1907787/1‑50 (MQ=255)
                                          cTGATTTGGCTTATTGCATACCAAAACAGCTTTGGCTACGTTGCTGGCTc    <  2:435394/50‑1 (MQ=255)
                                          cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc    <  1:94155/50‑1 (MQ=255)
                                          cTGATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTc    <  2:1188129/50‑1 (MQ=255)
                                            gATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt  >  1:3496786/1‑50 (MQ=255)
                                            gATTTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGt  >  2:177459/1‑50 (MQ=255)
                                               |                                              
CAGCTGGTTGCGAAATAACCATCTGAACTATCAGGAACTTTCCTGATCTGGCTGATTGCATACCAAAACAGCTTTCGCTACGTTGCTGGCTCGT  >  NC_000913/3090235‑3090328

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