Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,222,671 G→C 10.4% intergenic (+127/+246) yohP → / ← dusC uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,222,6710GC10.4% 260.2 / 6.0 111intergenic (+127/+246)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base G (9/84);  new base C (8/2);  total (25/86)
Fisher's exact test for biased strand distribution p-value = 3.86e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 155 total aligned reads displayed.

GCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATC  >  CP009273/2222606‑2222738
                                                                       |                                                                    
gcTCCTGTCATGTGCATTTCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATcgcc                                                                    <  1:44997805/68‑1 (MQ=255)
gcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATcgcc                                                                    <  1:43347757/68‑1 (MQ=255)
gcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATcgcc                                                                    <  1:23909641/68‑1 (MQ=255)
gcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATcgcc                                                                    <  1:48022774/68‑1 (MQ=255)
  tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑‑CGGGGGGGGGCCAAcgccg                                                                   >  1:6282074/1‑67 (MQ=25)
   ccTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATcgccgcc                                                                 <  1:51653182/68‑1 (MQ=255)
     tGTCCTGTGCCTTGCTTCATATACTCTCTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATcgccgccgc                                                               <  1:14931022/68‑1 (MQ=18)
     tGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATcgccgccgc                                                               <  1:38191745/68‑1 (MQ=255)
     tGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATcgccgccgc                                                               <  1:48994421/68‑1 (MQ=255)
     tGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATcgccgccgc                                                               <  1:53753435/68‑1 (MQ=255)
      gTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCgccgccgcc                                                              <  1:23930083/68‑1 (MQ=255)
      gTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCgccgccgcc                                                              <  1:4916283/68‑1 (MQ=255)
        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGcccc                                                            <  1:57461167/68‑1 (MQ=255)
        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGcccc                                                            <  1:3365499/68‑1 (MQ=255)
        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGccc                                                             <  1:24676092/67‑1 (MQ=255)
        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGccc                                                             <  1:56821442/67‑1 (MQ=255)
        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGccc                                                             <  1:10469799/67‑1 (MQ=255)
        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGccc                                                             <  1:40224675/67‑1 (MQ=255)
        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGccc                                                             <  1:49100376/67‑1 (MQ=255)
        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGccc                                                             <  1:43111201/67‑1 (MQ=255)
        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGccc                                                             <  1:14141919/67‑1 (MQ=255)
        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGccc                                                             <  1:13662489/67‑1 (MQ=255)
        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGccc                                                             <  1:6460446/67‑1 (MQ=255)
        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGccc                                                             <  1:8639963/67‑1 (MQ=255)
        cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGccc                                                             <  1:1246604/67‑1 (MQ=255)
        cATGTGCATTGATTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGccc                                                             <  1:56186392/67‑1 (MQ=255)
          tgtgCATTGCTTCATATAATCACTGGCGCAAGGCGCGCGG‑‑GCGGGGGCGGCCCAAcccccccccccc                                                            >  1:47492133/1‑67 (MQ=255)
          tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑‑CGGGGGGGGGCCCAAcccccccccccc                                                            >  1:28407694/1‑67 (MQ=19)
          tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑‑CGGGGGGGGGCCAACcccccccccccc                                                            >  1:9127732/1‑67 (MQ=21)
          tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑GGGGGGGCGGGCCAAAcccccccccccc                                                            >  1:49918469/1‑67 (MQ=21)
          tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑GCCGGGGGGGCCCAAAGCccccccccccc                                                          >  1:20943921/1‑67 (MQ=19)
          tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑GCCGGGGCGGGCCAAAGCccccccccccc                                                          >  1:14506666/1‑67 (MQ=19)
          tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCt                                                          <  1:46204961/68‑1 (MQ=255)
          tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCt                                                          <  1:17814696/68‑1 (MQ=255)
          tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCt                                                          <  1:12730945/68‑1 (MQ=255)
          tgtgCATTGCTGCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCt                                                          <  1:20716281/68‑1 (MQ=255)
             gCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGG‑‑CGGGGGGGGGCCAAccccccccccccccc                                                          >  1:18633335/1‑66 (MQ=16)
             gCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑GCCCGGGGGGCGCCCAAccccccccccccccc                                                          >  1:29936407/1‑66 (MQ=16)
             gCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑GCCGGGGGGGGGCCAAACcccccccccccccc                                                         >  1:41995420/1‑65 (MQ=16)
             gCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑GAGGGGGGGGGCCAAAGCccccccccccacc                                                        >  1:26230534/1‑64 (MQ=16)
             gCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGGGGCCAATCgcccccgcccccacc                                                        >  1:28678889/1‑64 (MQ=18)
               tttGCTTCAAATAATCACTGGCCCAAGCAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAcc                                                    <  1:35007641/67‑1 (MQ=17)
               aTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑GGGGGGGGGGGCCCACGCCCCCGGCCCCCGccc                                                     >  1:5290171/1‑65 (MQ=16)
               aTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGGGCCCAAACCCCccgccccccgcccc                                                    >  1:27631117/1‑67 (MQ=18)
               aTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCACTCGCCGCCGCCCCCTGCAcc                                                    <  1:54058249/68‑1 (MQ=255)
               aTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAcc                                                    <  1:51798147/68‑1 (MQ=255)
               aTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAcc                                                    <  1:36743664/68‑1 (MQ=255)
               aTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAcc                                                    <  1:42743237/68‑1 (MQ=255)
               aTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAcc                                                    <  1:17909468/68‑1 (MQ=255)
                 tGCTTCATATAATCCCTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAcccc                                                  <  1:40769736/68‑1 (MQ=255)
                 tGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAcccc                                                  <  1:15996950/68‑1 (MQ=255)
                 tGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCGGCAcccc                                                  <  1:31031612/68‑1 (MQ=255)
                 tGCTTAATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CTGGGGGCGGCAAATCGCCGCCGCCCCCTGCAcccc                                                  <  1:55734704/68‑1 (MQ=17)
                 tGCTTAATAAAATCACTGGCCCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCAcccc                                                  <  1:18167135/68‑1 (MQ=14)
                  gCTTCATATAATCACTGGCGCAAGGAGCGCGG‑‑CGGGGGGGGGCCAAAACCGCCCCCCCCCGcccccc                                                  >  1:57218695/1‑67 (MQ=12)
                  gCTTCATATAATCACTGGCGCAAGGAGCGC‑‑T‑‑CAGGGGGGGGCCCATCGACCCCCCCCCCAGGAcccc                                                  >  1:7128551/1‑67 (MQ=17)
                  gCTTCATATAATCACTGGCGCAAGGAGCGC‑‑T‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAccccc                                                 >  1:38410502/1‑68 (MQ=255)
                  gCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑GCGGGGGGGGCCCAAAACCCCCCCCCCCCGcccccc                                                  >  1:25652955/1‑67 (MQ=12)
                    ttCATATAATCACTGGCGCACGGAGCGC‑‑G‑‑CAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCgg                                               <  1:18945715/68‑1 (MQ=17)
                    ttCATATAATCACTGGCGCAAGGAGCGCGG‑‑GGGGGGGGGGCCAAACCCCCCCCCCCCCCCAccccccc                                               >  1:1354996/1‑66 (MQ=12)
                    ttCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCgg                                               <  1:32571285/68‑1 (MQ=255)
                      cATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGc                                             <  1:32102786/68‑1 (MQ=255)
                         ataATCACTGGCGCAAGGAGCGCGGG‑‑GGGGGGGGGGCCAACCCCCCCCCCCCCCCCCCCCCCCggggg                                            >  1:56472388/1‑66 (MQ=9)
                         ataATCACTGGCGCAAGGAGCGCGGG‑‑GGGGGGGGGGCCAACCCCCCCCCCCCCCCCCCCCCCCGGGcg                                            >  1:16675692/1‑67 (MQ=9)
                         ataATCACTGGCGCAAGGAGCGCGG‑‑GGGGGGGGGGCCAAAGCCCCCCCCCCCCCCCCCCCCCggggcg                                           >  1:39260314/1‑65 (MQ=9)
                         ataATCACTGGCGCAAGGAGCGCGG‑‑CGGGGGGGGGCCCAAGCCCCCCCCCCCCCCCCCCCCCCggccg                                            >  1:37465600/1‑65 (MQ=9)
                         ataATCACTGGCGCAAGGAGCGCGG‑‑CGGGGGGGGGCCAAACCCCCCCCCCCCCCCACCCCCCGGcgca                                          >  1:13044263/1‑65 (MQ=9)
                         ataATCACTGGCGCAAGGAGCGCGG‑‑CGGGGGGGGGCCAAACCCCCCCCCCCCAGCCCCCCCCGGcgcc                                          >  1:9681212/1‑65 (MQ=9)
                         ataATCACTGGCGCAAGGAGCGC‑‑G‑‑GCGGGGGGGGGCCCGCGCCGCCCCCCCCCGCCGCCCCCgggggc                                           >  1:46986302/1‑65 (MQ=9)
                         ataATCACTGGCGCAAGGAGCGC‑‑G‑‑GCGGGGGGGGGCCAGAACCCCCCCCCCCCGCCCCCCCgcgcggg                                          >  1:32123218/1‑65 (MQ=9)
                            aTCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCTCCGGGCTCTGGc                                       >  1:19991321/1‑68 (MQ=18)
                                                                 ccAAACCCCTCCGCCCCCGGCACCCCCGGGCTATGGCGAACAACATCGCCGCTTCGCGGTGCCCTCgg        <  1:13142580/68‑1 (MQ=11)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg        <  1:416818/67‑1 (MQ=31)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg        <  1:9440284/67‑1 (MQ=31)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg        <  1:36640120/67‑1 (MQ=31)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg        <  1:6155064/67‑1 (MQ=31)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg        <  1:35502671/67‑1 (MQ=31)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg        <  1:33130932/67‑1 (MQ=31)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg        <  1:53752248/67‑1 (MQ=31)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg        <  1:24493621/67‑1 (MQ=31)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg        <  1:2724906/67‑1 (MQ=31)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg        <  1:8032660/67‑1 (MQ=31)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg        <  1:29460679/67‑1 (MQ=31)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg        <  1:27580573/67‑1 (MQ=31)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc       <  1:1241044/68‑1 (MQ=33)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc       <  1:44800177/68‑1 (MQ=33)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc       <  1:45396358/68‑1 (MQ=33)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc       <  1:30395434/68‑1 (MQ=33)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc       <  1:44534212/68‑1 (MQ=33)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc       <  1:12876634/68‑1 (MQ=33)
                                                                  cAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc       <  1:10273968/68‑1 (MQ=33)
                                                                  cAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg        <  1:50818807/67‑1 (MQ=14)
                                                                  cAATCGCCGCCGCCCCCGGCAACCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc       <  1:54421804/68‑1 (MQ=255)
                                                                  cAATCCCGGCCGCCCCCGGCACCCACTGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc       <  1:27388975/68‑1 (MQ=11)
                                                                      cgccgccgcCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAt   <  1:19637571/68‑1 (MQ=35)
                                                                      cgccgccgcCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAt   <  1:740269/68‑1 (MQ=35)
                                                                      cgccgccgcCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAt   <  1:40994360/68‑1 (MQ=35)
                                                                       gccgccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAt   >  1:32516254/1‑67 (MQ=35)
                                                                       gccgccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAt   >  1:45078250/1‑67 (MQ=35)
                                                                       gccgccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTAt   >  1:45656766/1‑67 (MQ=35)
                                                                       |                                                                    
GCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑‑G‑‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATC  >  CP009273/2222606‑2222738

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