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.9% 457.9 / ‑2.3 92intergenic (+89/‑6)yggI/endAconserved protein YggI/DNA‑specific endonuclease I
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (50/41);  minor base A (1/0);  total (51/41)
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

ACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAGC  >  NC_000913/3090296‑3090385
                                             |                                            
aCCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTg                                         >  1:2644222/1‑50 (MQ=255)
aCCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTg                                         >  1:2221779/1‑50 (MQ=255)
 ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACATGGAGAAGTGc                                        >  1:3657674‑M1/1‑49 (MQ=255)
 ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTg                                         <  1:5628364/49‑1 (MQ=255)
 ccAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGc                                        >  1:534567‑M1/1‑49 (MQ=255)
  cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:7078249‑M1/1‑48 (MQ=255)
  cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       <  1:1995103‑M1/50‑3 (MQ=255)
  cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       <  1:5708920‑M1/50‑3 (MQ=255)
  cAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:4601650‑M1/1‑48 (MQ=255)
   aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:1717056‑M1/1‑47 (MQ=255)
   aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:2658240‑M1/1‑47 (MQ=255)
   aaaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcc                                       >  1:6955186‑M1/1‑47 (MQ=255)
    aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                     <  1:5985780‑M1/50‑5 (MQ=255)
    aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                     <  1:1908408‑M1/50‑5 (MQ=255)
    aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                     >  1:2190082‑M1/1‑46 (MQ=255)
    aaaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccca                                     >  1:6682756‑M1/1‑46 (MQ=255)
     aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctac                                    <  1:2174946‑M1/50‑6 (MQ=255)
     aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctac                                    <  1:6751979‑M1/50‑6 (MQ=255)
     aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctac                                    <  1:1657817‑M1/50‑6 (MQ=255)
     aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                     >  1:1297058‑M1/1‑45 (MQ=255)
     aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                     >  1:307016‑M1/1‑45 (MQ=255)
     aaCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGccta                                     >  1:6648843‑M1/1‑45 (MQ=255)
      aCAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctaca                                   <  1:6506845‑M1/50‑7 (MQ=255)
        aGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacact                                 <  1:873051‑M1/50‑9 (MQ=255)
         gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacacta                                <  1:5209711‑M1/50‑10 (MQ=255)
         gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacacta                                <  1:5746464‑M1/50‑10 (MQ=255)
         gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacacta                                <  1:2939587‑M1/50‑10 (MQ=255)
         gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                  >  1:4286275‑M1/1‑41 (MQ=255)
         gCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacac                                  >  1:6018553‑M1/1‑41 (MQ=255)
          cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                               <  1:4142898‑M1/50‑11 (MQ=255)
          cTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactag                               <  1:5167678‑M1/50‑11 (MQ=255)
           tttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              >  1:7117734‑M1/1‑39 (MQ=255)
           tttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              >  1:5993671‑M1/1‑39 (MQ=255)
           tttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              >  1:3016553‑M1/1‑39 (MQ=255)
           tttCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              >  1:6927531‑M1/1‑39 (MQ=255)
             tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcg                             <  1:3085113‑M1/49‑13 (MQ=255)
             tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              <  1:5116068‑M1/48‑12 (MQ=255)
             tCGCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagc                              <  1:1457946‑M1/48‑12 (MQ=255)
               gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                          >  1:2953165‑M1/1‑35 (MQ=255)
               gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                          >  1:4693007‑M1/1‑35 (MQ=255)
               gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                          >  1:101731‑M1/1‑35 (MQ=255)
               gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                          >  1:4504376‑M1/1‑35 (MQ=255)
               gCTACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcggga                          >  1:5428852‑M1/1‑35 (MQ=255)
                cTACGTTGCTGGCTCGTTTTAACACGGAAAAGTGcctacactagcgggat                         >  1:7156266‑M1/1‑34 (MQ=255)
                 tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                        >  1:216546‑M1/1‑33 (MQ=255)
                 tACGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                        >  1:6138998‑M1/1‑33 (MQ=255)
                  aCGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattc                       >  1:1306482‑M1/1‑32 (MQ=255)
                   cGTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggatt                        <  1:4130110‑M1/48‑18 (MQ=255)
                    gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     <  1:7379645‑M1/50‑21 (MQ=255)
                    gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     <  1:5658706‑M1/50‑21 (MQ=255)
                    gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     >  1:6277384‑M1/1‑30 (MQ=255)
                    gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     <  1:989875‑M1/50‑21 (MQ=255)
                    gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     >  1:5938981‑M1/1‑30 (MQ=255)
                    gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     <  1:242048‑M1/50‑21 (MQ=255)
                    gTTGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctt                     <  1:2961248‑M1/50‑21 (MQ=255)
                      tGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                   >  1:1012114‑M1/1‑28 (MQ=255)
                      tGCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                   >  1:549802‑M1/1‑28 (MQ=255)
                       gCTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttt                   <  1:250266‑M1/49‑23 (MQ=255)
                       gCTGGCTCGTTTTAACACGGAGAAGCGcctacactagcgggattcttttt                  >  1:2437684‑M1/1‑27 (MQ=255)
                        cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                  <  1:1434798‑M1/49‑24 (MQ=255)
                        cTGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattcttttt                  <  1:2275573‑M1/49‑24 (MQ=255)
                         tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggcttctttttgt                >  1:2560353‑M1/1‑25 (MQ=255)
                         tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                >  1:5332134‑M1/1‑25 (MQ=255)
                         tGGCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                >  1:5695726‑M1/1‑25 (MQ=255)
                          ggCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                >  1:563245‑M1/1‑24 (MQ=255)
                           gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt               <  1:732560‑M1/49‑27 (MQ=255)
                           gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt               <  1:1949240‑M1/49‑27 (MQ=255)
                           gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgtt               >  1:3130959‑M1/1‑23 (MQ=255)
                           gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                <  1:6160948‑M1/48‑26 (MQ=255)
                           gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                <  1:5504148‑M1/48‑26 (MQ=255)
                           gCTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgt                <  1:4279685‑M1/48‑26 (MQ=255)
                            cTCGTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaa             >  1:2344210‑M1/1‑22 (MQ=255)
                               gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaaccc          >  1:6937670‑M1/1‑19 (MQ=255)
                               gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaaccc          >  1:122598‑M1/1‑19 (MQ=255)
                               gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaaccc          >  1:1005588‑M1/1‑19 (MQ=255)
                               gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaaccc          <  1:5267167‑M1/50‑32 (MQ=255)
                               gTTTTAACACGGAGAAGTGcctacactagcgggattctttttgttaaccc          >  1:393330‑M1/1‑19 (MQ=255)
                                ttttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc         >  1:2617404‑M1/1‑18 (MQ=255)
                                ttttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc         <  1:3348852‑M1/50‑33 (MQ=255)
                                ttttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc         >  1:4984297‑M1/1‑18 (MQ=255)
                                 tttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccc         >  1:286698‑M1/1‑17 (MQ=255)
                                  ttAACACGGAGAAGTGcctacactagcgggattctttttgttaaccccta       >  1:6085937‑M1/1‑16 (MQ=255)
                                  ttAACACGGAGAAGTGcctacactagcgggattctttttgttaaccccta       >  1:5400924‑M1/1‑16 (MQ=255)
                                  ttAACACGGAGAAGTGcctacactagcgggattctttttgttaaccccta       >  1:534798‑M1/1‑16 (MQ=255)
                                  ttAACACGGAGAAGTGcctacactagcgggattctttttgttaacccct        <  1:2767759‑M1/49‑34 (MQ=255)
                                    aaCACGGAGAAGTGcctacactagcgggattctttttgttaacccctacc     >  1:226356‑M1/1‑14 (MQ=255)
                                    aaCACGGAGAAGTGcctacactagcgggattctttttgttaacccctac      <  1:2799960‑M1/49‑36 (MQ=255)
                                    aaCACGGAGAAGTGcctacactagcgggattctttttgttaacccctac      <  1:4080011‑M1/49‑36 (MQ=255)
                                     acacGGAGAAGTGcctacactagcgggattctttttgttaacccctacc     <  1:767281‑M1/49‑37 (MQ=255)
                                     acacGGAGAAGTGcctacactagcgggattctttttgttaacccctacc     <  1:5237009‑M1/49‑37 (MQ=255)
                                      cacGGAGAAGTGcctacactagcgggattctttttgttaacccctacccc   <  1:1525028‑M1/50‑39 (MQ=255)
                                       acGGAGAAGTGcctacactagcgggattctttttgttaacccctacccca  <  1:5777585‑M1/50‑40 (MQ=255)
                                             |                                            
ACCAAAACAGCTTTCGCTACGTTGCTGGCTCGTTTTAACACGGAGTAAGTGATGTACCGTTATTTGTCTATTGCTGCGGTGGTACTGAGC  >  NC_000913/3090296‑3090385

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