Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,019,857 C→T intergenic (+56/‑28) 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,8570CT98.3% 158.1 / NA 59intergenic (+56/‑28)b3836/b3838hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (25/33);  total (26/33)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGCCGCCTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGC  >  NC_000913/4019787‑4019919
                                                                      |                                                               
aaaGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCCTCATCGGCCtcgtcgttt                                                            >  2:312295/1‑73 (MQ=255)
      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:3792348/76‑1 (MQ=255)
       aGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGcc                                                   >  2:3671578/1‑76 (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:3838523/76‑1 (MQ=255)
            tAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAAc                                              <  1:298040/76‑1 (MQ=255)
               tCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGAc                                           <  2:4858978/76‑1 (MQ=255)
               tCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGAc                                           >  2:636122/1‑76 (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:4959346/1‑76 (MQ=255)
                           atCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgg                              >  1:4944210/1‑76 (MQ=255)
                           atCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGTCG‑‑TTCTGGGGCCGCAACGACTGCCTGTGGCgg                              >  1:318160/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)
                                                                      |                                                               
AAAGAGCAGGTGTAATCCGTGTTTGATATCGGTTTTAGCGAACTGCTATTGGTGTTCATCATCGGCCTCGCCGCCTTCTGGGGCCGCAACGACTGCCTGTGGCGGTAAAAACGGTAGCGGGCTGGATTCGCGC  >  NC_000913/4019787‑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: