Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,194,577 2 bp→GC coding (97‑98/774 nt) yjaD → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,194,5770CG100.0% 220.4 / NA 84R33A (CGG→GCG) yjaDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (51/33);  total (51/33)
*NC_0009134,194,5780GC97.6% 205.2 / NA 83R33A (CGG→GCG) yjaDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base C (48/33);  total (50/33)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAACGGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGG  >  NC_000913/4194507‑4194645
                                                                      ||                                                                   
tCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCa                                                                 >  1:909920/1‑76 (MQ=255)
 cACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGGAGCaa                                                                >  1:4126419/1‑76 (MQ=255)
  aCGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCCTATGGGGAAGCGGCaaa                                                               >  2:4517742/1‑76 (MQ=255)
    ggCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAAtt                                                             <  1:567027/76‑1 (MQ=255)
    ggCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAAtt                                                             >  2:3130002/1‑76 (MQ=255)
    ggCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAAtt                                                             >  1:3105351/1‑76 (MQ=255)
     gCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAAttt                                                            >  2:1653316/1‑76 (MQ=255)
      cTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTc                                                           >  1:576626/1‑76 (MQ=255)
      cTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTc                                                           <  1:363746/76‑1 (MQ=255)
       tggtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCg                                                          >  2:1774006/1‑76 (MQ=255)
       tggtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCg                                                          >  1:713128/1‑76 (MQ=255)
       tggtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCg                                                          >  1:4212320/1‑76 (MQ=255)
       tggtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCg                                                          <  2:3657933/76‑1 (MQ=255)
         gtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGAt                                                        >  1:4484411/1‑76 (MQ=255)
         gtggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGAt                                                        <  1:3846063/76‑1 (MQ=255)
          tggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATc                                                       >  2:4050487/1‑76 (MQ=255)
          tggGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATc                                                       >  1:999226/1‑76 (MQ=255)
          tcggTCGTCAGCCCTGAACAAAAATTATGGCTGCCCGGCGGAGAATTGCCATATGCCGAAGGGGCAAAGTTCGATc                                                       >  2:2558345/3‑76 (MQ=255)
           gggTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCGATATGGCGGAGCGGCAAATTTCGATCt                                                      >  1:2324291/1‑76 (MQ=255)
            ggTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCtt                                                     <  1:4301023/76‑1 (MQ=255)
             gtcgtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtg                                                    <  1:3782674/76‑1 (MQ=255)
             gtcgtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtg                                                    >  1:1556988/1‑76 (MQ=255)
              tcgtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTtgt                                                   <  1:497700/76‑1 (MQ=255)
                gtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTgg                                                 <  1:3680477/76‑1 (MQ=255)
                gtcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTgg                                                 <  2:2801320/76‑1 (MQ=255)
                 tcAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTggg                                                >  2:3744513/1‑76 (MQ=255)
                  cAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGt                                               >  1:2261962/1‑76 (MQ=255)
                  cAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGt                                               <  2:2796173/76‑1 (MQ=255)
                   aGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTc                                              <  1:957314/76‑1 (MQ=255)
                   aGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTc                                              >  1:2968030/1‑76 (MQ=255)
                   aGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTc                                              <  2:2433560/76‑1 (MQ=255)
                    gCCATGAACAAAAATTATGGTTGCCGCAGGGAGAATGGCTATATGGAGAAGCGGCAAATTTCGGTCTTGTGCGTCa                                             >  2:4721683/1‑76 (MQ=255)
                    gCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCa                                             <  2:3285373/76‑1 (MQ=255)
                    gCCATGAACAAAAATTATGGCTGCCGAAGGGAGAATTGCCATACGGCGTAGCGGCAGACGTCGATCCTGTGGGACa                                             >  2:4623726/1‑76 (MQ=255)
                     ccATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAg                                            >  1:4076387/1‑76 (MQ=255)
                      cATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgc                                           <  2:3273926/76‑1 (MQ=255)
                      cATGAACAAAAATTATGGGTGCCGAAGGGAGAATTGCCATATGGCGCAGCGGAAAATTTCGATCTTGTCGGTCAgc                                           >  2:3192383/1‑76 (MQ=255)
                        tGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgc                                         >  2:629259/1‑76 (MQ=255)
                        tGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgc                                         <  2:3695545/76‑1 (MQ=255)
                        tGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgc                                         <  1:2527087/76‑1 (MQ=255)
                          aaCAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgcgc                                       >  1:3606532/1‑76 (MQ=255)
                          aaCAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAgcgcgc                                       <  2:263020/76‑1 (MQ=255)
                            cAAAAATTATGGTTGCCGAAGGGAGAATCACCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCAc                                     >  1:3879218/1‑76 (MQ=255)
                             aaaaaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACt                                    >  2:4067652/1‑76 (MQ=255)
                              aaaaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAACTTTCGATCTTGTGGGTCAGCGCGCCCtc                                   >  2:57896/1‑75 (MQ=255)
                              aaaaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTa                                   >  1:9158/1‑76 (MQ=255)
                               aaaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTAc                                  <  2:4558070/76‑1 (MQ=255)
                               aaaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTAc                                  <  1:836585/76‑1 (MQ=255)
                                aaTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACa                                 >  2:3668758/1‑76 (MQ=255)
                                 aTTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACaa                                >  1:2245440/1‑76 (MQ=255)
                                  ttATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACaaa                               >  2:2327567/1‑76 (MQ=255)
                                   tATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAAt                              >  1:3951196/1‑76 (MQ=255)
                                   tATGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAAt                              <  1:2344154/76‑1 (MQ=255)
                                    aTGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATc                             <  1:3793721/76‑1 (MQ=255)
                                     tGGTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAAtag                            >  2:1301176/1‑74 (MQ=255)
                                      ggTTGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAACTTCGATCTTGTGGGTCCGCGCGCACTACAAATCgg                           >  2:992115/1‑76 (MQ=255)
                                         tGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGa                        >  2:752020/1‑76 (MQ=255)
                                         tGCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGa                        >  1:4068503/1‑76 (MQ=255)
                                          gCCGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGaa                       <  2:3791360/76‑1 (MQ=255)
                                           ccGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAAt                      >  2:2313187/1‑76 (MQ=255)
                                           ccGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAAt                      <  1:3888600/76‑1 (MQ=255)
                                            cGAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATg                     >  1:137479/1‑76 (MQ=255)
                                             gAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATgg                    <  2:1655504/76‑1 (MQ=255)
                                             gAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATgg                    >  1:1898082/1‑76 (MQ=255)
                                             gAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATgg                    >  1:2083704/1‑76 (MQ=255)
                                             gAAGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATgg                    >  2:4401738/1‑76 (MQ=255)
                                               aGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCa                  <  2:1184948/76‑1 (MQ=255)
                                               aGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCa                  <  1:763948/76‑1 (MQ=255)
                                               aGGGAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCa                  >  2:2095586/1‑76 (MQ=255)
                                                 ggAGAATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgg                <  1:278915/76‑1 (MQ=255)
                                                 gaagaCTTGCCATATGGCGAAGCGGCACATTAGGAACGTGTAGGTCAGCGCGCACTGCAAATCGGCGCATGGCAgg                <  1:2558345/74‑1 (MQ=255)
                                                  gagaATTGCCATATGGCGAAGCGGCAAATGTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAggg               >  2:1722101/1‑76 (MQ=255)
                                                   agaATTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAgggg              <  2:2485481/76‑1 (MQ=255)
                                                    gaATTGCCATATGGCGAAGCGCCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAggggg             <  2:3254032/76‑1 (MQ=255)
                                                    cacttGCCATACGGCGCAGCGGCAAATTTCGACCTTGTGGCCCAGCGCGCACTACACATCGGCGCATGTCAggggg             >  2:180312/4‑76 (MQ=255)
                                                     aaTTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAACTCGGCGAATGGCAGGGGGa            >  2:4729441/1‑76 (MQ=255)
                                                      aTTGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGaa           <  2:2810557/76‑1 (MQ=255)
                                                        tGCCATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAAcc         <  1:4289328/76‑1 (MQ=255)
                                                          ccATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCCGGGGGAACCTg       >  2:706129/1‑76 (MQ=255)
                                                           cATATGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGt      >  2:2087945/1‑76 (MQ=255)
                                                             tatGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGttt    >  2:4423743/1‑76 (MQ=255)
                                                             tatGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGttt    >  1:1094906/1‑76 (MQ=255)
                                                              atGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTg   >  2:2928045/1‑76 (MQ=255)
                                                               tGGCGAAGCGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTgg  <  1:2744578/76‑1 (MQ=255)
                                                                      ||                                                                   
TCACGGCTGGTGGGTCGTCAGCCATGAACAAAAATTATGGTTGCCGAAGGGAGAATTGCCATATGGCGAACGGGCAAATTTCGATCTTGTGGGTCAGCGCGCACTACAAATCGGCGAATGGCAGGGGGAACCTGTTTGG  >  NC_000913/4194507‑4194645

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