Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,944,676 2 bp→GC intergenic (+20/+75) trpT → / ← yifDA tRNA‑Trp/transcription regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,944,6760CG98.6% 185.4 / ‑3.9 71intergenic (+20/+76)trpT/yifDAtRNA‑Trp/transcription regulator
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (25/45);  total (25/46)
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.80e-01
*NC_0009133,944,6770GC98.6% 181.5 / ‑4.9 71intergenic (+21/+75)trpT/yifDAtRNA‑Trp/transcription regulator
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (25/45);  total (25/46)
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.86e-01

GGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAACGTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCAT  >  NC_000913/3944609‑3944746
                                                                   ||                                                                      
ggTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAt                                                                <  1:741283/76‑1 (MQ=255)
ggTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAt                                                                >  1:780127/1‑76 (MQ=255)
ggTCTCCAAAACCGGGAGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAt                                                                >  2:1794386/1‑76 (MQ=255)
 gTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtt                                                               <  2:2715200/76‑1 (MQ=255)
  tctcCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAttt                                                              <  2:2541967/76‑1 (MQ=255)
  tctcCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCCGCCAGAAATCATCCTTAGCGAAAGCTAAGGActt                                                              >  2:3227477/1‑76 (MQ=255)
   ctcCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttt                                                             <  1:4739569/76‑1 (MQ=255)
   ctcCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttt                                                             >  2:4399363/1‑76 (MQ=255)
     ccAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttttt                                                           >  1:2357259/1‑76 (MQ=255)
      cAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAttttttt                                                          <  1:4740084/76‑1 (MQ=255)
       aaaaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGAtttttttt                                                         <  2:847063/76‑1 (MQ=255)
        aaaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTa                                                        <  1:3233413/76‑1 (MQ=255)
        aaaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAACCATCCTTAGCGAAAGCTAAGGATTTTTTTTa                                                        <  2:4076399/76‑1 (MQ=255)
         aaCCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTAt                                                       >  1:265195/1‑76 (MQ=255)
          acggggTGTTGTGATTTCGAGTCGCTCGCCCCCTGCCAGAAATCATCCTTAGCTAACGCTAATGATTTTTTTtctc                                                      <  2:1313379/73‑1 (MQ=255)
           ccGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCt                                                     <  2:1464688/76‑1 (MQ=255)
             gggTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTCTCTGa                                                   >  2:861941/1‑76 (MQ=255)
             gggTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGa                                                   >  2:955288/1‑76 (MQ=255)
              ggTGTTGGGGGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGaa                                                  <  2:4380402/76‑1 (MQ=255)
              ggTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGaa                                                  <  2:2094430/76‑1 (MQ=255)
              ggTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGaa                                                  <  2:4366122/76‑1 (MQ=255)
               gtgtTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGaaa                                                 <  2:3938300/76‑1 (MQ=255)
               gtgtTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGaaa                                                 >  2:2333925/1‑76 (MQ=255)
                 gtTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGACATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATa                                               <  1:3036523/76‑1 (MQ=255)
                 gtTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTCTCTGAAATa                                               >  1:4789725/1‑76 (MQ=255)
                  ttGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATaa                                              <  2:3939198/76‑1 (MQ=255)
                    gggAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAAcc                                            >  2:1748506/1‑76 (MQ=255)
                     ggAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAAccc                                           <  1:3990733/76‑1 (MQ=255)
                      gAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCct                                          <  1:2673447/76‑1 (MQ=255)
                       aGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCctc                                         <  2:1406732/76‑1 (MQ=255)
                       aGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCctc                                         >  2:4782054/1‑76 (MQ=255)
                        gTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCctct                                        <  2:2056307/76‑1 (MQ=255)
                         tgcGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCtctc                                       <  1:1795322/74‑1 (MQ=255)
                          tCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTcc                                      >  2:4213641/1‑76 (MQ=255)
                          tCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTcc                                      <  2:2215950/76‑1 (MQ=255)
                          tCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTcc                                      <  2:3639685/76‑1 (MQ=255)
                          tCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTcc                                      <  1:244475/76‑1 (MQ=255)
                          tCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTcc                                      <  1:451971/76‑1 (MQ=255)
                           ggAGTCGCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCg                                     <  2:1528386/75‑1 (MQ=255)
                           cGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCg                                     <  2:606002/76‑1 (MQ=255)
                             aGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGaa                                   >  2:1613907/1‑76 (MQ=255)
                             aGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGaa                                   <  1:3476808/76‑1 (MQ=255)
                              ttctctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAg                                  >  2:2858069/2‑76 (MQ=255)
                              gTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAg                                  <  2:1291432/76‑1 (MQ=255)
                              gTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAg                                  <  2:3019197/76‑1 (MQ=255)
                                 tctcCGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTaa                               <  1:275402/76‑1 (MQ=255)
                                   cccGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAAt                             <  2:524560/75‑1 (MQ=255)
                                    ccGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATc                            >  1:2268625/1‑76 (MQ=255)
                                    ccGCCCCTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATc                            <  1:980655/76‑1 (MQ=255)
                                       ccccTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCtt                         >  2:4348947/1‑76 (MQ=255)
                                        cccTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTc                        >  2:1844549/1‑76 (MQ=255)
                                        cccTGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTc                        <  2:4056008/76‑1 (MQ=255)
                                           tGCCAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTAc                     <  2:1182333/76‑1 (MQ=255)
                                             ccAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCg                   >  1:2796889/1‑76 (MQ=255)
                                             ccAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCg                   >  2:4404019/1‑76 (MQ=255)
                                              cAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCgg                  >  2:2701762/1‑76 (MQ=255)
                                              cAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCgg                  <  1:3580248/76‑1 (MQ=255)
                                              cAGAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCgg                  <  1:869652/76‑1 (MQ=255)
                                                gAAATCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCa                <  2:780127/76‑1 (MQ=255)
                                                 aaaTCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCAt               >  2:3615819/1‑76 (MQ=255)
                                                  aaTCATCCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCAGAAGTAAATCCTTCTACCGGCATc              >  1:4797634/1‑76 (MQ=255)
                                                    taatcCTAGGCGAAAGCTAGGGTTTTTTTTAATCGGAAATAACCCTCCCCGAATTAATCCTTCTACCGGCATCCtt           <  2:3794749/74‑1 (MQ=255)
                                                     catcCTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCtt           >  2:2144199/1‑76 (MQ=255)
                                                        ccTTAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgcc        <  1:3554265/76‑1 (MQ=255)
                                                          ttAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccga      <  1:4399363/76‑3 (MQ=255)
                                                          ttAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccag      <  1:2404261/76‑1 (MQ=255)
                                                           tAGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccagc     <  1:2478929/76‑1 (MQ=255)
                                                            aGCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccagcc    >  1:3552124/1‑76 (MQ=255)
                                                             gCGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTgccagcca   <  2:445702/76‑1 (MQ=255)
                                                              cGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCAt  <  2:4423405/76‑1 (MQ=255)
                                                              cGAAAGCTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCAt  <  2:288345/76‑1 (MQ=255)
                                                                   ||                                                                      
GGTCTCCAAAACCGGGTGTTGGGAGTTCGAGTCTCTCCGCCCCTGCCAGAAATCATCCTTAGCGAAACGTAAGGATTTTTTTTATCTGAAATAACCCTCTCCGAAGTAAATCCTTCTACCGGCATCCTTGCCAGCCAT  >  NC_000913/3944609‑3944746

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