Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,175,949 +G intergenic (+82/‑76) nusG → / → rplK NusG/50S ribosomal subunit protein L11

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,175,9471.G96.9% 255.6 / ‑0.5 99intergenic (+80/‑78)nusG/rplKNusG/50S ribosomal subunit protein L11
Reads supporting (aligned to +/- strand):  ref base . (1/1);  major base G (45/48);  minor base A (1/2);  total (48/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 106 total aligned reads displayed.

GCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGA  >  NC_000913/4175872‑4176017
                                                                            |                                                                        
gCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTAt                                                                           >  1:4436356/1‑76 (MQ=255)
  gATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgg                                                                       >  1:4253387/1‑76 (MQ=255)
  gATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgg                                                                       >  1:4409413/1‑76 (MQ=255)
    tCAAAAAATCGGCGATTTAATCGTTGCACAAGTCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgggc                                                                     <  1:476865/76‑3 (MQ=255)
    tCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgggc                                                                     >  2:2601226/1‑74 (MQ=255)
     cAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATgggcc                                                                    >  1:4112313/1‑73 (MQ=255)
     cAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGAGCGATTGGAATACAATTTCGCGACTTTTGTTTTTATgggcc                                                                    >  2:4338357/1‑73 (MQ=255)
        aaaaGCGGCGATTTAATCGTTGCACAAGGCGGGAGATTGGAATACAATTTCGCGCCGTTTGTTTTTATGGGCCTTg                                                                 <  1:3935474/76‑1 (MQ=255)
         aaaGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGc                                                                >  2:5740352/1‑76 (MQ=255)
           aGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGccc                                                              <  1:3565319/76‑1 (MQ=255)
           aGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGccc                                                              <  2:4409412/76‑1 (MQ=255)
           aGCGGCGATTTAATCGGTGCACAAGGCGTGAGATTGGAAGACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGccc                                                              >  2:920438/1‑76 (MQ=255)
             cggcgATTTAATCGTTGCACAAGGCGTGAGATTGGAACACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGt                                                            <  1:2920326/76‑1 (MQ=255)
              ggcgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTa                                                           <  1:5079818/76‑1 (MQ=255)
              ggcgAGTTAATCGTGGCACAAGGCGGGCGATTGGAATACAATGTCGCGCATGTTGGTTTTATCGGCCTTGCCCGTa                                                           >  1:2626532/1‑76 (MQ=255)
                cgATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaaa                                                         <  2:430762/76‑1 (MQ=255)
                 gATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaaaa                                                        >  2:3359185/1‑76 (MQ=255)
                 gATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaaaa                                                        >  1:5625418/1‑76 (MQ=255)
                 gATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTaaaa                                                        >  1:2501478/1‑76 (MQ=255)
                   tttAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACg                                                      >  2:1891360/1‑76 (MQ=255)
                   tttAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACg                                                      >  2:2422500/1‑76 (MQ=255)
                   tttAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATAGGCCTTGCCCGTAAAACg                                                      >  2:1876620/1‑76 (MQ=255)
                    ttAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGa                                                     <  1:2229128/76‑1 (MQ=255)
                     tAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAt                                                    >  1:1347728/1‑76 (MQ=255)
                      aaTCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAtt                                                   >  1:3375350/1‑76 (MQ=255)
                      aaTCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAtt                                                   >  1:1564355/1‑76 (MQ=255)
                      aaTCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGTCCTTGCCCGTAAAACGAtt                                                   <  1:1619709/76‑1 (MQ=255)
                         cGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAttttt                                                <  1:2275128/76‑1 (MQ=255)
                          gTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAtttttt                                               <  2:5450476/76‑1 (MQ=255)
                          gTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAtttttt                                               <  1:4376940/76‑1 (MQ=255)
                          gTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGAtttttt                                               >  2:1693323/1‑76 (MQ=255)
                           ttGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTta                                              >  2:70119/1‑76 (MQ=255)
                            tGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTtat                                             >  1:4324237/1‑76 (MQ=255)
                            tGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTtat                                             >  2:1998729/1‑76 (MQ=255)
                              cacaAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTatat                                           >  1:2004811/1‑76 (MQ=255)
                               acaAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATc                                          >  1:5636473/1‑76 (MQ=255)
                                caAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCa                                         <  1:3927022/76‑1 (MQ=255)
                                caAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCa                                         >  1:2659294/1‑76 (MQ=255)
                                 aaGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCAc                                        <  1:5740353/76‑1 (MQ=255)
                                  aGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACg                                       >  2:4269492/1‑76 (MQ=255)
                                    gCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTCTATCACggg                                     <  1:4173084/76‑1 (MQ=255)
                                    gCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACggg                                     <  2:3931642/76‑1 (MQ=255)
                                     cGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACgggg                                    >  1:109497/1‑76 (MQ=255)
                                      gTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGa                                   <  1:2233919/76‑1 (MQ=255)
                                      gTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGa                                   <  2:3375350/76‑1 (MQ=255)
                                      gTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGa                                   <  1:920438/76‑1 (MQ=255)
                                      gTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGa                                   >  1:3546/1‑76 (MQ=255)
                                      gTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGa                                   >  1:2803489/1‑76 (MQ=255)
                                        gagaTTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGc                                 <  2:1507147/76‑1 (MQ=255)
                                           aTTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCctc                              <  2:3637908/76‑1 (MQ=255)
                                           aTTGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCctc                              >  1:4319416/1‑76 (MQ=255)
                                            ttGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGCTTTTTTATATCACGGGGAGCctct                             >  2:104939/1‑76 (MQ=255)
                                            ttGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCctct                             <  1:5759162/76‑1 (MQ=255)
                                            ttGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCctct                             >  1:4097857/1‑76 (MQ=255)
                                            ttGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCctct                             <  2:2004811/76‑1 (MQ=255)
                                             tGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCtctc                            <  1:5517456/76‑1 (MQ=255)
                                             tGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCtctc                            >  1:748927/1‑76 (MQ=255)
                                             tGGAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCtctc                            <  1:3030376/76‑1 (MQ=255)
                                              ggAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCTCGGGGAGCCTCTCa                           >  2:2941760/1‑76 (MQ=255)
                                              ggAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCa                           <  2:2569512/76‑1 (MQ=255)
                                               gAATCCAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCag                          <  2:3955169/76‑1 (MQ=255)
                                               gAATACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCag                          >  2:4553670/1‑76 (MQ=255)
                                                aaTACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCaga                         <  2:562965/76‑1 (MQ=255)
                                                aaTACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCaga                         >  1:145982/1‑76 (MQ=255)
                                                aaTACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCaga                         <  1:1693323/76‑1 (MQ=255)
                                                  tACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAgg                       >  2:3521485/1‑76 (MQ=255)
                                                  tACAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAgg                       <  2:1237620/76‑1 (MQ=255)
                                                   aCAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGc                      <  1:1371642/76‑1 (MQ=255)
                                                   aCAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGc                      <  1:3049725/76‑1 (MQ=255)
                                                   aCAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGc                      >  1:5355458/1‑76 (MQ=255)
                                                   aCAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGc                      <  1:1245853/76‑1 (MQ=255)
                                                   aCAATGTCGCGCCTTTTGTTTTTATGCGCCTTGCCCGTAAAACGAGTTTTTATATCACGGGGAGCCTCTCAGAGGc                      >  2:4783130/1‑76 (MQ=255)
                                                    cAATTTCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCg                     >  1:3244516/1‑76 (MQ=255)
                                                        ttCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGctat                 >  1:41337/1‑76 (MQ=255)
                                                         tCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGttatt                <  2:1814723/76‑1 (MQ=255)
                                                         tCGCGCCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGttatt                >  2:2061837/1‑76 (MQ=255)
                                                          tgcgcCTTTTGTTTTTATGGGCCTAGCCCGTAAAACGATGTTTTATATCACGGGGAGCCTCTCAGAGGCGttatta               <  2:79855/75‑1 (MQ=255)
                                                            cgcCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTAcc             <  2:3546/76‑1 (MQ=255)
                                                             gcCTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTAccc            <  2:2501478/76‑1 (MQ=255)
                                                              ccTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCa           >  1:231492/1‑76 (MQ=255)
                                                              ccTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCa           <  1:3569690/76‑1 (MQ=255)
                                                              ccTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCa           <  2:5240915/76‑1 (MQ=255)
                                                              ccTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTAACCa           >  1:4924606/1‑76 (MQ=255)
                                                              ccTTTTGTTTTGATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCa           >  1:1838603/1‑76 (MQ=255)
                                                               cTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCaa          >  1:1060826/1‑76 (MQ=255)
                                                               cTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCaa          <  1:1330312/76‑1 (MQ=255)
                                                               cTTTTGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCaa          >  2:3906247/1‑76 (MQ=255)
                                                               cTTTTGTTTTTATAGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCaa          <  1:1876620/76‑1 (MQ=255)
                                                                ttttGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAAc         >  2:3786137/1‑76 (MQ=255)
                                                                  ttGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACtt       <  1:3906248/76‑1 (MQ=255)
                                                                  ttGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACtt       <  1:4553671/76‑1 (MQ=255)
                                                                   tGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATGTCACGGGGAGTCTCTCAGAGGCGTTATTACCCAACTTg      <  2:4968412/76‑1 (MQ=255)
                                                                   tGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTg      >  2:4236172/1‑76 (MQ=255)
                                                                   tGTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTg      >  2:2386931/1‑76 (MQ=255)
                                                                    gTTTTTATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGa     <  1:3519787/76‑1 (MQ=255)
                                                                     tttttATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAg    >  2:5824077/1‑76 (MQ=255)
                                                                     tttttATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAg    <  1:2422500/76‑1 (MQ=255)
                                                                      ttttATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAgg   <  2:4112312/76‑1 (MQ=255)
                                                                      ttttATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAgg   <  2:5713952/76‑1 (MQ=255)
                                                                       tttATGGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACCTGAGGa  >  1:750094/1‑76 (MQ=255)
                                                                            |                                                                        
GCGATCAAAAAAGCGGCGATTTAATCGTTGCACAAGGCGTGAGATTGGAATACAATTTCGCGCCTTTTGTTTTTATGGCCTTGCCCGTAAAACGATTTTTTATATCACGGGGAGCCTCTCAGAGGCGTTATTACCCAACTTGAGGA  >  NC_000913/4175872‑4176017

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