Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,019,860 Δ2 bp intergenic (+59/‑24) b3836 → / → b3838 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,019,8600C.98.3% 258.6 / ‑4.1 59intergenic (+59/‑25)b3836/b3838hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (26/32);  minor base G (0/1);  total (26/33)
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
*NC_0009134,019,8610C.98.3% 257.9 / ‑3.6 59intergenic (+60/‑24)b3836/b3838hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (26/32);  total (26/33)
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

CAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGCCGCCTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGC  >  NC_000913/4019793‑4019919
                                                                   ||                                                           
cAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGc                                                    >  1:2132888/1‑76 (MQ=255)
cAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGc                                                    >  2:5363908/1‑76 (MQ=255)
cAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGc                                                    <  2:5174520/76‑1 (MQ=255)
cAGGTGTAATCCGTGTTTGATATCGGGTTTAGCGCACTGCTATTGGTGTTCATCCTCGGCCTCGTCG‑‑TTCTGGGGc                                                    >  2:5274854/1‑76 (MQ=255)
cAGGTGTAAGCCGGGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGc                                                    <  2:2161970/76‑1 (MQ=255)
 aGGTGTAGTCCGTGTTGGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGcc                                                   <  2:2130199/76‑1 (MQ=255)
 aGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGcc                                                   >  2:3671578/1‑76 (MQ=255)
 aGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGcc                                                   <  2:3792348/76‑1 (MQ=255)
  ggTGTAACCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCg                                                  <  2:3182830/76‑1 (MQ=255)
   gtgtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGc                                                 >  1:1435068/1‑76 (MQ=255)
   gtgtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGc                                                 >  2:2267447/1‑76 (MQ=255)
   gtgtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGc                                                 >  2:4823168/1‑76 (MQ=255)
     gtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCaa                                               <  1:4853023/76‑1 (MQ=255)
     gtAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCaa                                               >  1:3372847/1‑76 (MQ=255)
      tAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAAc                                              <  1:298040/76‑1 (MQ=255)
      tAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAAc                                              <  1:3838523/76‑1 (MQ=255)
         tCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGAc                                           >  2:636122/1‑76 (MQ=255)
         tCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGAc                                           <  2:4858978/76‑1 (MQ=255)
         tCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCAACATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGAc                                           <  1:2555780/76‑1 (MQ=255)
          ccGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACt                                          >  2:3608668/1‑76 (MQ=255)
          ccGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACt                                          >  1:1715880/1‑76 (MQ=255)
          ccGTGTTTGATATCGGTTTTAGCGAACTGCTATTAGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACt                                          <  1:18821/76‑1 (MQ=255)
           cGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTg                                         >  2:2421858/1‑76 (MQ=255)
             tgtTTGATATCGGTGTTAGAGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGCGCCGCAACGACTGcc                                       >  1:1464259/1‑76 (MQ=255)
              gtTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCt                                      <  1:4537235/76‑1 (MQ=255)
                  gATATCGGTTTTAGCGATCTGTGATGGGGGTTCATCCCCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTgg                                  <  2:4482844/76‑1 (MQ=255)
                  gATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTgg                                  >  2:3130510/1‑76 (MQ=255)
                  gATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTgg                                  >  1:2852725/1‑76 (MQ=255)
                  gATATCGGTTTGAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTgg                                  <  2:2839809/76‑1 (MQ=255)
                    tttCGGTTTTAGCGGACTGCTATGGGTGTTCATCATCGGCTTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCg                               <  2:600621/74‑1 (MQ=255)
                    tatCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCg                               >  2:218415/1‑76 (MQ=255)
                     atCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgg                              >  1:4944210/1‑76 (MQ=255)
                     atCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgg                              >  1:318160/1‑76 (MQ=255)
                     atCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgg                              >  1:4959346/1‑76 (MQ=255)
                        ggTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgggga                           >  2:3983744/1‑73 (MQ=255)
                        ggTTTTAGCGAACTGCTATTGGTGTTCAGCATGGGCCTCGTCG‑‑TTCGGGGGCCGCAACGACTGCCTGTGGCGGTaa                           <  2:4207792/76‑1 (MQ=255)
                            ttAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAAc                       <  1:1189127/76‑1 (MQ=255)
                             tAGCGACCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACg                      <  2:1382968/76‑1 (MQ=255)
                             tAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACg                      <  2:4551315/76‑1 (MQ=255)
                              aGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGGAAAAACg                      >  2:3400098/1‑76 (MQ=255)
                                cGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTa                   <  1:984761/76‑1 (MQ=255)
                                 gaccTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAg                  <  1:1234255/73‑1 (MQ=255)
                                 gAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGGAAAACGGTAg                  >  1:251592/1‑76 (MQ=255)
                                 gAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGGAAAAAGGGAg                  >  1:5633330/1‑76 (MQ=255)
                                  aaCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGc                 <  2:4944209/76‑1 (MQ=255)
                                  aaCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGc                 <  2:1102278/76‑1 (MQ=255)
                                   aCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGGGAAAACAGTAACg                >  2:818831/1‑76 (MQ=255)
                                   aCTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGAGAAAACGGGAGCg                >  1:2572996/1‑76 (MQ=255)
                                       cTATTGTGGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCt            <  1:433331/76‑1 (MQ=255)
                                       cTATTGGTTTTCATCACCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCt            <  2:4190530/76‑1 (MQ=255)
                                        tgttGGTGTTCAGCATCGGCCTCGTCG‑‑TGCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTg           <  2:5066464/74‑1 (MQ=255)
                                        tATTGGTGTTCATCATGGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTg           <  2:1911750/76‑1 (MQ=255)
                                        tATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTg           <  1:3130510/76‑1 (MQ=255)
                                           tGGTGTTCATCAGCGGCCTCGTCG‑‑TTCGGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGAt        <  2:5401464/76‑1 (MQ=255)
                                              tgtTCATCATCGGCCTCGTCGGC‑‑CTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcg     <  2:1151220/76‑1 (MQ=255)
                                               tttCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgc    <  1:768610/75‑1 (MQ=255)
                                               gtTCTTCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgc    <  2:4570495/76‑1 (MQ=255)
                                                 tcatcaTCGGCCTTGCCG‑‑TTCTGCGGCCGCAACGACTGCCTGTGGCGGGAAAAAAGGGCGCAGGGTGGTGTcgcgg  >  1:5718328/1‑75 (MQ=255)
                                                 tcatcaTCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTcgcgc  <  1:2152868/76‑1 (MQ=255)
                                                                   ||                                                           
CAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGCCGCCTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGC  >  NC_000913/4019793‑4019919

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: