Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,090,341 Δ1 bp intergenic (+89/‑6) 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,3410T.98.3% 288.1 / ‑1.5 60intergenic (+89/‑6)yggI/endAconserved protein YggI/DNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (31/28);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAGC  >  NC_000913/3090295‑3090385
                                              |                                            
tACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTagt                                           <  1:1828665/50‑3 (MQ=255)
  ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGc                                        <  1:4911777‑M1/50‑2 (MQ=255)
  ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGc                                        >  1:396796‑M1/1‑49 (MQ=255)
   cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:965067‑M1/1‑48 (MQ=255)
   cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:5316815‑M1/1‑48 (MQ=255)
   cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:3836861‑M1/1‑48 (MQ=255)
    aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcct                                      <  1:3757895‑M1/50‑4 (MQ=255)
    aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcct                                      <  1:4658577‑M1/50‑4 (MQ=255)
    aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:3408019‑M1/1‑47 (MQ=255)
        cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                  <  1:43535‑M1/50‑8 (MQ=255)
        cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                  <  1:5022048‑M1/50‑8 (MQ=255)
        cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                  <  1:4657353‑M1/50‑8 (MQ=255)
        cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                  >  1:4969641‑M1/1‑43 (MQ=255)
         aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacact                                 <  1:198720‑M1/50‑9 (MQ=255)
         aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacact                                 <  1:3938693‑M1/50‑9 (MQ=255)
         aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                  >  1:3662200‑M1/1‑42 (MQ=255)
         aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                  >  1:3159688‑M1/1‑42 (MQ=255)
          gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacacta                                <  1:5249048‑M1/50‑10 (MQ=255)
          gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                  >  1:663016‑M1/1‑41 (MQ=255)
           cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                               >  1:5359298‑M1/1‑40 (MQ=255)
            tttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              >  1:501795‑M1/1‑39 (MQ=255)
            tttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              >  1:1955465‑M1/1‑39 (MQ=255)
              tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                            >  1:3361247‑M1/1‑37 (MQ=255)
               cGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                            <  1:2747638‑M1/49‑14 (MQ=255)
                gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                          >  1:2746159‑M1/1‑35 (MQ=255)
                gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                          >  1:31519‑M1/1‑35 (MQ=255)
                gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggg                           <  1:1812375‑M1/49‑15 (MQ=255)
                gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggg                           <  1:2159592‑M1/49‑15 (MQ=255)
                   aCGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattc                       >  1:4253819‑M1/1‑32 (MQ=255)
                    cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattct                      <  1:3011685‑M1/50‑20 (MQ=255)
                    cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattct                      <  1:871292‑M1/50‑20 (MQ=255)
                    cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattct                      >  1:4036884‑M1/1‑31 (MQ=255)
                     gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     <  1:3673905‑M1/50‑21 (MQ=255)
                     gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     >  1:4521774‑M1/1‑30 (MQ=255)
                      ttGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttt                    <  1:1942158‑M1/50‑22 (MQ=255)
                      ttGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttt                    <  1:581942‑M1/50‑22 (MQ=255)
                       tGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                   >  1:1078937‑M1/1‑28 (MQ=255)
                       tGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                   >  1:2553178‑M1/1‑28 (MQ=255)
                         cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                  <  1:2473787‑M1/49‑24 (MQ=255)
                          tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                >  1:3969669‑M1/1‑25 (MQ=255)
                           ggCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                >  1:917800‑M1/1‑24 (MQ=255)
                           ggCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                >  1:4576226‑M1/1‑24 (MQ=255)
                            gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtta              <  1:1793412‑M1/50‑28 (MQ=255)
                            gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt               <  1:3931763‑M1/49‑27 (MQ=255)
                            gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt               <  1:3409142‑M1/49‑27 (MQ=255)
                             cTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaa             >  1:4845718‑M1/1‑22 (MQ=255)
                              tCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaac            >  1:244182‑M1/1‑21 (MQ=255)
                               cGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaacc           >  1:5534915‑M1/1‑20 (MQ=255)
                               cGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaacc           >  1:3677831‑M1/1‑20 (MQ=255)
                                gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaaccc          <  1:950213‑M1/50‑32 (MQ=255)
                                   ttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccct        <  1:2271233‑M1/49‑34 (MQ=255)
                                   ttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccct        <  1:2161164‑M1/49‑34 (MQ=255)
                                   ttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccct        <  1:2634065‑M1/49‑34 (MQ=255)
                                   ttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccct        <  1:334095‑M1/49‑34 (MQ=255)
                                     aaCACGGAGAAGTGcctacactagcgggattctttttgttaacccctacc     >  1:1774607‑M1/1‑14 (MQ=255)
                                      acacGGAGAAGTGcctacactagcgggattctttttgttaacccctaccc    <  1:1910793‑M1/50‑38 (MQ=255)
                                      acacGGAGAAGTGcctacactagcgggattctttttgttaacccctacc     <  1:4957035‑M1/49‑37 (MQ=255)
                                        acGGAGAAGTGcctacactagcgggattctttttgttaacccctacccca  >  1:1341532‑M1/1‑11 (MQ=255)
                                        acGGAGAAGTGcctacactagcgggattctttttgttaacccctacccca  >  1:5385508‑M1/1‑11 (MQ=255)
                                        acGGAGAAGTGcctacactagcgggattctttttgttaacccctacccca  >  1:1112157‑M1/1‑11 (MQ=255)
                                              |                                            
TACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAGC  >  NC_000913/3090295‑3090385

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