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.99.3% 685.6 / ‑1.9 141intergenic (+89/‑6)yggI/endAconserved protein YggI/DNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (73/66);  total (74/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

ACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAG  >  NC_000913/3090296‑3090384
                                             |                                           
aCCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTg                                        >  1:348233/1‑50 (MQ=255)
aCCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTg                                        >  1:2549538/1‑50 (MQ=255)
 ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGc                                       >  1:1332552‑M1/1‑49 (MQ=255)
  cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                      <  1:3269919‑M1/50‑3 (MQ=255)
  cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                      <  1:4887086‑M1/50‑3 (MQ=255)
  cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                      >  1:2584897‑M1/1‑48 (MQ=255)
  cAAAACAGCTTTCGCTACGTTGCTGGATCGTTTTAACACGGAGAAGTGcc                                      >  1:1546178‑M1/1‑48 (MQ=255)
    aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                    <  1:4285754‑M1/50‑5 (MQ=255)
    aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                    >  2:1552251‑M1/1‑46 (MQ=255)
    aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGATAAGTGccta                                    >  1:1747545‑M1/1‑46 (MQ=255)
      aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctaca                                  <  1:1173620‑M1/50‑7 (MQ=255)
      aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctaca                                  <  1:4216318‑M1/50‑7 (MQ=255)
      aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctaca                                  <  2:1720052‑M1/50‑7 (MQ=255)
      aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctaca                                  >  2:2680833‑M1/1‑44 (MQ=255)
      aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctaca                                  <  2:3261415‑M1/50‑7 (MQ=255)
      aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctaca                                  <  2:3272350‑M1/50‑7 (MQ=255)
       cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                 >  2:4792338‑M1/1‑43 (MQ=255)
       cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                 >  2:642504‑M1/1‑43 (MQ=255)
       cAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                 <  1:4002846‑M1/50‑8 (MQ=255)
        aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccttcact                                <  2:112206‑M1/50‑9 (MQ=255)
        aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacact                                <  1:4240455‑M1/50‑9 (MQ=255)
        aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                 >  2:2730292‑M1/1‑42 (MQ=255)
         gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacacta                               <  2:3499668‑M1/50‑10 (MQ=255)
         gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacacta                               <  2:4571669‑M1/50‑10 (MQ=255)
          cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                              <  1:602069‑M1/50‑11 (MQ=255)
          cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                              <  2:3718238‑M1/50‑11 (MQ=255)
          cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                              <  2:5125920‑M1/50‑11 (MQ=255)
          cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacaatag                              <  2:551155‑M1/50‑11 (MQ=255)
          cTTTCGCTAAGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                              >  2:3977909‑M1/1‑40 (MQ=255)
           tttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                             >  2:216628‑M1/1‑39 (MQ=255)
            ttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcg                            >  2:2370993‑M1/1‑38 (MQ=255)
             tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgg                           >  2:667174‑M1/1‑37 (MQ=255)
             tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcg                            <  2:3970155‑M1/49‑13 (MQ=255)
             tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcg                            <  2:2882496‑M1/49‑13 (MQ=255)
             tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcg                            <  1:4557032‑M1/49‑13 (MQ=255)
             tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                             <  1:2844239‑M1/48‑12 (MQ=255)
               gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                         >  2:5537624‑M1/1‑35 (MQ=255)
               gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                         >  1:2672257‑M1/1‑35 (MQ=255)
               gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                         <  1:3780036‑M1/50‑16 (MQ=255)
               gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                         >  1:4273280‑M1/1‑35 (MQ=255)
               gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggg                          <  2:4838505‑M1/49‑15 (MQ=255)
                cTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggat                        >  2:2082858‑M1/1‑34 (MQ=255)
                cTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggat                        <  2:1866632‑M1/50‑17 (MQ=255)
                 tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                       >  2:4315249‑M1/1‑33 (MQ=255)
                 tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                       >  2:4903883‑M1/1‑33 (MQ=255)
                 tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                       >  2:4006231‑M1/1‑33 (MQ=255)
                 tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                       >  1:352939‑M1/1‑33 (MQ=255)
                  aCGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattc                      >  2:5631946‑M1/1‑32 (MQ=255)
                   cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattct                     <  2:2473546‑M1/50‑20 (MQ=255)
                   cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattc                      <  1:4429787‑M1/49‑19 (MQ=255)
                   cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                       <  1:3973607‑M1/48‑18 (MQ=255)
                    gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                    >  1:5155103‑M1/1‑30 (MQ=255)
                     ttGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttt                   >  1:2691199‑M1/1‑29 (MQ=255)
                     ttGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttt                   <  1:1338601‑M1/50‑22 (MQ=255)
                      tGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                  >  2:3857393‑M1/1‑28 (MQ=255)
                      tGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                  >  2:4530705‑M1/1‑28 (MQ=255)
                      tGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                  >  2:3014289‑M1/1‑28 (MQ=255)
                       gCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                 >  2:1515890‑M1/1‑27 (MQ=255)
                       gCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                  <  1:4557383‑M1/49‑23 (MQ=255)
                       gCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                  <  2:2879887‑M1/49‑23 (MQ=255)
                       gCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                  <  2:5147544‑M1/49‑23 (MQ=255)
                       gCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                  <  1:5765545‑M1/49‑23 (MQ=255)
                       gCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                  <  1:5625120‑M1/49‑23 (MQ=255)
                        cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttg                <  2:3496563‑M1/50‑25 (MQ=255)
                        cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttg                >  2:5602356‑M1/1‑26 (MQ=255)
                        cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                 <  1:4157851‑M1/49‑24 (MQ=255)
                         tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt               >  2:2471716‑M1/1‑25 (MQ=255)
                         tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt               >  2:4075706‑M1/1‑25 (MQ=255)
                         tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt               >  2:1048373‑M1/1‑25 (MQ=255)
                         tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt               >  1:5160356‑M1/1‑25 (MQ=255)
                          ggCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt               >  2:4818464‑M1/1‑24 (MQ=255)
                          ggCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt               >  2:1313914‑M1/1‑24 (MQ=255)
                          ggCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt               >  1:1304875‑M1/1‑24 (MQ=255)
                           gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtta             <  1:4030067‑M1/50‑28 (MQ=255)
                           gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt              >  1:1743605‑M1/1‑23 (MQ=255)
                           gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt              <  1:4903526‑M1/49‑27 (MQ=255)
                           gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt              >  1:4207225‑M1/1‑23 (MQ=255)
                            cTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaa            <  2:4879190‑M1/50‑29 (MQ=255)
                            cTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaa            <  1:4184149‑M1/50‑29 (MQ=255)
                             tCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaac           >  1:2439561‑M1/1‑21 (MQ=255)
                              cGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaacc          >  1:5365592‑M1/1‑20 (MQ=255)
                               gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaaccc         >  2:708950‑M1/1‑19 (MQ=255)
                               gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaaccc         >  1:3709090‑M1/1‑19 (MQ=255)
                               gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaac           <  2:2564466‑M1/48‑30 (MQ=255)
                                ttttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc        >  1:1147387‑M1/1‑18 (MQ=255)
                                ttttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc        <  2:5178021‑M1/50‑33 (MQ=255)
                                 tttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc        >  2:625245‑M1/1‑17 (MQ=255)
                                 tttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc        >  1:2838842‑M1/1‑17 (MQ=255)
                                  ttAACACGGAGAAGTGcctacactagcgggattctttttgttaaccccta      >  2:5729446‑M1/1‑16 (MQ=255)
                                  ttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccct       <  2:2052557‑M1/49‑34 (MQ=255)
                                   tAACACGGAGAAGTGcctacactagcgggattctttttgttaacccctac     >  2:3504184‑M1/1‑15 (MQ=255)
                                   tAACACGGAGAAGTGcctacactagcgggattctttttgttaaccccta      <  1:671550‑M1/49‑35 (MQ=255)
                                   tAACACGGAGAAGTGcctacactagcgggattctttttgttaaccccta      <  1:3222370‑M1/49‑35 (MQ=255)
                                   tAACACGGAGAAGTGcctacactagcgggattctttttgttaaccccta      <  1:533783‑M1/49‑35 (MQ=255)
                                    aaCACGGAGAAGTGcctacactagcgggcttctttttgttaacccctacc    >  1:5244136‑M1/1‑14 (MQ=255)
                                    aaCACGGAGAAGTGcctacactagcgggattctttttgttaacccctac     <  1:5533983‑M1/49‑36 (MQ=255)
                                    aaCACGGGGAAGTGcctacactagcgggattctttttgttaacccctac     <  1:5705082‑M1/49‑36 (MQ=255)
                                     acacGGAGAAGTGcctacactagcgggattctttttgttaacccctacc    <  1:817800‑M1/49‑37 (MQ=255)
                                      cacGGAGAAGTGcctacactagcgggattctttttgttaacccctacccc  <  2:1662282‑M1/50‑39 (MQ=255)
                                      cacGGAGAAGTGcctacactagcgggattctttttgttaacccctacccc  <  1:2912948‑M1/50‑39 (MQ=255)
                                             |                                           
ACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAG  >  NC_000913/3090296‑3090384

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